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THÉORIE 

DE 

L'ÉCOULEMENT  TOURBILLONNANT  ET  TUMULTUEUX  DES  LIQUIDES 

DANS  LES  LITS  RECTlLIGNES  A  GRANDE  SECTION. 


ÉTUDE  DES   RÉGIMES   GRADUELLEMENT   VARIÉS   (l 


«  1.  La  publication,  par  M.  Bazin,  de  ses  Expériences  sur  la  distribution 
des  vitesses  dans  les  tuyaux  (2),  fournit  en  ce  moment  aux  hvdrauliciens  les 
premières  données  précises,  acquises  à  la  Science,  touchant  l'établissement 
du  régime  uniforme  à  l'entrée  et  dans  la  première  partie  amont  des  lits 
cylindriques  allongés  où  ce  régime  existe.  En  effet,  dans  ces  expériences, 
faites  sur  un  tuyau  en  ciment  lissé  de  8orI"  de  diamètre  et  8o,a  de  longueur, 
le  m  es  u  rage  des  vitesses  des  filets  fluides,  à  travers  les  trois  sections  situées 
au  quart,  au  milieu  et  aux  trois  quarts  de  la  longueur,  a  permis  de  recon- 
naître que  ces  vitesses  u  étaient  pareilles  aux  points  homologues  sur  les 


(')  Ce   Mémoire   fait   suite   à   celui   que  j'ai  publié   l'année  dernière  sous  Je  titre  : 

Théorie  de  l'écoulement  tourbillonnant  et  tumultueux  des  liquides  dans  les  lits 
rectilignes  à  grande  section,  et  qui  était  presque  entièrement  consacré  au  régime 
uniforme. 

(  '2  )  Mémoires  présentés  par  divers  savants  à  V Académie  des  Sciences  de  l  Institut 
de  France,  t.  XXXII,  n°  G;  voir  la  paye  14  du  Mémoire. 
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(6  ) 

deux  dernières  sections,  entre  lesquelles  régnait  par  conséquent  l'unifor- 
mité du  régime,  mais  qu'elles  étaient  notablement  moins  inégales,  du 
centre  au  contour,  sur  la  première.  Dans  celle-ci",  leur  rapport  à  la  vitesse 
moyenne  U,  aux  distances/-  de  l'axe  exprimées  par  les  fractions  suivantes 
du  rayon  R, 

.L  —  n  1  1  1  i  2.  «  1  i;  . 

|>     — u7  8'  8'  8'  8'  8'  8'  8'  16'  '  ' 

avait  respectivement  les  valeurs 

yT=i,i2o5,  1,1190,   i,io6j,   1,0970,   1,082."),   1,0070,    1,0120,  0,929.5,  0,8600,  0,7000, 

.  lundis  que,  |)lus  en  aval,  dans  la  seconde  moitié  du  tuyau,  ce  rapport  était 
(  t.iï  mpyjennfc,  avec  des  écarts  de  sens  divers,  pouvant  al'cr,  entre  les  deux 
'  sections  considérées  de  celte  moitié,  jusqu'à  0,0 1 2  5)  : 

JT  =1,1675,   i,i6o5,   1 , 1 4 7 -> ?   i,i258,   1,0923,   1,0473,   1,0008,  0,9220,  o,8^65,  0,7^15. 

»  L'uniformité  emploie  donc  à  s'établir  une  assez  grande  longueur, 
supérieure,  comme  on  voit,  à  cinquante  fois  le  rayon,  dans  le  cas  d'un 
tuyau  circulaire  à  parois  polies,  puisque  les  vitesses  n'y  étaient  pas  réglées 
après  un  parcours  de  vingl-cinq  diamètres;  et  il  y  a  lieu  d'étudier  le  régime 
graduellement  varié  qui  règne,  sur  celte  longueur,  entre  la  première  section 
amont  où  les  filets  fluides,  avant  terminé  leur  rapide  épanouissement  consé- 
cutif à  la  contraction  de  l'entrée,  sont  désormais  presque  parallèles,  sans 
courbure  sensible,  mais  encore  beaucoup  trop  rapides  près  de  la  paroi, 
beaucoup  trop  lents  sur  l'axe  pour  se  conserver  tels,  et  la  section,  relative- 
ment très  distante,  où  même  le  filet  le  plus  central  a  pris  enfin  toute  sa  vi- 
tesse, après  s'être  accéléré  jusque-là  à  mesure  que  se  ralentissait  le  fluide 
extérieur  plus  retenu  par  la  paroi. 

»  Tel  est  le  genre  de  régime  graduellement  varié,  non  considéré  jusqu'ici 
dans  les  grandes  sections  ou  dans  les  mouvements  tourbillonnants,  qui 
fera  le  principal  objet  de  la  présente  étude.  Toutefois,  je  rattacherai  cette 
étude  à  la  théorie  générale,  que  je  reprendrai  d'abord,  des  autres  modes 
d'écoulement  graduellement  variés,  permanents  ou  non  permanents,  bien 
plus  fréquents  dans  les  cours  d'eau  découverts  :  savoir,  de  ceux  où  il  se  pro- 
duit soit  d'un  point  à  l'autre  de  la  longueur,  soit  sur  place,  des  change- 
ments de  section  fluide  et  de  vitesse  moyenne,  assez  bien  amenés  pour 
laisser  subsister  partout,  avec  quelque  approximation,  la  distribution  des 


(  7  ) 
vitesses  caractéristique  du  régime  uniforme,  altérée  seulement  dans  une 
mesure  comparable  à  ces  changements  eux-mêmes.  Leur  théorie,  très 
facilitée  par  l'emploi  immédiat,  qu'elle  comporte,  de  la  méthode  des 
approximations  successives  à  partir  des  formules  du  régime  uniforme,  a  été 
indiquée,  il  est  vrai,  dès  1 87 1  (')  et  publiée  in  extenso  dans  un  Mémoire  de 
l'année  suivante  1872  (2);  mais  je  lui  ai  trouvé  depuis  des  simplifications 
et  des  compléments  qui  permettent  d'en  abréger  beaucoup  et  d'en  mieux 
synthétiser  l'exposition. 


§  II.  —  Équations  fondamentales  de  l'écoulement  graduellement  varié. 

»  2.  J'ai  déjà  donné  la  définition  et  les  équations  fondamentales  des 
régimes  graduellement  variés,  dans  mon  étude  de  l'année  dernière  sur 
l'écoulement  tourbillonnant  et  tumultueux  (:î  ).  Les  composantes  transver- 
sales v,  w  de  la  vitesse  moyenne  locale  V  y  sont  assez  petites,  ainsi  que  les 
dérivées,  par  rapport  à  l'abscisse  x  ou  au  temps  /,  soit  de  la  composante 
longitudinale  u,  soit  de  la  section  normale  fluide  n,  soit  de  la  vitesse 
moyenne  de  débit  U  à  travers  cette  section,  pour  que  les  produits  de  toutes 
ces  quantités,  entre  elles  ou  par  le  petit  coefficient  z  du  flottement  inté- 
rieur, soient  négligeables  dans  les  équations  du  mouvement.  De  plus,  les 
dérivées  successives  en  x  ou  t  des  mêmes  petites  quantités  sont  d'ordres 
de  grandeur  de  moins  en  moins  sensibles,  en  sorte  que  l'on  doit  n'en  tenir 
compte  qu'à  des  degrés  d'approximation  de  plus  en  plus  élevés;  et  même 
la  dérivée  de  1  en  x,  comparable  au  produit  de  s  par  certaines  des  petites 
dérivées  précédentes,  peut  être  supposée  nulle. 

»  5.  Cela  posé,  les  deux  équations  indéfinies  du  mouvement,  où 
figurent  les  deux  très  petites  accélérations  latérales  /,  w' ',  se  trouveront 


(')   Comptes  rendus,  t.  LXXIII,  p.  3^  et  101  ;  3  et  10  juillet  1871. 

(2)  Essai  sur  la  théorie  des  eaux  courantes,  aux  t.  XXIII  et  XXIV  du  Recueil  des 
savants  étrangers  ou  des  Mémoires  présentés  etc  ,  §§  VI,  IX,  X,  XI,  XII,  XXVI, 
XXXVI.  Noir  surtout  le  §  XL,  rédigé  en  i8j3,  et,  au  Tome  XXIV,  une  Addition,  p.  5o,. 

(8)  Comptes  rendus,  22  et  29  juin  1896;  t.  CXXII,  p.  1 449  el  iJ'7t  ou  n°  20,  p.  17. 
du  précédent  Mémoire  consacré  surtout  au  régime  uniforme  et  intitulé  :  Théorie  de 
l'écoulement  tourbillonnant  et  tumultueux  des  liquides  dans  les  lits  rectilignes  à 
grande  section. 


(  8  ) 
débarrassées,  comme  on  a  vu  par  mon  étude  de  l'année  dernière  (1),  des 
termes  en  s;  et  elles  exprimeront  que  la  pression  moyenne  p  varie,  à  l'in- 
térieur de  chaque  section  a  normale  à  l'axe  des  x,  comme  dans  un  fluide 
sans  frottements,  sa  différentielle  suivant  un  petit  chemin  sjdy'1  H-  riz'1  y 
égalant  sa  valeur  statique  moins  le  produit  ?(v'  dy  -+-  w  riz). 

»  Celui-ci,  par  unité  de  longueur  du  chemin  suivi,  est,  comme  v'  eUr',  du 
second  ordre  de  petitesse;  on  pourra  donc,  à  peu  près  toujours,  le  né- 
gliger. Mais  sa  présence  dans  l'expression  de  rip  entraîne  évidemment  la 
condition  d'intégrabilité 


(0 


dv>        dw' 
dz  dy 


»  Ce  sera  une  des  équations  indéfinies  du  problème;  et  il  faudra,  par 
conséquent,  malgré  l'extrême  petitesse  de  ses  termes,  y  avoir  recours,  sauf 
dans  les  àeux  cas  simples  de  sections  rectangulaires  d'une  grande  largeur 
constante  et  de  sections  circulaires,  où  d'évidentes  considérations  de  symé- 
trie en  tiendront  lieu.  Elle  servira  à  déterminer  les  petites  composantes 
transversales  v,  w  de  la  vitesse,  concurremment  avec  la  condition  de  con- 
tinuité ou  de  conservation  des  volumes  fluides 

.  dv         dw  du 

\     '  dy         dz  dx 

»  i.  L'expression  de  p,  ainsi  formée  pour  l'époque  actuelle  t  et  à 
partir  de  Yaxe  hydraulique  actuel  où  elle  est  appelée  p»,  c'est-à-dire  à 
partir  du  point  où  cet  axe  perce  la  section  g  comprenant  le  point  considéré 
(x,  y,  s),  aura  donc  la  partie  non  hydrostatique  —  f  p  (V  dy  -f-  w'  riz),  duc 
aux  forces  centrifuges  du  fluide  ( — oe',  —  on/  par  unité  de  volume)  au 
moins  dans  la  mesure  où  la  vitesse  V  est  normale  au  plan  des  yz  ou  de  g. 
Mais  cetlc  partie  sera  du  second  ordre  de  petitesse,  et,  dans  la  dérivée 
de  p  en  x  figurant  au  premier  membre  de  l'équation  indéfinie  en  u' ,  elle  ne 
fournira,  par  suite,  que  des  termes  du  troisième  ordre,  négligeables,  même 
quand  on  poussera  jusqu'à  la  seconde  approximation  le  calcul  du  régime 
graduellement  varié.  L'équation  indéfinie  principale  de  ce  régime  sera 
donc  la  relation  (^3)  de  ma  Note  du  29  juin  1896  (2),  ou,  par  suite,  en 
prenant  comme  variables  n,  '(  les  coordonnées  homologues  de  y,  z  dans 


(  ')  Comptes  rendus,  t.  CXXII,  p.  1620,  ou  n°  27,  p.  23,  du  précédent  Mémoire. 
(2)  Comptes  rendus,  t.  CXXII,  p.  1021,  ou  bien  p.  i\  du  Mémoire  précédent. 


(9) 
une  section  semblable  à  <r,  mais  de  rayon  moyen  i,  la  formule  (26)  de  la 
môme  Note,  savoir 


LF^>y+!LF(^>y 


Sit; 


(3)  iK^J+iKo^J+Aïi-Ax. 

»  Dans  cette  formule,  I  désigne  la  pente  motrice  actuelle,  F(r,,  '()  une 
fonction  censée  donnée  pour  chaque  forme  de  la  section  et  u0  la  vitesse  u 
sur  une  génératrice  convenue  du  fond,  par  exemple,  sur  la  génératrice 
équidistante  des  deux  bords.  Il  s'y  joint  la  condition  (26)  spéciale  au  con- 
tour de  <t  et  qui  est 

dl-L 
(4  )  (sur  le  contour)     F(u,  l) -~  =  -  ksjb0f(r„  Ç), 

dv  désignant  toujours  une  petite  normale  menée  au  contour  dans  la  sec- 
tion de  rayon  moyen  1,  à  partir  de  chaque  point  intérieur  (r,,  ()  infiniment 
voisin,  et  /(r,,  '()  une  fonction  donnée  sur  tout  le  contour,  généralement 
peu  différente  de  1  le  long  du  contour  mouillé  y,  mais  nulle  aux  surfaces 
libres. 

»  On  a  vu  d'ailleurs  qu'il  résulte  de  ces  deux  équations,  pour  détermi- 
ner la  vitesse  u0  au  milieu  du  fond,  la  formule  notée  (27),  qui,  en  appe- 
lant Oïl/la  moyenne  des  valeurs  de /('<%  '()  le  long  du  contour  mouillé  y  et, 
finalement,  Dïlu'  la  valeur  moyenne  de  u'  sur  toute  l'aire  n,  est 

»  5.  Enfin,  il  a  été  démontré  au  même  endroit:  d'une  part,  que  ces 
équations  sont  indépendantes  du  choix  des  axes  rectangles  des  y  et  des  z 
dans  un  plan  parallèle  à  la  section  n  •  d'autre  part,  qu'elles  déterminent 
complètement  les  vitesses  u  à  travers  celle-ci,  dès  qu'on  y  donne  la  pente 
motrice  actuelle  I  et,  en  chaque  point  {y,  z)  ou  (n,  £),  l'accélération  lon- 
gitudinale actuelle  u'.  Or  cette  dernière  peut  s'évaluer,  sauf  erreur  de 
l'ordre  des  carrés  ou  produits  qu'on  néglige,  en  remplaçant  la  vitesse  lon- 
gitudinale u  par  la  vitesse  totale  V  =  \ju2  +  v2  -+■  w'1 ,  ou  en  prenant  ainsi 
X accélération  tangentielle,  dérivée  complète  par  rapport  au  temps  (c'est- 
à-dire  pour  une  même  particule  suivie  dans  son  mouvement  moyen  local) 

de  la  vitesse  V;  et  l'on  voit  alors  qu'elle  est  tout  aussi  indépendante  des 
B.  2 


(  io) 

axes  coordonnés  que  la  pente  motrice  elle-même, 


d  (o     ,     Po 


(6)  l  =  ~  ds  \~°    '     ?g 

où  60  est  Yallilude  et  p0  la  pression  moyenne  sur  Yaxe  hydraulique  actuel, 
dont  on  considère  l'élément  ds  compris  entre  la  section  g  proposée,  qui 
lui  a  été  menée  normale,  et  une  section  voisine  parallèle. 

»  En  faisant  idéalement  la  même  substitution  permise  de  V  et  V0  à  a 
et  u0  dans  les  premiers  membres  des  formules  (3),  (4)  et  (5),  ces  équa- 
tions auront,  dans  toutes  leurs  pari  ies,  une  signification  géométrique  ou 
physique,  indépendante  des  axes  généraux  de  coordonnées  que  l'on  pourrait 
choisir  ensuite  pour  étudier  dans  chaque  cas  le  mouvement  à  toutes  les 
époques;  caries  y,  z,  our,,£,  qui  y  figureront,  seront  des  coordonnées 
locales  propres  à  la  section  fluide  g  menée,  à  l'époque  /,  normalement  en 
un  point  donné  quelconque  de  Yaxe  hydraulique  tel  qu'il  est  à  ce  moment, 
c'est-à-dire  de  l'axe  du  tuyau  quand  il  s'agit  d'un  tuyau  plein,  et  d'une 
coupe  longitudinale  de  la  surface  libre  actuelle  quand  il  s'agit  d'un  cours  d'eau 
découvert. 

»  Enfin,  les  vitesses  et  les  accélérations,  variant  lentement  d'un  point  à 
l'autre  dans  les  sens  parallèles  ou  presque  parallèles  à  la  direction  moyenne 
de  l'écoulement,  ont  sensiblement  les  mêmes  valeurs  et  les  mêmes  dérivées 
en  sens  divers  aux  points  correspondants  de  sections  voisines  g,  presque 
pareilles  d'ailleurs,  menées  en  un  même  endroit  dans  le  fluide  suivant  des 
orientations  un  peu  différentes,  mais  toutes  à  peu  près  normales  aux  vi- 
tesses V.  Et  l'on  peut  même  y  confondre  des  dérivées  de  u  suivant  des 
sens  transversaux  (des y  ou  des  z)  un  peu  différents,  faisant  entre  eux  des 
angles  du  premier  ordre  de  petitesse,  pourvu  que  ce  soit  dans  des  termes, 
provenant  des  frottements,  où  figurait  d'abord  le  petit  facteur  s,  comme 
sont  ceux  en  F  (y;,  "Ç)  de  (2)  et  de  (3). 

»  Donc,  si  l'on  prend,  pour  toutes  les  époques,  un  axe  fixe  des  x  sui- 
vant la  direction  moyenne  de  l'écoulement,  et  des  axes  des  y  et  des  z  per- 
pendiculaires, les  quantités  u,  ç,  w,  u',  v',  w'  pourront,  également,  se 
rapporter  à  ces  axes  fixes  et  communs  dans  les  formules  (1),  (2),  (3),  (4) 
et  (5),  où  la  pente  motrice  I  sera  la  dérivée  en  x,  changée  de  signe,  de 

l'expression  60  -h  ^,  mesurée  le  long  d'un  axe  hydraulique  généralement 

variable  d'un  instant  à  l'autre,  mais  dont  les  éléments  ds  ne  feront  avec 
l'axe  des  x  que  des  angles  ayant  leurs  carrés  négligeables. 


(  "  ) 


§  III.  —  Équations  qui  déterminent  le  mode  de  distribution  des  vitesses 

dans  l'écoulement  varié. 

»  6.  Si,  après  avoir  divisé  (3)  et  (4)  par  k\j'B0,  l'on  élimine  I  de  (3) 
par  la  relation  (5),  il  vient,  pour  déterminer  le  mode  de  distribution  des 
vitesses  à  travers  la   section  t,   le  système  suivant  d'équations,  où  nous 

avons  pris  comme  inconnue  la  même  fonction  — —  ( i  )  que  dans  le 

cas  du  régime  uniforme  (')  : 


O. 


dll  —  Uo  r  , .1.     „    i       <V  _1\    U  —  Ut 


(au  contour)  F-j =  = — /,  (au  milieu  du  fond) 


»   Posons 

(8)  iif-|f  =  F'^-)  +  ^.7F!(-,,r). 

«oAry'Bo  6Do«0   /. 

F,  étant  la  fonction,  déjà  considérée,  qui  définit  pour  chaque  forme  de 
section  la  distribution  des  vitesses  dans  le  régime  uniforme  et  que  déter- 
mine le  système 


^9)      1  d¥ 

[  (au  contour)  F  —p  =  —  f,         (au  milieu  du  fond)  F,  =  o. 

»   Les  formules  (7)  deviendront  évidemment,  en  F2, 

3      i  d¥9 

I    (au  contour)  F  -—  =  o,  (au  milieu  du  fond)  F2  =  o  ; 


et  il  est  clair  que  celles-ci,  une  fois  l'accélération  u'  connue  aux  divers 
points  (y,  z)ou  (y],  Z,)àea,  détermineront  complètement  la  fonction  F2,  qui 


(!)  Comptes  rendus,  6  juillet  1896,  t.  CXXIII,  p.  7,  ou  §  VIII  du  Mémoire  précé- 
dent, p.  26. 


(  >2  ) 

dépendra  non  seulement  de  r,,  £,  mais  aussi  des  autres  quantités  entrant 
dans  u' . 

»  Quand  la  fonction  F2  sera  ainsi  trouvée,  la  relation  (8)  donnera,  pour 
exprimer  dans  toutes  les  sections  semblables  le  rapport  de  la  vitesse  u  en 
un  point  quelconque  à  la  vitesse  u0  au  milieu  du  fond,  la  formule 


§  IV.  —  Relation  entre  la  vitesse  moyenne  et  la  vitesse  au  fond. 

»  7.  Le  rapport  de  la  vitesse  moyenne  U  à  la  vitesse  w0  au  milieu 
du  fond  s'obtiendra,  par  suite,  en  prenant  la  moyenne  des  valeurs  du 
second  membre  de  (i  i)  sur  toute  l'aire  g  de  la  section  fluide;  ce  qui  donne, 
comme  généralisation  de  notre  formule  (3i)  de  régime  uniforme  ('),  si 
DlcF2  désigne  la  valeur  moyenne  de  la  fonction  F2  (r,,  £), 

(  1 2)  £  =  i  +  (fcm-FO  v'b;  4-  — 4=-  £  ^F2. 

»  On  voit  qu'il  suffirait  de  connaître  0ïLF2  pour  pouvoir  tirer  de  (12)  la 
vitesse  u0,  au  milieu  du  fond,  en  fonction  de  la  vitesse  moyenne  ou  de 
débit  U;  après  quoi,  la  substitution  de  cette  valeur  de  u0  dans  l'équa- 
tion ('5)  donnerait,  entre  la  vitesse  moyenne,  le  rayon  moyen,  la  pente 
motrice  et  l'accélération  moyenne  DKu' ,  une  relation,  propre  à  jouer  dans 
les  régimes  graduellement  variés  le  rôle  capital    de    l'équation    usuelle 

-  I  =  b\J2  dans  le  régime   uniforme.   Or  l'expression  désirée  de  0ïlF2  se 

déduit  aisément  des  équations  (9)  et  (10)  définissant  F,  et  F2,  sans  qu'on 
ait,  à  beaucoup  près,  besoin  de  les  intégrer. 

•»  8.  Ajoutons,  en  effet,  les  premières  équations  (9)  et  (10),  respecti- 
vement multipliées  parF2^et  par  —  ¥{dn;  et  observons  que  da,  ou  dydz, 
est  le  produit  du  carré  du  rayon  moyen  par  l'élément  d'aire  d-rid'C,  dans  la 
section  semblable  de  rayon  moyen   r.  Puis  intégrons  les  résultats  dans 


(')   Comptes  rendus,  6  juillet  1896,  t.  CXXIII,  p.  8,  ou  p.  27  du  Mémoire  précé- 
dent. 


(  '3  ) 
toute  l'étendue  de  celle-ci,  après  avoir  remplacé  les  différences 

Par    ^   F2.F-^--F1.F^-    ,-   F2.F-7r--F1.F-7F   •    Les    termes    ou 


ûrr,  \    2  '       û?7)  '  '       <r/r,  J     dt\    -'       dX,  '  "       <iw 

paraissent  ces  différences  se  transformeront,  à  la  manière  ordinaire,  en 
intégrales  de  contour,  que  les  secondes  relations  (9),  (10)  simplifieront  et 
réduiront  à  la  partie  mouillée  du  contour.  Revenons  enfin  à  la  section 
effective  c  et  à  son  contour  mouillé  y,  en  multipliant  les  différentielles 
sous  les  signes/  par  les  facteurs  convenables.  Alors,  si  on  (F,  m')  désigne 
la  valeur  moyenne  du  produit  F,  (r,,  'Ç)id  dans  tonte  l'étendue  a,  nous 
aurons,  après  avoir  divisé  par  a, 

(i3)  —   / /F2  -^  -f-  0H/.3TLFj  =  ,m  F4.3]W  —  Dft  (F,  a'). 

Jy  7. 

»  Dans  cette  relation,  le  premier  terme  égale  évidemment  le  produit 
de  DKf  par  une  valeur  de  F2  intermédiaire  entre  la  plus  petite  et  la  plus 
grande  que  prenne  cette  fonction  le  long  du  contour  mouillé  /.  Or,  si  les 
vitesses,  aux  divers  points  de  la  paroi,  étaient  réparties  dans  le  mouve- 
ment varié  comme  dans  le  mouvement  uniforme,  on  y  aurait,  d'après  (8), 
F2=o,  le  premier  membre  de  (8)  s'y  réduisant  à  F, (?;,£).  Sans  avoir 
besoin  d'admettre  qu'il  en  soit  rigoureusement  ainsi,  il  est  clair,  par  ana- 
logie avec  ce  qui  a  lieu  dans  le  régime  uniforme,  que  les  écarts  relatifs  de 
vitesse,  propres  au  mouvement  varié,  seront  bien  moindres  le  long  du 
contour  mouillé  que  dans  tout  l'intérieur  de  la  section.  Autrement  dit,  la 
fonction  F2  se  maintiendra,  le  long  de  y,  beaucoup  plus  voisine  que  dans 
l'aire  a  de  sa  valeur  zéro  réalisée  au  milieu  du  fond.  Donc  le  premier 
terme  de  (i3)  est  négligeable  devant  le  deuxième,  et  cette  relation  donne 

(14)  =DnFiXlu<-XliFiu>), 

ÔWf 

»   9.  Telle  est  l'expression  de  D1oF2  à  substituer  dans  (12).  Remplaçons 

ensuite  le  binôme  1 -h  ^v^^ô^F,  parsa  valeur,  1  /    °'.  L' ,  contenant  le  coef- 
ficient usuel  b  qui  entre  dans  la  formule  du  régime  uniforme  (  '  ),  et  mettons 


(')   Comptes  rendus,   6  juillet  1896,   t.  CXXIII,   p.  8,   formules  (32)  et   (33);  ou 
p.  27  du  Mémoire  précédent. 


(  >4) 

d'ailleurs  cette  valeur  en  facteur  commun  au  second  membre.  Nous  aurons 
.    „.  U  /B^ïÇ/l      ,    ,  A  /~r  a  DïlFi3ïlu'-DK(Flu')~] 

(I3)  ^  =  \~îT~V^I{\^~fyr    -JB^ÏÔT/-    -J' 

expression  où  le  terme  qui  suit  l'unité  dans  la  parenthèse  sera  très  petit 
et  aura  son  carré  négligeable,  puisque  nous  admettons  ici  un  mode  de 
distribution  des  vitesses  voisin  de  celui  du  régime  uniforme.  Comme  on 
veut  avoir  u'I  en  fonction  de  U2,  il  reste  à  renverser  cette  valeur  du  quo- 
tient de  U  par  u0  et  à  l'élever  au  carré,  en  employant  d'ailleurs  la  formule 
du  binôme  et  en  substituant  à  B0u20DTU/,  dans  le  terme  en  u,  sa  valeur  de 
première  approximation  6U2.  Il  vient,  pour  représenter  (au  facteur  près  pg) 
le  frottement  extérieur  moyen  B0u20DTl/  par  unité  d'aire,  la  formule 

(ib)  B0^/=&U2-2*y/^- ' * '-■ 

»   On  donne  au  second  membre  une  signification  plus  intuitive  en  obser- 
vant que  F,   et  At/^-. reviennent,  d'après  (8)  et  (i5),   dans  les  petits 

/  r  .      U  —  Mo         ?  «ovHn  J       u    "  —  ll0 

termes  ou  u  est  en  facteur,  a -=L>  A— ,T— >  et  ont  pour  produit  — r 

Alors  '/0  s'élimine  du  dernier  terme  de(i6);  et,  en  indiquant  finalement  par 
(u2)'  la  dérivée  complète  de  u2  relative  au  temps,  c'est-à-dire  sa  dérivée 
liai'  prise  en  suivant  une  même  particule,  par  3ïl(m2)'  la  valeur  moyenne 
de  (u~)' ,  il  vient 

(17)  B.ll^Il/^  OU2 -*---[  — ^  -2DIL//J. 

Telle  est  la  valeur  de  B0 wj 3Tb/ qu'il  faudra  porter  dans  l'équation  (5),  où 
figure  la  pente  motrice  I.  Résolue  par  rapport  à  I,  cette  équation  sera 

(.8)  l  =  6U'ï  +  ipi^  -au/|. 


§  V.  —  Formule  générale  pour  la  valeur  moyenne,  sur  une  section, 
de  toute  dérivée  complète  par  rapport  au  temps. 

»  10.  Il  manque  encore  à  son  dernier  terme  et  à  celui  de  (17)  d'être 
reliés  le  plus  simplement  possible  à  la  vitesse  moyenne  U  ou  à  ses  dérivées 
partielles  en  x  et  t. 


(  '5  ) 
»  Pour  y  parvenir,  démontrons  d'abord  la  formule  générale  suivante, 
où  il  s'agit  de  tout  courant  fluide,  permanent  ou  non  permanent,  dont  les 
particules,  d'une  densité  p  constante  ou  variable,  possèdent  des  vitesses  V 
ayant,  à  l'époque  /,  les  composantes  u,  v,  w  suivant  des  oc,  y,  z  fixes  et 
où,  U  étant  la  vitesse  moyenne  (de  débit)  suivant  l'axe  des  x,  à  travers 
les  sections  a  normales  à  cet  axe  et  fonctions,  comme  elle-même,  de  x  et 
de  t,  t  désigne  d'ailleurs  toute  fonction  continue  des  quatre  variables  /,  x, 
y,  z,  enfin,  t',  sa  dérivée  complète  par  rapport  au  temps,  ou  dérivée  prise 
en  suivant  durant  l'instant  dt  la  particule  venue  en  (;r,y,z)à  l'époque/  : 

(19)  f^'dn  =  ±  fj.dc  +  ±  fjUTdc. 

»   Considérons,  en  effet,  à  l'époque  t,  la  somme  2/tit',  pour  toutes   les 

particules  m  =  zdndx  comprises  dans  le  volume  fluide  /  dx  i  dn,  que 

limitent  les  sections  t,  n  ayant  deux  abscisses  voisines  et  constantes  x, 
x  -+-  Ix;  et  soient  t,  la  valeur,  à  l'époque  t  -h  dl,  de  la  fonction  t  pour  la 
particule  m,  ns  la  section  fluide  correspondant,  pour  la  même  époque  /  4-  dt, 
à  chaque  abscisse  intermédiaire  entre  x  et  x  -+-  Ax.  Il  est  clair  que 

(20)  2mT'  =  irS~m~>  ~  Iwt)' 

et,  d'autre  part,  que  Imz' ,  Imz  ont  les  deux  expressions  respectives 

dx  (  z-'  dn     et       /  dx  f  z-.dn. 

»  Pour  évaluer  Im-r,,  observons  que  la  masse  1m  comprend,  à 
l'époque  t  -f-  dl,  la  tranche  fluide  limitée  par  les  deux  sections  g, 
d'abcisses  x,x  -+-  Ix,  où  les  particules  m  se  grouperont  en  éléments  de 
volume  dxdnK  donnant  les  éléments  d'intégrale  dxz-dn{,  moins  le  fluide, 
p(udt)da  à  fort  peu  près,  entré  par  chaque  élément  de  la  première  sec- 
tion durant  l'instant  dt,  et  donnant  l'élément  d'intégrale  —  dtzu-dn,  plus 
enfin  le  fluide  analogue  o(adt)da  sorti  dans  le  même  instant  par  chaque 
élément  de  la  dernière  section  1'  et  fournissant  à  l'intégrale  l'élément 
dtzu-.dn.  La  somme  lm-7t  sera  donc 

»   Divisons  par  dt  son  excédent  sur  l'expression  de  1m-,  et  nous  aurons 


(  '6) 
évidemment,  d'après  (20), 

C   d*ffd°  =  l      dx{iJr<h  +  £fj"T'h} 

ce  qui,  en  supposant  àx  infiniment  petit,  revient  bien  à  la  formule  (19). 


£  VI.       Applications  de  cette  formule,  notamment  à  l'équation  de  continuité 
du  fluide  pour  toute  l'étendue  des  sections,  etc. 

»  11.  Une  première  application,  indispensable,  de  (19)  s'obtient  en  po- 
sant t  =  1,  de  manière  à  exprimer  la  conservation  de  la  masse  fluide  ïm. 
Dans  le  cas  auquel  nous  nous  bornons  d'un  liquide,  il  vient  ainsi,  après 
suppression  du  facteur  alors  constant  p,  et  en  observant  que  U  est  la 
valeur  moyenne  de  u  sur  toute  l'aire  a,  l'équation  de  continuité  en  CJ  et  g, 
savoir 

(2I)  dt  +  -ïïr=0' 

»  Elle  signifie  que  la  tranche  fluide,  gAx,  comprise  entre  deux  sections 
voisines  a,  g\  à  abscisses  constantes,  et  à  travers  lesquelles  les  vitesses 
moyennes  sont  U   et  U',   décroît,    par  unité   de   temps,   d'une   quantité 

ch 

-  \x,  égale  à    la  différence  de   leurs  débits,   qui    est   t'U' — n\J,  ou 


dt 

d. -A 


!±x. 


dx 
»   Posons  maintenant,  dans  (19),  t  =  soit  u,  soit  m2;  et  faisons  d'ailleurs 

où  r,  désigne  ainsi  (  ')  l'excédent  sur  l'unité,  toujours  positif,  du  rapport  du 
carré  moyen  des  vitesses  u  à  travers  une  section  au  carré  de  leur  moyenne  U, 
et  où  a,  peu  différent,  comme  on  sait,  de  1  -+-3t;  (2),  est,  suivant  l'usage  des 


(  '  )  Il  est  vrai  que  r,  désigne  déjà  le  rapport,  dans  chaque  section  a,  de  la  coordon- 
née y  à  une  certaine  ligne  de  la  section.  Mais  il  ne  résultera  de  ce  double  emploi  de  r, 
aucune  confusion,  les  deux  rôles  n'ayant  rien  de  commun. 

(2)  En  effet,  si  l'on  pose  77  — i  -+-  A,  la  valeur  moyenne  de  A  sur  toute  l'aire  a  de  la 


(  >7) 

hydrauliciens,  le  rapport  analogue  du  cube  moyen  des  vitesses  u  au  cube 
de  la  vitesse  moyenne.  La  formule  (19),  divisée  par  p<j,  donnera 

(23)  sw('=-|-ar+ ^ j,        sM8.y=_^_i_v_+_a5_j 

Dédoublons  les  termes  où  figure  une  dérivée  en  x,  en  y  considérant 
(1  -+-  Y))U2ff,  <xU3c  comme  produits  de  Uc  par  (1  H-  tq)TJ  ou  par  <xU2;  puis 
éliminons,  grâce  à  (21),  la  dérivée  de  Ut  en  a?.  Il  vient,  après  quelques 
réductions  évidentes, 

/  , dV  rf.(i4-r,)L  U  d<s 

\  dt  dx  '  1   dt 

(2.4)  / 

,    ,. -, ,       a. (ih-t))U2        TTc/.aL2        /  v  U2  rf<r 

=»-("•)'=      dt      ■tD-as--<«-,-*>-rar 

»  Enfin  ces  valeurs,  portées  dans  (1  8)  et  (17),  donneront  aisément  les 
expressions  désirées  de  la  pente  motrice  et  du  frottement  extérieur  moyen 
par  unité  d'aire  (au  facteur  près  og)  : 


(25) 


(26) 


t       7  rT9  y       /  \   d  II2 

I=;/Jlj-  -4  -{-(2  7.  —  I  —  7}  )  -5 

7         v  y  rAr  2^ 

1  +2  7]  c/L  a  —  1  —  2  7)  U  ofcr     _     U  /'.,  rf.a  —  T)  rfr 


t/f  "•  7  <://,  sr  \  dx  dt 


! : 1 U  ; -+-    -y1 

g  /  dt  y  g  \  dx  dt 


section  sera  nulle.  Donc  les  carré  et  cube  moyens  1  4-7),  a  de  =-tj  savoir 

;))I(l  +  2A  +  A2)      et      31L(i4-  3A+3A2-+- A3), 

seront  respectivement 

i  +  01c(A2)      et      1  -+-33)t(A2)4-;)rc(A3). 

Ainsi,  d'une  part,  tj  =  31L(A2),  et  les  valeurs  absolues  de  A  sont,  en  moyenne >  de 
l'ordre  de  \frt;  d'autre  part,  Oïl  (A3)  étant  dès  lors  comparable  à  r.y/r,,  c'est-à-dire,  en 
général,   beaucoup  plus   petit   que  3  011  (A2)  ou   3?,,  %   se  réduit  bien,  sensiblement, 


(  i8  ) 


§  vil.  —  Équation  générale  du  mouvement  et  formule  du  frottement  extérieur, 

à  une  première  approximation. 

»  12.  Les  coefficients  2x  —  i  —  y),  14- 2Y),  a  —  i  —  2y),  2(01  —  i— yi),2yj, 
i  -h  3  Y]  —  a,  calculés  par  les  relations  (22)  qui  définissent  r,  et  a,  pourront 
être  réduits  à  leurs  valeurs  sensiblement  constantes  de  régime  uniforme, 
dans  tous  les  écoulements  assez  graduellement  variés  pour  que  le  mode 
de  distribution  des  vitesses  diffère  peu  de  ce  qu'il  est  dans  ce  régime;  car 
ces  coefficients  multiplient  des  dérivées  de  U  ou  de  a  petites  du  premier 
ordre,  et  les  parties  de  7),  et  ajoutées  par  de  pareilles  variations  de  régime 
n'apporteraient  aux  termes  considérés  que  des  corrections  non  linéaires, 
supposées  négligeables. 

»  De  plus,  dans  les  écoulements  graduellement  variés  auxquels  nous 
voulons  nous  borner  d'abord,  et  où  se  feront  assez  lentement  les  change- 
ments de/orme  de  a  influant  sur  les  valeurs  de  régime  uniforme  de  y]  et  a, 
les  petites  parties  variables  de  ces  coefficients  seront,  comme  celles  mômes 
que  contiendra  le  rapport  de  u  à  U  et  d'où  elles  proviendront,  de  l'ordre  des 
dérivées  partielles  premières  de  U  ou  de  a;  et  leurs  dérivées  en  x  ou  en  l 
atteindront,  par  suite,  comme  les  dérivées  secondes  de  U  ou  de  a,  le 
deuxième  ordre  de  petitesse.  C'est  dire  qu'à  une  première  approximation,  le 
dernier  terme,  double,  de  chacune  des  équations  (25)  et  (26)  sera  négligeable. 

»  Les  deux  équations  (25)  et  (26)  auront  ainsi  leurs  seconds  membres 
réduits  aux  quatre  premiers  termes;  et  les  quatrièmes,  affectés  des  coeffi- 
cients a  —  1  —  2ri,  1  -+-  3-n  —  a,  très  petits  par  rapport  aux  coefficients  pré- 
cédents, seront  môme  peu  sensibles.  On  pourra  dire,  en  particulier,  que  la 
pente  motrice  I  se  divise  en  trois  parties  principales,  employées  respectivement , 

l'une,    b[j2-,   à  vaincre    le  frottement  extérieur  de   régime   uniforme;   la 

deuxième,  (20c  —  1  —  ^)-f-(— ~)>  a  accélérer  d'amont  en  aval  le  mouvement, 
en  y  accroissant  la  hauteur  due  à  la  vitesse  moyenne  U;  enfin,  la  troisième, 
— ■ — -  —7-1  à  accélérer  le  mouvement  sur  place. 

g  Cil 

»  D'après  la  formule  (2G),  le  frottement  extérieur  moyen  par  unité 
d'aire  comprend  pareillement  trois  parties  principales,  dont  les  deux  der- 
nières, dépendant  des  mômes  variations  du  mouvement,  montrent  que,  à 
égalité  de  vitesse  moyenne  U,  la  vitesse  au  fond  u0  croit  quand  le  mouvement 


(  i9  ) 
s'accélère  ainsi  soit  d'amont  en  aval,  soit  sur  place.  Ces  accélérations  ten- 
dent donc  à  égaliser  les  vitesses  à  travers  chaque  section,  comme  Dupuit 
en  avait  fait  quelque  part  la  remarque  pour  le  cas  du  régime  permanent,  il 
y  a  un  demi-siècle,  dans  ses  Études  théoriques  et  pratiques  sur  le  mouvement 
des  eaux  courantes  ( 1  ). 

»  15.  L'hypothèse,  faite  ici,  d'un  mode  de  distribution  des  vitesses 
peu  différent  de  celui  du  régime  uniforme,  astreint  évidemment,  dans  la 
formule  (26),  les  termes  qui  suivent  b\J2  à  être  notablement  moindres  que 
bU2;  sans  quoi  le  rapport  de  u0  à  U  en  serait  trop  altéré.  Mais,  heureusement, 
les  coefficients  2(a  —  1  —  •/]),  2r„  1  -h  3r,  —  a  de  ces  termes  sont  de  petites 
fractions  des  coefficients  correspondants  iv.  —  1  —  yi,  i  -4-  iv\,  ce  —  1  —  ir{ 
dans  la  formule  (25);  car  r,  ne  dépasse  guère  0,02  ou  o,o3  (sauf  dans  le 
cas  de  parois  très  rugueuses)  et  x  égale  1  4-3y),  à  un  écart  près  de  l'ordre 
den\/n.  Aussi  le  terme  bU2  pourra-t-il,  dans  (26),  être,  comme  on  l'admet, 
très  supérieur  à  ceux  qui  le  suivent,  sans  que,  dans  (23),  les  termes  corres- 
pondant à  ceux-ci,  ou  dus  à  la  variation  du  mouvement,  soient  tenus  d'être 
moindres  que  le  terme  en  h.  Autrement  dit,  grâce  aux  inégalités  modérées 
des  vitesses  à  travers  chaque  section  dans  les  écoulements  tourbûlonnanls ,  le 
régime  peut  y  être  graduellement  varié,  avec  vitesses  distribuées  approximati- 
vement comme  dans  le  régime  uniforme,  tout  en  différant  beaucoup  d'un 
régime  uniforme.  De  là,  le  champ  étendu  d'application  et  l'utilité  de  l'équa- 
tion (23)  du  mouvement. 


(*)  11  observe  que  le  principe  de  Daniel  Bernoulli,  appliqué,  entre  deux,  mêmes 
sections  normales  d'un  courant  permanent  sensiblement  rectiligne,  aux.  divers  filets 
lluides  de  ce  courant,  impose  à  tous,  de  la  section  amont  à  la  section  aval,  la  même 

variation  de  la  hauteur  —  ou  —  due  à  la  vitesse.  Si  donc,  dit-il  à  peu  près,  les  fi- 

lets  s'accélèrent,  les  carrés  u-  conserveront  entre  eux,  tout  en  grandissant,  leurs  diffé- 
rences primitives  (du  moins  le  long  de  parcours  assez  modérés  pour  qu'on  puisse  y 
négliger  l'action  des  frottements),  et  les  différences  entre  leurs  racines  u  s'atténueront. 
C'est  bien  ce  qu'ont  vérifié,  toujours  pour  l'état  permanent,  les  remous  d'abaisse- 
ment observés  par  M.  Bazin.  Les  valeurs  de  r,  qu'il  y  a  constatées  sont  un  peu  moindres 
que  celles  de  régime  uniforme  {Recherches  hydrauliques;  ire  Partie,  p.  262).  Or  r, 
ou  OÏL  (A2)  mesure,  en  quelque  sorte,  par  sa  racine  carrée,  l'inégalité  relative  moyenne 
de  vitesse  des  filets  fluides  composant  le  courant;  et  la  diminution  de  tj  indique  bien 
une  tendance  de  ceux-ci  vers  l'égalisation  de  leurs  rapidités,  aux  endroits  dont  il 
s'asrit  où  le  fluide  s'accélère. 


(    20) 

»  Au  contraire,  dans  les  mouvements  bien  continus,  où  s'annule  la  vitesse 
u0  à  la  paroi,  yj,  a  —  i  sont  comparables  à  l'unité  ;  et  tout  régime  graduelle- 
ment varié,  avec  vitesses  distribuées  approximativement  comme  dans  un  régime 
uniforme,  est  quasi  uniforme,  quant  à  la  relation  existant  entre  la  pente 
motrice,  la  vitesse  moyenne  et  le  rayon  moyen.  L'étude  d'un  tel  régime 
n'offre  donc  qu'un  intérêt  restreint  (*  ). 


(')  J'en  ai,  toutefois,  donné  la  théorie  au  §  I  d'un  Mémoire  publié  en  1878  au  Jour- 
nal de  Lioiwille. 

Un  cas  particulier,  celui  qui  concerne  l'établissement  du  régime  uniforme  dans  la 
première  partie  amont  d'un  long  tube,  et  que  j'étudie  pour  l'écoulement  tourbillon- 
nant, c'est-à-dire  pour  les  tuyaux,  aux  §§  XXI  à  XXX  du  présent  Mémoire,  offre  un 
certain  intérêt  aux  physiciens;  car  il  permet,  par  le  calcul  de  la  fraction  de  la  charge 
(ou  hauteur  motrice)  qui  est  employée  à  établir  le  régime  uniforme,  d'expliquer  les 
débits  obtenus  par  Poiseuille  dans  celles  de  ses  expériences  sur  l'écoulement  le  long 
des  tubes  fins  qui  constituent  sa  seconde  série,  où  la  longueur  était  insuffisante  pour 
que  la  fraction  de  charge  dont  il  s'agit  ici  pût  être  négligée.  Une  étude  surtout  expé- 
rimentale de  M.  Maurice  Couette  (Thèse  de  doctorat  es  Sciences  physiques,  Sur  le 
frottement  des  liquides,  Paris,  1890),  ayant  appelé  mon  attention  sur  ce  cas  spécial, 
je  l'ai  abordé  dans  les  Comptes  rendus  de  l'Académie  des  Sciences,  d'abord  les  9  et 
16  juin  1890  (t.  CX,  p.  1160  et  1238),  puis,  plus  complètement  pour  les  tubes  circu- 
laires, les  6  et  i3  juillet  1891  (t.  CXIII,  p.  9  et  49)-  Un  premier  aperçu  y  démontre 
que  la  hauteur  de  charge  employée  à  établir  le  régime  uniforme,  dans  un  tel  tube  à 

entrée  évasée,  est  environ  — >  U  désignant  la  vitesse  moyenne.  Après  quoi,  la  mise  en 

compte,  par  des  intégrations  en  série  compliquées,  des  variations  de  vitesse  des  filets 

après  l'entrée,  m'a  donné   1,12  —  >  avec   une  longueur   minima  0,26  — R2U   environ 

A'  £ 

(où  £  =  0,0000001 3 p^  pour  l'eau  à  io°),  sans  laquelle  le  régime  uniforme  n'existerait 
pas,  même  à  l'extrémité  aval,  à  des  écarts  relatifs  près  de  0,01  sur  les  vitesses. 

M.  Jules  Delemer  a  repris  à  fond  les  mêmes  calculs  dans  sa  Thèse  de  doctorat 
es  Sciences  mathématiques,  Sur  le  mouvement  varié  de  l'eau,  etc.  (Paris,  1896),  et 
en  a  fait  l'application  détaillée   à  la   seconde  série   d'expériences  de  Poiseuille.  Il  a 

trouvé  notamment,  pour  la  charge  employée  à  établir  le  régime  uniforme,  i,i23 — > 

o 
quand  on  se  borne,  comme   j'avais  fait,  à  la  première  solution  simple,  ou  au  terme 

U2 

fondamental  de  1  intégrale  générale,  mais  i,i346  —  >  quand  on  prend  les  deux  pre- 

o 

miôres  solutions  simples. 

On  verra  plus  loin,  dans  deux  Notes  annexées  aux  §§  XXV  et  XXVIII  du  présent 
Mémoire,  comment  la  théorie  de  ces  phénomènes  pourrait  se  déduire,  comme  cas 
limite,  de  celle  des  écoulements  tourbillonnants,  convenablement  modifiée. 


(    21     ) 


§  VIII.  —  Usages  de  l'équation  générale  du  mouvement  graduellement  varié. 

»  14.  La  pente  motrice  lest  liée,  dans  un  tuyau  plein,  à  la  dérivée  en  x 
de  la  pression  p0  le  long  de  l'axe  :  dans  un  canal  découvert,  elle  l'est  à 
la  dérivée  analogue  de  la  section  fluide  <y,  car  elle  se  confond  alors  avec  la 
pente  de  superficie,  dont  l'excédent  sur  la  pente  donnée  du  lit  est  le  quo- 
tient, par  dx,  de  la  hauteur  de  la  bande  supérieure,  — dç,  qui  manque 
actuellement  à  la  section  fluide  d'abscisse  x  -+-  dx  pour  égaler  la  section  a 
d'abscisse  x.  En  joignant  la  condition  (21)  de  conservation  des  volumes 
fluides  à  la  relation  (25),  on  aura  donc  deux  équations  aux  dérivées  par- 
tielles, en  U  et/?0  dans  le  cas  du  tuyau  où  la  section  a  est  connue,  en  U  et  n 
dans  le  cas  du  canal  découvert  où  p0,  pression  à  la  surface  libre,  est  ou 
constante,  ou,  du  moins,  donnée.  S'il  s'agissait  d'un  tuyau  élastique,  à  sec- 
tion g  variable  avec/?0,  cas  dans  lequel  ces  deux  équations  (21),  (25)  con- 
tiendraient les  trois  fonctions  inconnues  distinctes  U,  a,  p0  de  x  et  de  l,  la 
théorie  de  l'élasticité  fournirait,  entre  n  etp0,  la  troisième  équation  néces- 
sa  i  re . 

»  On  aura  donc  les  formules  indispensables  pour  rattacher,  dans  le  cas 
général  d'un  régime  non  permanent,  les  états  successifs  du  courant  fluide 
à  son  état  initial,  et  pour  déterminer,  dans  le  cas  plus  particulier  d'un  ré- 
gime permanent,  où  les  équations  deviendront  simplement  différentielles 
en  x,  les  variations,  d'amont  en  aval  ou  d'aval  en  amont,  soit  de  la  pres- 
sion p0  sur  l'axe,  soit  de  la  section  fluide  t. 

»  15.  Ces  équations  expliquent  facilement,  comme  on  peut  le  voir  dans 
mon  Mémoire  cité  de  1872  ('),  les  principales  circonstances  qu'offrent  les 
cours  d'eau  découverts,  soit  parvenus  à  un  état  sensiblement  permanent  et 
étudiés  dans  leurs  longues  parties  à  variations  graduelles  de  chaque  sec- 
tion aux  suivantes,  soit  considérés  dans  des  états  de  crue  ou  de  décrue 
survenant,  les  uns,  assez  vite,  les  autres,  avec  une  certaine  lenteur.  Mais 
l'observation  de  tels  phénomènes  ne  comporte  guère  le  degré  de  précision 
qu'il  faudrait  pour  contrôler,  dans  l'équalion  (25),  les  coefficients 
2x  —  1  —  Y],   1  -f-  2 Y],   a,  —  1  —  2T|  des  termes  dus  à  la  non-uniformité  du 

(')■  Essai  sur  la  théorie  des  eaux  courantes,  §§  XIII  à  XVI,  XXVII,  et  XXXVI  à 
XXXIX. 


(    ") 

régime,  en  tant  qu'ils  diffèrent,  le  premier,  de  l'unité,  qu'il  excède  sensi- 
blement de  5y),  le  deuxième,  aussi  de  l'unité,  et  le  troisième,  de  zéro,  qu'il 
surpasse  de  n  environ.  Or  il  n'en  est  pas  tout  à  fait  de  même  de  la  propa- 
gation de  Yonde  ou  intumescence  que  produit  une  variation  assez  rapide, 
mais  momentanée,  de  la  hauteur  d'eau  et  de  la  vitesse  moyenne,  à  une  des 
deux  extrémités  du  cours  d'eau,  onde  descendante,  ou  dirigée  suivant  le 
courant,  quand  elle  survient  à  Ventrée  ou  extrémité  amont,  onde  ascen- 
dante, allant  contre  le  courant,  quand  elle  survient  à  X embouchure  ou 
extrémité  aval,  et  que  la  vitesse  U  du  courant  est  insuffisante  pour  arrêter 
sa  progression  vers  l'amont.  Dans  ces  deux  cas,  la  vitesse  de  propagation, 
calculable  par  les  équations  (21)  et  (^5),  se  trouve  dépendre  assez  des 
petites  parties  des  coefficients  en  question,  pour  que  sa  mesure  effective  les 
mette  en  évidence,  du  moins  dans  les  cours  d'eau  torrentueux  et  à  fond 
non  poli  où  la  rapidité  um  des  filets  superficiels  excède  fortement  la  vitesse 
moyenne  U.  Il  y  a  donc  lieu  d'évaluer  ici  cette  vitesse  de  propagation,  en 
vue  de  vérifier  expérimentalement  l'équation  générale  (  25)  des  écoulements 
graduellement  variés. 


§  IX.  —  Son  emploi  dans  le  calcul  de  la  propagation  des  ondes  ou  des  remous 

le  long  d'un  courant. 

»  16.  Supposons  qu'il  s'agisse  d'un  canal  rectangulaire  de  pente  con- 
stante et  d'une  très  grande  largeur  également  constante,  où  se  trouve 
établi,  avant  la  perturbation  constituant  Fonde  étudiée,  un  régime  uni- 
forme, à  vitesse  moyenne  U  pour  la  profondeur  d'eau,  également  donnée,  H. 
Soient  U  -h  U'  et  II  -4-  h  les  nouvelles  valeurs  de  ces  quantités  lors  du  pas- 
sage de  l'onde,  ou  U'  et  h  les  petites  variations,  fonctions  de  x  et  de  /, 
qu'ont  subies  celles-ci  à  partir  des  valeurs  primitives  constantes  U,  H. 

»  La  pente  initiale  de  superficie,  donnée  par  la  formule  du  régime  uni- 
forme ou  par  l'équation  (23)  réduite  à  ses  deux  premiers  termes,  est  le 
quotient  de  bXJ2  par  H,  vu  la  valeur  prirpitive  H  du  rayon  moyen.  A  l'état 
non  permanent,  chaque  section  se  relevant  fie  h,  cette  pente  I  devient  évi- 
demment b-jj-  —  y-;  et,   en  même  temps,  le  second   membre   de   (25), 

réduit,  comme  il  a   été   dit,   par  la   suppression   de   son    dernier  terme 
(double),  prend  l'expression 


(U  +  U')«  xU-t-U'rfU'       1 

11  +  /'  v  '       s        dx 


2r.  rfU'        a-  1  —  27)  U  +  U'gf/i 


dx  2         dt  g  II  -+-  h  dt 


(  »3) 

»  A  une  première  approximation,  l'on  pourra  négliger  les  termes  non 
linéaires  en  U'  ou  h  et  même,  à  cause  de  la  petitesse  du  coefficient  b,  les 
produits  de  b  par  U'  ou  par  h.  L'équation  (25)  deviendra  donc 


(27)        — 


rf* 


dx 


dh 
dx 


i  —  2  r,  U  <iA 
"7  H  ~dt 


=  o. 


»  D'autre  part,  l'équation  (21) donnera,  au  même  degré  d'approximation, 
(28) 


dh 
~dt 


dW  dh_ 

dx  dx 


O. 


»  17.  On  peut  tirer  de  (28)  la  valeur  de  la  dérivée  de  U'  en  x,  pour 
la  substituer  dans  la  relation  (27),  différentiée  préalablement  par  rapport 
à  x.  Il  vient  ainsi  l'équation  en  h, 


,         x       drh 
(29)       -^F+2 


I  + 


3( 


l)   5  7] 


2  (14-  2T,) 


L 


d2h  gtt  —  (2a  —  1  —  r()  U2  rf* h 

dx  dt  1  +  2  r,  Gfo;2 


=  O. 


»  Décomposons  son  premier  membre  en  deux  facteurs  symboliques  du 
premier  degré,  ou,  autrement  dit,  adoptant  comme  inconnue  provisoire  la 
fonction  auxiliaire 


(3o) 


dh 

dt 


dh 
dx 


mettons  l'équation  (29)  sous  la  forme 


(3i) 


d± 
dt 


,  d'I 
o)  -~  =  o. 
dx 


»   Les  deux  constantes  o,  o/  auront  respectivement  pour  somme  et  pour 
produit  les  coefficients  des  second  et  troisième  termes  de  (29);  d'où 


(32) 


I  + 


3(a  — 

I  )  5  T, 

2(. 

+  27,) 

3  (a  — 

I  )   57, 

2(l   +2  7,) 

R  désignant  le  nombre  positif  dont  le  carré  est 

"3  (a  —  1)  —  57, 


±  ki/551, 

V     1  +  27, 
V     1+27) 


u 

Uq 


(33) 


R5 


I   +      O.  —   I   —   2T, 


U2 


1  -h  2  r,  \  g-  H 


»   Nous  prendrons  le  radical,  dans  (32),  avec  les  signes  supérieurs  s'il 
s'agit  d'ondes  descendantes,  avec  les  signes  inférieurs  s'il  s'agit  d'ondes 


(  *4  ) 

ascendantes.  Dans  le  premier  cas,  en  plaçant  l'origine  des  x  à  l'entrée  du 
canal,  on  n'aura  le  long  de  celui-ci  que  des  abscisses  x  positives;  et  l'on 
pourra  convenir  de  compter  le  temps  à  partir  d'un  moment  où  l'expres- 
sion (3o)  de  <j/  sera  encore  nulle  pour  x^>o.  Dans  le  second  cas,  en  pla- 
çant l'origine  des  x  à  l'embouchure,  l'on  n'aura,  au  contraire,  que  des 
abscisses  x  négatives;  et  l'on  comptera  de  môme  le  temps  t  à  partir  d'une 
époque  où  -|  =  o  tout  le  long  du  canal,  savoir  pour  x  <  o. 

»  18.  Cela  posé,  formons  les  intégrales  du  problème  pour  la  région  du 
cours  d'eau  située  en  avant  du  plan  x  —  o//,  c'est-à-dire  ayant  des  abscisses 
plus  grandes  que  co'/,  dans  le  cas  d'ondes  descendantes  où  l'on  a  to'<[co, 
mais  plus  petites  que  <m'  t  dans  le  cas  d'ondes  ascendantes,  ou  0/  est  ^>  co. 
Cette  région  comprendra  évidemment  tout  le  canal  si  o  et  u>'  ont  signes 
contraires  ;  ce  qui  arrivera  dans  les  cours  d'eau  franchement  tranquilles  (non 
torrentiels),  où  \j gM  excède  notablement  U.  Même  dans  un  cours  d'eau 
torrentiel,  elle  sera  assez  étendue  pour  que  la  différence  x  —  tôt,  inférieure 
ou  supérieure  à  x  —  ia' t  suivant  que  les  ondes  descendent  ou  remontont 
le  courant,  y  varie  de  part  et  d'autre  de  zéro,  et  finalement  dans  d'aussi 
larges  limites  que  l'on  voudra,  une  fois  /  devenu  assez  grand. 

»  La  région  considérée  donnant,  suivant  les  cas,  x  —  u  l  ^>  o  ou 
x  —  o//<^o,  l'intégrale  de  (3i),  savoir  6=  une  fonction  de  x  —  o//,  y 
devient  simplement  ^  =  o  à  raison  de  la  donnée  d'état  initial,  qui  annule  h 
et  <\>,  quand  t  s'annule,  pour  toutes  les  valeurs  respectivement  positives  ou 
négatives  de  x.  Dès  lors,  l'équation  (3o)  a  elle  môme  pour  intégrale 
h  =  F(x  —  mt),  où  la  fonction  F,  nulle  (au  moins  sensiblement)  pour  les 
valeurs  soit  positives,  soit  négatives  de  sa  variable,  en  vertu  de  la  môme 
condition  d'état  initial,  reste  arbitraire  pour  les  valeurs  respectivement 
négatives  ou  positives  de  cette  variable  x  —  ut. 

»   L'équation  (28),  ainsi  devenue  -7-  [HLT —  (co  —  U)A]  =  o,  montre 

que  l'expression  II U'  —  (co  —  V)h  est  nulle  partout,  comme  aux  points  du 
canal  que  l'onde  n'a  pas  encore  atteints;  et  il  vient,  en  définitive,  pour  la 
solution  cherchée,  vérifiant  d'ailleurs  (27)  non  moins  que  (28), 

(34)  h  =  F(x  —  co/),  U'=^^lh. 


(  «5) 


§  X.  —  Calcul  théorique  de  l'influence  qu'a  la  déformation  de  la  masse  fluide 
sur  leur  vitesse  de  propagation. 


(38) 


»  19.  D'après  ces  formules,  la  constante  a>,  définie  par  (3s)  et  (33),  est 
la  célérité  (ou  vitesse  de  propagation)  des  ondes.  Or,  en  négligeant  des 
termes  très  petits  de  l'ordre  de  r,-,  les  relations  (33)  et  (32)  donnent 


(35) 
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»  D'autre  part,  les  formules  (3^)  de  mon  Étude  de  l'année  dernière  ('  ) 
donnent  de  leur  coté,  pour  exprimer  la  distribution  des  vitesses  u  aux 
diverses  profondeurs  relatives  'Ç,  dans  notre  cours  d'eau  très  large, 


(36) 


u 

û  =  l 


2  Y     \  3 


t  -+- 


Ç-i  Ki-3?). 


En  élevant  le  troisième  membre  soit  au  carré,  soit  au  cube,  puis  multipliant 
par  r/'C,  intégrant  de  £  =  o  à  £  =  i,  et  retranchant  enfin  l'unité,  nous  aurons 


(37) 


4  /  Um 

5V  U 


—  i 


r  =  3r, 


7 


Substituons  dans  le  dernier  membre  de  (35)  cette  valeur  de  %  —  i,  puis 
celle  der,,  et  il  viendra,  comme  expression  générale  de  la  célérité  oj  des 
petites  ondes  de  translation,  le  long  d'un  courant  large  à  régime  uniforme, 


=  u±Noiï 


rjU 


*sjg\l 


u 


7  V  \fgH 


(«I 


U)2 


T 


AA±-L 

35  V  V>H 


("„t-li)3 
U2 


»  20.  Celte  expression  dépend  des  trois  quantités  U,  ±y/^H,  um — U, 
qui  définissent,  l'une,  U,  la  vitesse  moyenne  du  cours  d'eau,  la  seconde, 
±  y^H,  la  célérité  des  ondes  telle  qu'elle  serait  par  rapport  à  la  masse 


(*)   Comptes  rendus,  t.  CXXIII,  p.  12  (ier  juillet  1896),  ou  p.  33  du  -Mémoire  pré- 
cédent. 

B.  \ 


(  26) 

fluide,  si  celle-ci  était  tout  entière  animée  de  la  vitesse  moyenne  U;  enfin, 
la  troisième,  um  —  U,  l'excédent  de  la  vitesse  maxima  um  à  la  surface  sur  la 
vitesse  moyenne  U,  excédent  caractérisant  l'inégalité  absolue  de  vitesse 
des  fdets  fluides.  La  différence  co  —  U  =j=  \jgll,  exprimée  en  majeure  partie 
par  le  troisième  terme  des  second  et  dernier  membres  de  (38),  mesure 
donc  l'influence,  sur  la  vitesse  de  propagation  w,  de  la  rapidité  avec  la- 
quelle se  déforme  sans  cesse  le  milieu  transmettant  les  ondes.  Ce  terme, 
vu  la  petitesse  ordinaire  de  r,,  est  peu  sensible,  comparativement  à  la  somme 
U  ±  yg'H  des  deux  précédents,  sauf  dans  les  deux  cas  :  i°  d'ondes  descen- 
dantes, le  long  d'un  cours  d'eau,  à  fond  rugueux  et  très  rapide,  rendant  à 
la  fois,  dans  ce  terme,  les  fonctions  yi,  U  relativement  considérables,  et 

notable  aussi  la  différence  — ■=•  —  v ^     (précédée  alors  du  siçne  -h)  qui 

2V4'H  u     Vl  b  /  1 

croit  avec  le  rapport  de  U  à  \jg\l\  2°  d'ondes  ascendantes  le  long  d'un 

courant  soit  presque  torrentiel,  soit  faiblement  torrentiel,  où  la  somme 

algébrique  U  —  y^H  s'approche  de  zéro,  par  suite  de  la  quasi-égalité  de 

U  et  de  vgH  ;  ce  qui  accroît  l'importance  du  terme  en  rt  considéré. 


§  XI.  —  Constatation  effective  de  cette  influence  et  vérification  précise 
de  l'équation  du  mouvement. 

»  21.  M.  Bazin  a  effectivement  constaté,  plusieurs  années  avant  l'exis- 
tence de  la  précédente  théorie,  que  la  célérité  w  était  bien,  dans  les  cas 
ordinaires,  U  ±  \/gli,  sauf  toutefois  pour  les  ondes  ascendantes,  le  long 
d'un  cours  d'eau  presque  torrentiel,  où  elle  se  trouvait  sensiblement  ré- 
duite en  valeur  absolue,  c'est-à-dire  plus  grande  que  U  —  s/gÛ  dans  le  sens 
d'amont  en  aval,  conformément  à  la  remarque  précédente  (').  Et  il  a, 
postérieurement  à  l'impression  de  ses  Recherches  hydrauliques,  après  avoir 
pris  connaissance  de  ma  formule  générale  (32),  publié  trois  observations 
faites  en  même  temps  que  les  précédentes,  mais  qu'il  avait  réservées  faute 
de  pouvoir  les  expliquer  par  la  formule  \J±\/gti  (2).  C'est  qu'elles  ren- 
traient justement  dans  le  premier  cas  exceptionnel,  signalé  ci-dessus,  où 


(')  Deuxième  Partie  (relative  aux  remous  et  à  la  propagation  des  ondes)  de  ses 
Recherches  hydrauliques,  p.  36  et  41. 

(2)  Comptes  rendus,  i5  juin  1 885,  t.  C,  p.  i49?-- 


(  27  ) 

les  deux  rapports  de  um  —  U  à  U  et  de  U  à  \/g'H  sont  relativement  considé- 
rables. On  y  avait,  en  effet,  comme  valeurs  observées  respectives  de  H,  U, 
um,  (o  (avec  des  erreurs  possibles  de  3  pour  ioo  environ  sur  co)  : 


H 

=  0,1 IO, 

o,i5o, 

o,235, 

U 

=  3,785, 

2,744, 

3,48i, 

Um 

=  5,5 1, 

3,49' 

4,55, 

<s> 

=  6,25, 

4,32, 

5,75. 

»  Or,  avec  ces  données  (et  g- =  9™, 809),  il  vient  seulement,  pour 
U  -\-\JgH,  les  valeurs  très  insuffisantes  : 

U+y/gH  =  4,824,         3,937,         4,999' 

tandis  qu'on  trouve,  par  les  formules  combinées  (37),  (33)  et  (32),  les 
valeurs  théoriques,  sensiblement  exactes,  comme  on  voit, 

u  =  6,190,         4,3ii,         5,555. 

»   La  formule  approchée  (38)  donnerait  les  résultats  plus  forts  : 
w  =  6,5ir,  4,327,  5,589. 

»  L'erreur,  insignifiante  sur  les  deux  derniers,  atteint  un  vingtième  sur 
le  premier;  c'est  que,  dans  l'expérience  correspondante,  le  nombre  d'or- 
dinaire très  petit  ri  avait  pris  la  valeur  énorme  0,166.  Ce  nombre,  encore 
fort  grand  dans  les  deux  autres  expériences,  s'y  réduisait  cependant,  respec- 
tivement, à  0,059  et  à  0,075. 


§  XII.  —  Calcul  de  l'accélération  longitudinale  u'  dans  un  écoulement 

graduellement  varié. 

»  22.  Notre  première  approximation  donne  la  même  célérité  10  à  toutes  les 
parties  d'une  intumescence  quelconque.  Donc  le  problème  important  de 
leur  lente  déformation  requiert  une  approximation  plus  élevée;  et  celle-ci 
exige  généralement  l'évaluation,  dans  (25)  (p.  17),  des  termes  de  r,  et 
x  —  ri  qui  dépendent  de  la  variation  du   mouvement.   Alors  l'intégration 


(  »8) 

du  système  (10)  et,  par  le  fait  même,  le  calcul  des  accélérations  u!  devien- 
nent inévitables. 

»  Nous  poserons,  pour  définir  le  mode  de  distribution  des  vitesses  aux 
divers  points  (y,  s),  ou  (r,,  '(),  de  la  section  a  d'abscisse  x, 

(3g)  TT  =  o  H-  ni,  ou         u  —  U(<p  -h  ct), 

o  étant  la  fonction  de  r,,  £  et  de  la  forme  de  g,  qui  exprimerait  le  rapport 
de  u  à  U  si  le  régime  uniforme  existait  à  la  traversée  de  cette  section, 
et  tô,  fonction  à  déterminer  de  r,,  £,  ,r  et  /,  représentant  le  petit  écart  dû  à 
la  variation  du  mouvement.  Si  nous  désignons  par  le  symbole  dc  la  diffé- 
rentielle complète  de  la  quantité  écrite  à  la  suite,  c'est-à-dire  sa  différen- 
tielle prise  en  suivant,  durant  l'instant  dl,  une  même  particule  fluide 
(dans  son  mouvement  moyen  local),  nous  aurons  successivement,  vu  que 
u  =  L  (o  -r-  ct),  que  y,  z  ne  figurent  pas  dans  U,  enfin  que  les  termes  non 
linéaires  par  rapport  à  e,  w,  ct  et  aux  dérivées  de  u,  U,  o  en  x  et  /  sont 
négligeables, 

'  ^(9  +  ^7/T  +  L dl 

(dU  dU\        rT  (do  do  do  do\         ,,  (dm  dm\ 

=  !™     +  U  f  o2  +  u  ($  +  u9É  +  ^  +  ,v£)  +  U  (£  H-  U?^ 

rf£    i  <7.r    '  \r/£  '  rt.r  a/  c/c/  \fl7  '  doc 

»  A  une  première  approximation  du  calcul  de  u',  on  pourra  même,  ?c/, 
supprimer  le  dernier  terme  double,  où  figurent  les  deux  dérivées  de  u  en  t 
etx.  En  effet,  dans  le  régime  graduellement  varié  que  nous  considérons 
actuellement,  la  perturbation,  n,  apportée  parla  variation  de  l'écoulement 
au  mode  de  distribution  des  vitesses,  est  censée  liée  aux  changements  de  a 
et  de  U  avec  x  et  /,  au  point  de  n'être  pas  d'un  ordre  de  grandeur  plus 
élevé  que  les  dérivées  premières  de  îouUenict/;  ce  qui  réduit  les  dé- 
rivées analogues  de  ts  à  l'ordre  de  petitesse  supérieur  des  dérivées  secondes 
de  cr  ou  de  U.  Ainsi  l'on  aura 

//    \  /        dU  tt  dU     .,        Tt(do        TT     do  do  do\ 


(  *9> 


§  XIII.  —  Formules  régissant  les  petites  composantes  transversales  de  la  vitesse. 

)>  23.  Le  calcul  de  u!  exigera  donc  l'emploi  non  seulement  de  la  fonction 
connue  <p,  mais  aussi  des  deux  petites  composantes  transversales  v,  w  de  la 
vitesse;  et  il  faut  préalablement  déterminer  celles-ci. 

»  Nous  avons  donné,  à  cet  effet  (p.  8),  la  condition  d'intégrabilité(i) 
et  l'équation  de  continuité  (2).  Il  faudra  y  joindre  la  relation  exprimant 
que  les  particules  situées  à  la  surface-limite  du  fluide  et  supposées  s'y 
mouvoir  avec  les  vitesses  moyen  nés  locales,  y  seront  encore  à  l'époque  t^c-dt; 
et  ce  sera  même  cette  relation  qui,  étudiée  la  première,  nous  suggérera  la 
marche  à  suivre  pour  traiter  la  question. 

»  Cette  partie  du  problème,  savoir,  la  recherche  de  v  et  w,  étant  très 
distincte  de  la  précédente,  nous  ne  nous  y  astreindrons  pas  à  considérer  des 
sections  1  dépendant  d'un  seul  paramètre  fonction  de  x  et  t,  dans  le  genre 
du  rayon  moven.  Nous  y  admettrons  deux  tels  paramètres  distincts  a,  h; 
ou,  autrement  dit,  nous  prendrons,  comme  équation  de  la  surface  limite 
du  fluide,  une  relation  de  la  forme 

(42)  t(^,Ç)  =  o,  avec  *  =  ^-°>  "  =  ^^> 


a,  h  étant  deux  fonctions  lentement  variables,  mais  quelconques,  dexelt, 
et  v0,  z0  les  coordonnées,  par  rapport  aux  axes  fixes  desy  et  z,  du  point  où 
l'axe  hydraulique  perce  la  section  a  d'abscisse  x,  coordonnées  que  nous 
ne  pourrons  pas  toujours  supposer  nulles  et  qui  seront,  en  général,  comme 
a  et  h,  des  fonctions  lentement  variables  de  x  et  t.  Ainsi,  les  coordonnées 
transversales  relatives  r,,  '(  seront  les  quotients  de  y  —  y0,  z  —  z0  par  deux 
longueurs  distinctes,  a,  h,  réductibles,  il  est  vrai,  au  rayon  moyen  dans  le 
cas  de  sections  toutes  semblables.  Néanmoins,  quand  il  s'agira,  même  dans 
ce  cas,  de  sections  rectangulaires  très  larges,  nous  appellerons  a  la  demi- 
largeur  qui,  alors,  deviendra  non  pas  égale,  mais  seulement  proportion- 
nelle au  rayon  moyen  h, 

»  Enfin,  nous  admettrons  que  la  fonction  <p  donnée,  exprimant  le  mode 
de  distribution  des  vitesses  dans  le  régime  uniforme,  soit  de  la  forme 
simple  o  (r,,  '();  ce  que  nous  savons  être  vrai  tout  au  moins  dans  le  cas  de 
sections  semblables,  avec  coefficient  de  frottement  extérieur  B  pareil  aux 
points  homologues,  et  aussi  dans  le  cas  de  sections  rectangulaires  très 
larges,  où  o  dépend  seulement  de  *£. 


(  3o  ) 

»  24.  Exprimons,  conformément  à  la  condition  énoncée,  que  <\>  ne  cesse 
pas  de  s'annuler  quand,  à  partir  de  valeurs  actuelles  de  /,  x,  y,  z  donnant 
•li  =  o,  on  fait  croître  t  de  dt  et  x,  y,  z  deudt,  vdt,  wdl,  ou,  très  sensible- 
ment, de  Uydt,  vdl,  wdt.  Il  viendra 

//o\  db       TT     <id/  d'I  d<\> 

(43)  ^  +  u9^.  +  ^+».^  =  o. 


»   Or,  les  formules  (42)  donnent  immédiatement 

,  .  ,\         d'I  d^\r\       da  i       dy0  d<\i 

'  **'    d(t,x)  ~~  ~~  dV,  L«  rf(^,.r)  +  â  dJT~xT) J  ~  dZ 

et  la  relation  (4^)  peut  s'écrire 


Ç        r/A  i       c/c, 

-h 


A  é/(  £,  .x)         h  d(t,  x) 


1    ûty 


(^£+U<p^-^^  +  Uo^] 


o. 


/ /  r \        I    u  ui,  ^  \  fl?£  '    tfta*  /  \dt  '  cfo y 

'43'      J  i   d'iV  fdz0       ,TT    ^-0\        ,,/ctà        TT     rf/A~| 

»  Cette  condition  au  contour  de  g  prendra  sa  forme  la  plus  simple,  si 
l'on  adopte  comme  fonctions  à  déterminer,  au  lieu  de  v  et  w,  les  deux 
expressions  entre  crochets.  Ou  mieux  encore,  nous  extrairons  de  ces  deux 

parenthèses  une   partie  y  avant  respectivement  la   forme  (a,  h) t.        '     , 

avec  deux  fonctions,  l'une,  /.,  de  x  et  t  seuls,  l'autre,  y,  de  -t]  et  (  seuls,  à 
déterminer  en  vue  des  plus  grandes  réductions  ultérieures  possibles.  En 
un  mot,  introduisant  deux  nouvelles  inconnues  auxiliaires  Val,  Uhpix  la 
place  de  v,  \v,  nous  poserons 

/  dvo        TT     dy0  (  da        TT     da\  c/y        tt 

j  dz0        TT     dz0         y  ( dh        T-     dh\  7  d-;        rT , 

«   La  relation  définie  (4->)  deviendra 
//    \  (dît  d{        db  d-(\        TT/rf«K         ^    \ 

»  2o.  Choisissons,  comme  condition  au  contour  destinée  à  déterminer 
la  fonction  jusqu'ici  disponible  y,  celle  qui  s'obtient  en  annulant  le  coef- 
ficient de  /.  dans  (47).  A  cet  effet,  considérant  y  dans  la  section  type 
(a  =  1,  k  =  1),  où  7),  s  seraient  les  vraies  coordonnées  absolues  y  —  y{>, 
z  —  z0,  supposons  qu'on  mène  sur  chaque  élément  du  contour  y',  à  partir 
d'un  point  intérieur  infiniment  voisin,  une  normale  dv.  Celle  qui  aboutira 


(  3i  ) 
ainsi  au  point  (r,,Ç)  du  contour  aura  ses  deux  projections  dr,,  d'C,  sur  les 
axes  proportionnelles  aux  deux  dérivées  de  6  en  r,,  Ç  figurant  dans  (47); 
en  sorte  que  l'annulation  du  coefficient  de  x  revient  à  poser 

(48)  -f-  =  o        (le  long  du  contour). 


Donc,  la  condition  (4?)  sera,  en  ).  et  {/., 


(49)  -r-  ^  -+-  -jt  [J-  =  o         (le  long  du  contour). 


»   26.  Transportons  actuellement  les  valeurs  (46)  de  v,  w  dans  l'équation 
de  continuité  (2),  où  le  second  membre  aura,  vu  l'expression  U9  -+-  Ugj 

de  u,   une  première  partie, j— '- ,  que  nous  transposerons  au  premier 

membre.  Si  nous  observons  que  la  dérivée  de  o  en  x  a  sa  formule  pareille 
à  (44)»  et>  de  plus,  que  l'aire  a  des  sections  est  proportionnelle  au  pro- 
duit ah,  d'une  part,  les  dérivées  de  o  s'élimineront  du  premier  membre, 
ainsi  que  a,  h,  et,  d'autre  part,  les  termes  en  9  s'y  réduiront  immédiatement, 

en  vertu  de  (21)  [p.  16],  à  —  5  -r-  En  appelant  enfin  A2y  la  somme  des 

deux  dérivées  secondes  directes  de  y  en  r,  et  en  '(,  il  viendra 


,  r    \  »  l   d~  /  \         tt  ( d\  d\x\  d.XJm  _. 

(00)         xa2Ï--;s(9-i)  +  U^  +  ^=--^-=-U 


dm 
dx 


»  Achevons  de  déterminer  les  fonctions  x  et  y  en  annulant  la  somme  des 
deux  premiers  termes.  Il  nous  faut,  pour  cela,  rappeler  l'expression  de 
o  —  1,  résultant  des  formules  (29)  et  (3o)  de  mon  Étude  de  l'année  der- 
nière ('), 

expression  qui  convient  tout  au  moins  pour  les  sections  semblables  et  pour 
celles  qui  sont  rectangulaires  très  larges.  Nous  devrons  donc  poser 

(52)  x  =  AV-Bo"         -  %  A,y  =  F, (il,  C)  -  .111F,. 

»   On  sait  que  cette  dernière  équation  indéfinie,  jointe  à  la  condition  au 


(  ')   Comptes  rendus,  t.  CXXIII,  p.  8,  ou  p.  26  du  Mémoire  précédent. 


(  32) 
contour  (48),  déterminera  complètement  en  n  et  £  la  fonction  auxiliaire  y, 
à  une  constante  près  qui  disparaît  des  dérivées  de  y  en  r,,  'C  et,  par  suite, 
des  expressions  (46)  de  t>,  w. 

»   Enfin,  l'équation  (5o)  de  continuité  sera,  en  1  et  a, 

(5:v  ^  +  rff  =  —  dx' 

»  27.  A  une  première  approximation,  la  dérivée  de  vs  en  x  est  ici  négli- 
geable, comme  se  trouvant  de  l'ordre  des  dérivées  secondes  de  U  et  c, 
tandis  que  v,  w  et,  par  suite,  les  parties  de  "X,  y.  à  évaluer  actuellement, 
sont  comparables  aux  dérivées  premières  des  mômes  quantités.  L'équa- 
tion (53)  reviendra  donc  à  poser,  en  introduisant  une  nouvelle  fonction 
auxiliaire  $  de  t,  x,  r,  et  'C, 

<s«  *  =  %•      v-  =  ~W 

»  Alors  la  condition  spéciale  (4<))>  ou  'es  deux  dérivées  partielles  de  6 
en  r,,  '(  sont  entre  elles  comme  les  deux  projections  respectives  d'Ç,  —'di\ 
d'un  élément  du  contour  dans  la  section  type  (a  ==  r,  h  =  i),  deviendra 
d$  =  o  (le  long  de  cet  élément)  et,  par  suite,  <I>=  une  constante  C  sur  tout 
le  contour  d'une  môme  section,  si  du  moins  l'on  admet,  comme  nous  le 
ferons,  que  ce  contour  soit  tout  d'une  pièce,  ou  que  la  section  a  n'ait  pas  à 
son  intérieur  de  lacune  non  occupée  par  le  fluide.  Or  l'on  peut  ajouter  —  C  à  <1» 
sans  modifier  les  dérivées  en  r„  C  de  cette  fonction,  ni,  par  suite,  X,  y.. 
Nous  pourrons  donc  astreindre  <I>,  dans  chaque  section  en  particulier,  à 
vérifier  la  relation 

(55)  <!>  —  o  (sur  le  contour  de  la  section  cr). 

))  D'ailleurs,  cette  petite  fonction  $  (de  l'ordre  de  grandeur  de  1,  a),  défi- 
nie dans  chaque  section  et  existant  par  suite  dans  toutes,  aura  sa  valeur 
variable  avec  x  et  /,  mais  assujettie  à  s'annuler,  comme  on  voit,  sur  toute 
la  surface  limite  du  cours  d'eau. 

»  28.  Avant  de  développer  la  condition  d'intégrabilité  (i),  qui  nous 
reste  comme  équation  indéfinie  pour  achever  de  déterminer  <]>,  observons 
que  cette  fonction  <î>  sera  identiquement  nulle  dans  les  deux  cas  simples  : 
i°  d'un  courant  à  section  rectangulaire  d'une  très  grande  largeur  con- 
stante i a  et  de  rayon  moyen  h,  cas  où  r0  est  une  constante  et  où  à,  F, 
dépendent  uniquement  de  '(;  2°  d'un  courant  rectiligne  à  section  circulaire 
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variable  (de  rayon  R),  cas  où  y0  =  z0  =  o,  a  =  h  =  ~R  el  où  r,,  'C  n'entrent 

dans  F,, -y  que  par  le  rappoi  t,  i  =  lyV-h  £*,  de  la  distance  r  à  l'axe  au 
ravon  R. 

)>  Alors,  en  effet,  le  mouvement  se  faisant  dans  les  plans  ou  normaux 
aux  y,  ou  menés  suivant  l'axe  des  x,  et  de  la  même  manière  dans  tous, 
v,  w  sont  les  produits  respectifs  d'une  même  fonction  soit  de  C  soit  det, 
par  zéro  et  i ,  ou  par  r,  et  "C;  et  ce  sont  par  conséquent,  dans  les  deux  cas, 
les  dérivées  en  r\  et  en  '(  d'une  même  fonction,  dépendant  soit  de  'C,  soit 
de  i.  En  outre,  les  équations  (02)  et  (7j8)  ne  font  aussi,  par  raison  de 
symétrie,  dépendre  y  que  de  'C  ou  de  x.  Dès  lors,  X,  ;a,  tirés  de  (4^>),  appa- 
raissent également  comme  étant  deux  telles  dérivées,  en  r\  et  '(,  d'une 
même  fonction  soit  de  'C,  soit  tlet.  Dans  le  cas  de  la  section  rectangulaire 
large,  a  étant  ainsi  nul,  la  première  équation  (54),  combinée  avec  la  con- 
dition (jd)  au  contour,  donne  (I>  =  o.  Dans  le  second  cas,  l'égalité  des 
deux  dérivées  respectives  de  a  en  "C  et  de  ;x  en  n  entraine,  vu  (54),  l'équa- 
tion A2$=o,  revenant  bien  à  <I>  =  o  par  suite  de  la  condition  au  con- 
tour (5  3). 

»  Dans  les  deux  cas  si  importants  dont  il  s'agit,  les  petites  composantes 
transversales  e,  w  de  la  vitesse  seront  donc  complètement  représentées, 
à  une  première  approximation,  par  les  formules  (46),  prises  avec  X,  pt,  nuls, 
et  où  /.,  y  auront  les  valeurs  définies  par  (48)  et  (02). 

.§  XIV.  —  Leur  emploi  dans  la  formation  de  l'équation  du  mouvement. 

»  29.  Même  dans  les  autres  cas,  les  parties  de  e,  w  dépendant  de  <I> 
ne  donneront,  dans  la  formule  de  l'accélération  longitudinale  u! ',  que  des 
termes  nuls  en  moyenne  sur  l'étendue  de  la  section  g,  et  v  restant  tels  après 
avoir  été  multipliés  par  une  puissance  quelconque  o'""1  de  la  fonction  <p: 
particularité  qui  les  élimine  des  équations  définitives  (17),  (18)  du  mou- 
vement et  réduit  bien  leur  importance. 

»  Formons,  en  effet,  l'expression  approchée  (41)  de  u  ,  en  substituant 
à  v,  w,  ).,  a  leurs  valeurs  respectives  (46).  (54)  et  aux  dérivées  de  o(rn  £) 
les  leurs,  résultant,  par  des  formules  comme  (44).  des  définitions  (42) 
de  7),  'Ç.  Il  viendra  immédiatement 

»   Or,  si  l'on  multiplie  le  dernier  terme,  dépendant  de   la  fonction  <P, 
B. 
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par  o'"-{  c/g,  et  qu'on  intègre  dans  toule  l'aire  c  ou  plutôt  dans  l'étendue  de 
la  seclion  type  ayant  r,,  *(  comme  coordonnées  y  — -  y0,  z  —  z0,  l'intégrale* 
obtenue  e^t(à  un  facteur  constant  près),  môme  pour  une  section  autre  que 
la  section  type, 

ff&r  %  -  *£S**-//K(t-$)  -  *(•-£)]** 

»   Transformons,   à  la  manière  ordinaire,  cette  intégrale  de  surface,  à 
parties  inlégrables  une  fois,  en  une  intégrale  prise  le  long  du  contour  -/'. 

I /élément  de  celle-ci  aura  évidemment  le  facteur  f-^rcosp  —  -r-sin  (3)  <Y/', 

identiqueà  la  différentielle  r/<I>  prise  le  long  du  contour  y'  et  nulle  en  vertu 
de  (55). 

»   Ainsi,  la  partie  de  w'  qui  dépend  de  la  fonction  <E>  a  bien  son  produit 
par  ç'""1  nul  en  moyenne. 

»   50.   Quant  à  la  partie  précédente,  où  figure  y,  la  valeur  moyenne  de 
son  produit  par  o'"-1    s'obtient  aisément,  en  observant  que  celte  valeur 

moyenne  est  celle  de — —  A2y.  En  effet,  le  produit  considéré,  ajouté  à 

l     KCS  '"  l     / 

— —  A2v»  donne,  au  facteur  près  — -, 


»  Or,  on  voit  que  cette  dernière  quantité,  multipliée  par  dridÇ  et  intégrée 
dans  toute  l'étendue  de  la  section  type,  se  convertit  en  une  intégrale  de 
contour,  à  élément  nul  en  vertu  de  (48).  Elle  a  donc  bien  zéro  pour  valeur 
moyenne;  et,  comme  il  résulte,  d'autre  part,  de  (5i)  et  (5a),  ou  de  l'an- 
nulation de  l'ensemble  des  deux  premiers  termes  de  (5o),  que 

la  valeur  moyenne  considérée  devient  celle  de 
(57)  _  JL*/  >*m         »\ 

»  Par  suite,  l'expression  (56)  de  u'  donne  immédiatement,  en  désignant 
encore  par  le  symbole  OR  la  valeur  moyenne,  sur  toute  l'étendue  c,  de  la 
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quantité  écrite  à  la  suite, 

(58)     ai(?-«iO  =  f^ll.f  +  U^a.f-  -  iL  ^(Oll.o—  -  3iL.?«). 

»  Dans  les  deux  cas  où  m  =i  et  où  t?2  =  2,  les  valeurs  moyennes  de  om 
et  o'"+t  sont  respectivement  celles  de  o,  9-  et  <ps,  c'est-à-dire,  pour  la  pre- 
mière, l'unité,  et,  pour  les  deux  dernières,  les  quantités  que  nous  avons 
appelées  i+r„  a,  mais  réduites  tout  de  suite  à  leurs  valeurs  de  régime 
uniforme.  Et  alors  les  équations  (17),  (18)  régissant  le  mouvement  gra- 
duellement varié,  dans  lesquelles  (ir )'  n'est  autre  que  11m  ou,  sauf 
erreur  négligeable,  i\Jou\  prennent  les  formes  (23),  (26)  simplifiées, 
auxquelles  nous  étions  parvenus  autrement. 


§  XV.  —  Mouvement  transversal  tournant,  dans  un  écoulement  permanent 

graduellement  varié. 

»  31.  L'importance  de  <I>  se  trouvant  ainsi  très  diminuée  par  le  double 
fait  que  cette  fonction  s'annule  dans  les  deux  formes  les  plus  intéressantes 
de  section  et  s'élimine  (à  une  première  approximation)  des  équations  défi- 
nitives du  mouvement  en  1,  7,  U,  u0,  nous  nous  contenterons  de  former 
l'équation  indéfinie  en  fI>  dans  l'hypothèse  d'un  régime  permanent,  avec 
axe  hydraulique  rectiligne,  c'est-à-dire  avec  y0  =  o  et  r()  =  o.  Alors  les 
accélérations  latérales  p',  \v'  se  réduisent  aux  produits  des  petites  dérivées 
de  c,  w  en  x  par  u  ou,  très  sensiblement,  par  Uo.  Or  les  expressions  (  46) 
de  v,  w  sont  elles-mêmes,  vu  (54)  et  la  première  (52), 

(39)         «■  =  i<  (asfi  +  «aj j.      <*  =  u  (35 ?c  -  *  7/- 

»    Sans  les  termes  en  <P,  le  rapport  —  égalerait  -, —  ou  T  -, ,  c  est-a-dire 

simplement  le  rapport  même  -  des  coordonnées  correspondantes  dans  le 

cas,  que  l'on  a  particulièrement  en  vue,  de  sections  toutes  semblables 
où  a  et  h  varient  avec  x  proportionnellement  à  leurs  valeurs.  La  vitesse 
transversale  résultant  fie  c,  kv  produirait  donc  un  mouvement  centrifuge 
ou  centripète  par  rapport  à  l'axe  hydraulique  (y  =  o,  z  =  o),  mais  nulle- 
ment rolaloire  autour  de  celui-ci.  Donc  la  fonction  $  exprime  ce  qu'on 
peut  appeler  le  mouvement  transversal  tournant  du  fluide.  Et,  en  effet,  les 
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parties  de  v,  w  qui  dépendent  de  <P,  elant  entre  elles  comme  —r->  —  -y-  (a 

la  traversée  de  sections  semblables),  représentent  des  vitesses  dirigées 
suivant  les  courbes  <ï>  =  const.,  qui  entourent  bien  l'axe  hydraulique, 
puisque  la  plus  extérieure  d'entre  elles,  celle  qui  a  l'équation  <P  =  o,  se 
confond  avec  le  contour  mouillé  de  la  section. 

»  La  différentiation  de  v,  w par  rapport  à  x  s'effectuera  en  ne  regardant 
comme  variable,  dans  chaque  terme  des  expressions  (5q),  que  son  petit 
facteur;  car  tout  autre  facteur  que  l'on  y  ferait  varier  introduirait  à  sa 
place  une  dérivée  très  petite,  dont  le  produit  par  le  facteur  déjà  petit  du 
terme  serait  de  Tordre  des  quantités  que  l'on  néglige.  Si,  pour  abréger, 
nous  désignons  au  moven  d'accents  les  dérivées  en  x  de  a,  h  et  <I>,  il  vien- 
dra, après  multiplication  par  L'o, 

/r     \  >  ri2/     „     2  (hv'\  >  TTo/i,/     OY  I       «?*' 

( bo)  v  =  \j-  (a  çr  r,  -h  a  o  -p-  I  '         w  =  U  -  (  h  <p-  Ç  —  ny  -j- 

»  L'équation  à  former,  exprimant  l'égalité  des  deux  dérivées  de  v'  e:i  r 
et  de  w'  en  y,  sera  après  suppression  du  facteur  U2, 

h    d   f    d<P'\         a    d  (    d<t>'\  /ik"  vdo         la"     d 


Y 


r 


(61)  a  dr,  [y  dr,  J    '     h  dX,  VY  dl  )  "  1  \    a    *  di\  h     "d 

Ce  sera  l'équation  cherchée  en  *!>',  si  Ton  v  substitue  à  <p  son  expression 
résultant  de(5ï).  Il  faudra  d'ailleurs  y  joindre  la  condition  définie,  dé- 
duite de  (55), 

(62)  $'=  o     (le  long  du  contour). 

En  effet,  <I>  s'annulant  sur  toute  la  surface  enveloppe  du  fluide,  l'on  a 
rtVp  =  o,  le  long  du  chemin  que  suit  toute  particule  de  la  couche  superfi- 
cielle, c'est-à-dire  quand,  à  partir  d'un  point  de  la  surface,  ce, y,  z  croissent 
deudt,  vdt,wdt.  Mais  les  produits  par  e,  w  des  dérivées  en  y,  z  de  la  petite 
quantité  <ï>  sont  de  l'ordre  des  termes  non  linéaires  que  Ton  supprime;  et 
l'équation  d<1>  =  o  rend  ainsi  négligeable  la  dérivée  (I>'  de  <I>  en  x  à  la  sur- 
face limite. 

»  52.  Les  équations  (Gi),  (62)  déterminent  complètement  <ï>'  quand  a, 
h  et,  par  suite,  le  second  membre  de  (61)  sont  censés  donnés  pour  toutes 
les  valeurs  de  x.  Car,  si  l'on  remplace  <p'  par  $'  -h  $', ,  il  vient,  pour  déter- 
miner <D'4,  la  condition  <!>',  =  o  au  contour,  avec  une  équation  indéfinie 
ayant  son  premier  membre  pareil  à  celui  de  (61),  mais  zéro  comme  second 
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membre.  Or  celle-ci,  multipliée  par  $',  dt\  dX,  puis  intégrée  dans  toute  la 
section  après  qu'on  a  remplacé 

.,   d  (    cW.\        . ,   d  (    d&.\  d  (,,       rf*'.\  (d^\\- 

donne  deux  termes,  intégrantes  une  fois,   que  réduit  à  zéro  l'annulation 
de  $'.  sur  le  contour,  avec  deux  autres  termes,  à  somme  dès  lors  nulle, 


<G3>  -// 


h  (d*\\*        a  /d4>\V 


o  dr,  dX,  —  o. 


Comme  on  a  évidemment  o>o  partout,  celle-ci,  (63),  oblige  d'annuler 
les  deux  dérivées  en  r,  et  'Ç  de  la  fonction  <!>,,  dès  lors  nulle  elle-même  à 
l'intérieur  non  moins  que  sur  le  contour. 

»  Quand  l'intégration  du  système  (61),  (62)  aura  fait  connaître  $',  il 
viendra  $  =  f$'  dx,  valeur  déterminée,  à  une  fonction  arbitraire  près 
de-/;,  (.  Or,  celle-ci  se  déterminera  elle-même  par  la  condition  que  <I>  s'an- 
nule aux  endroits  où  a,  A,  et,  par  suite,  g,  U  deviennent  constants  :  car  le 
régime  qu'on  étudie  est  supposé  devenir  uniforme  dès  que  a,  U  ne  varient 
plus;  et  v,  w,  <ï>  môme,  d'après  (39)  et  (55),  se  réduisent  alors  à  zéro. 

»  55.  Les  sections  où  $'  s'annulera  dans  l'état  permanent,  les  seules 
où,  par  suite,  $  puisse  s'annuler,  seront,  d'après  (61),  celles  où  l'on  aura 

rrr\  1      d'i  l      d'î 

(04)  ■ n 7 7TT7  ZJ?  —  O. 

x       J  aa  y,  dt\         lilrt  dC 

»  L'intégration,  immédiate,  de  cette  équation  aux  dérivées  partielles  en 
rt,  C»  montre  que  o  est  alors,  dans  chaque  section  en  particulier,  fonction 
de  r,,  '(  par  l'intermédiaire  de  la  variable  unique  aa'r,2  -+-  kh"'C,2.  Supposons 
d'abord  que  celte  fonction  ne  se  réduise  pas  à  une  constante,  c'est-à-dire, 
d'après  (5i),  que  le  coefficient  de  frottement  extérieur  B0  ne  soit  pas  infi- 
niment petit.  Alors  les  courbes  o  =  const.,  d'égale  vitesse  dans  le  régime 
uniforme,  devront  donc  être  des  coniques  semblables  et  semblablement 
placées,  à  centre  commun  situé  sur  l'axe  hydraulique  (7)  =  o,  £  =  o)  du 
courant,  ellipses  ou  hyperboles  suivant  que  les  deux  dérivées  secondes  de 
a  et  h  en  x  auront  ou  n'auront  pas  même  signe. 

»  On  ne  connaît,  pour  les  écoulements  tourbillonnants  étudiés  ici,  que 
les  deux  cas  de  sections  circulaires  (ou  demi-circulaires)  et  de  sections  rec- 
tangulaires larges,  où  les  courbes  d'égale  vitesse  dans  le  régime  uniforme 
soient  ainsi  des  coniques;  et  même,  dans  le  second  cas  où  o  ne  dépend 
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effectivement  que  de  £  ou  de  Ç2,  la  variable  de  <p  exigée  par  (64)»  qu'on 

aa" 
peut  supposer  être  'Ç2  -h  7777  v?  ne  se  réduit  à  Ç2  que  si  le  produit  aa"  égale 

une  fraction  négligeable  de  M",  c'est-à-dire  seulement  dans  l'hypothèse 
d'une  largeur  2«  constante  (en  écartant  celle  d'une  largeur  uniformément 
croissante  ou  décroissante,  qui  excluerait  partout  la  possibilité  d'un  régime 
uniforme).  Ainsi,  dans  un  canal  à  largeur  très  grande,  mais  variable,  les 
expressions  (5r>)  de  v,  w  ne  se  réduisent  pas,  en  général,  à  leurs  premiers 
termes,  simples,  comme  on  aurait  pu  l'espérer  ('). 

»  Mais  il  reste  le  cas  extrême  de  parois  assez  polies  pour  qu'on  puisse 
écrire  approximativement  o  =  i,  c'est-à-dire  supposer  tous  les  filets  fluides 
à  peu  près  également  rapides;  cas  où  l'équation  (61),  ayant  son  second 
membre  négligeable,  admet  pour  solution  <I>''  =  o  et,  par  suite,  $  =  o. 
Donc,  alors,  pour  toute  forme  de  section,  c,  w  dépendent  linéairement, 
d'après  (5g),  des  coordonnées  transversales  rt,  C.  C'est  ce  que  nous  verrons 
bientôt  plus  complètement  (2). 

»  oi.  Bornons-nous  maintenant  à  la  supposition  de  sections  5  toutes 
semblables,  pour  laquelle  seule,  sauf  le  cas  particulier  ci-dessus  de  rec- 
tangles larges,  a  été  établie  la  formule  (5i)  de  <p  —  1 .  Alors,  afin  d'em- 
brasser aussi  ce  cas  particulier  dans  lequel  a  représente  la  demi-largeur, 
évitons  de  poser  a  =  h,  et  admettons  seulement  que  a  soit  le  produit  du 
rayon  moyen  h  par  une  constante.  Le  second  membre  de  (61)  deviendra 
plus  symétrique  en  y  faisant,  comme  on  le  peut  évidemment, 

ia"  __  i_a_  a  _   (<£_  __  h"\  a  <xh"  __   fa"  __  //"\  h 

h  a     h  \  a  h  j  h  a  \  a  h   !  a 


(')  Toutefois,  si  les  courbures  des  filets  fluides,  ou  les  accélérations  latérales  v',  w' , 
devenaient  assez  petites  pour  que  leur  influence  s'abaissât  à  l'ordre  de  petitesse  des 
termes  non  linéaires  négligés  dans  notre  analyse,  tous  les  résultats  basés  sur  nos 
équations  (1)  et  surtout  (61)  seraient  évidemment  illusoires,  la  vraie  équation  en  * 
devenant  beaucoup  plus  compliquée. 

(2)  Dans  le  cas  opposé  de  parois  assez  rugueuses  (ou  d'une  valeur  de  B0  assez 
grande)  pour  annuler  u  et  tp  à  la  paroi  comme  dans  les  mouvements  bien  continus, 
léquation  (61),  développée  en  effectuant  les  diflérentiations  indiquées  à  son  premier 
membre,  ne  laisse  subsister  aux  parois  que  les  deux  termes  où  figurent  les  dérivées 
premières  de  ç  en  tq,  Ç,  multipliées  par  celles  de  <!>'.  L'annulation  d'une  de  ces  der- 
nières entraîne  donc  celle  de  l'autre  ;  et  la  condition  <!''—  o,  d&  —  o,  le  long  du  contour, 
y  donne  en  tous  sens  £/*!>'  =  o,  puis  <ï>  =  o,  au  sein  de  la  couche  mouillant  la  paroi. 
Autrement  dit,  et  vu  les  formules  (5g),  les  deux  conditions  distinctes  v  =  o,  w  =  o  se 
trouvent  d'elles-mêmes  vérifiées  dans  cette  couche,  comme  il  le  fallait  bien  physiquement. 
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»   Portons  en  même  temps  dans  (61)  l'expression  de  ç  fournie  par(5i) 
et  posons  enfin 

(65) 


$'  = 


,,..,    v  II         Cl 

(bO)    --r- 


i  +  A'v'Bo^F. 
»   L'équation  indéfinie  (61)  deviendra,  en  T, 

dv 


AV'B0 


H- F, 


rfr. 


a  cl 
h  dt 


vÂv/U 


I", 


rfS 


A  i/B, 


-h  F 


h    y    dF 


a      dFj 
h  r'  7/T 


relation  où  a,  h  n'entrent  que  par  leurs  rapports,  indépendants  de  x;  de 
sorte  que  a?  n'y  figure  pas  plus  que  dans  la  condition  spéciale  au  contour, 
devenue  r  =  o.  La  nouvelle  fonction  T  dépend  donc  uniquement  de  r„  £, 
du  paramètre  k\B0  et  de  la  forme  de  la  section.  C'est  ce  que  nous  spéci- 
fierons, en  écrivant  T  (y,,  '(,£v150)  au  lieu  de  T. 

»  Enfin  l'équation  (65),  multipliée  par  dx  et  intégrée,  donnera  pour  <I> 
la  même  formule  que  pour  $',  avec  simple  remplacement  de  a",  h"  par 
a  ,  h':  car,  dans  la  différentiation  de  celle-ci  par  rapport  à  x,  la  mise  en 
compte  de  la  variation  des  dénominateurs  a,  h  n'introduirait  que  des 
termes  non  linéaires  et  négligeables.  Il  n'y  aura  pas,  d'ailleurs,  à  ajouter 
une  fonction  arbitraire  der,,  'C,  puisque  $  doit  tendre  vers  zéro  aux  endroits 

où  y  tendraient  les  dérivées  a' ,  h'.  Et  si  nous  observons  que f-  -r-  =  -  -ï- > 

^       a        h         z  dx 

nous  aurons  la  valeur  définitive  cherchée 


(6?) 


«I»  = 


ky/ÏÏt 


1   di 


1  H-^s/BoDIcF,  "  dx 


rM,ksfïï0 


§  XVI.  —  Mouvement  transversal,  dans  l'écoulement  à  travers  des  sections 
ou  rectangulaires  d'une  grande  largeur  constante,  ou  circulaires. 

»  55.  Laissant  de  coté  la  section  rectangulaire  à  largeur  la  variable, 
calculons  e,  w,  dans  l'hypothèse  générale  d'un  mouvement  non  perma- 
nent, pour  les  trois  seuls  cas  où  nous  connaissions  o,  savoir,  ceux  de  la 
section  soit  rectangulaire  large(mais  avec  ia  constant),  soit  circulaire  ou 
demi-circulaire,  et  celui  d'une  section  de  forme  quelconque,  mais  alors 
avec  parois  assez  polies  pour  que  la  différence  z>  —  1  soit  de  l'ordre  des 
quantités  dont  nous  négligeons  les  carrés  et  produits. 

»    Dans  les  deux  premiers  cas,  (î>,  a,  y.  étant  nuls,  il  suffit  d'obtenir  y  par 


(%) 
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l'équation  indéfinie  (52),  A2y  =  F,  —  ;TtF,,  où  F,  dépend  seulement  soit 
de  '(,  soit  du  quotient  x  =  .7  \rl-  -t-  '(2  de  la  distance  r  à  l'axe  par  le  rayon  R, 
et  où,  par  suite,  y  ne  dépendra  aussi  que  de  (  ou  de  x.  On  aura  donc 


soit         ^=Ff-3lLF|f         soit         __^J=Fl-,mF,. 

»   Une  intégration  presque  immédiate  en  déduit,   vu  la  condition  (4$) 
annulant  la  dérivée  première  de  y  à  la  limite  '(  =  1  ou  t  =  1, 

(68)     g=f  (F,-D1LF,)^         ou         £  =  ;■    f  (F(-;)iiF,)^A; 

T  £ 

et  les  produits  de  cette  expression  soit  par  zéro  et  1 ,  soit  par  7A  et  -^,  seront 

les  dérivées  dey  en  r,  et  '(,  à  porter  dans  les  formules  (\6)  de  v,  <r,  ou  dans 
celle,  (56),  de  u\  où  9  est  donné  par  (5i).  Tl  vient,  par  exemple,  vu  la 
valeur  (01)  de  •/., 


TT  — 
dx 


"'^-dTV  +  Vdx-V2 


i  +  k\/B0DTl-Fj     ?  dt 


F',(  r  (F,  -OltF,)^,     ou     --  r(Fi—;-)KYi)tdx 


»  Il  ne  reste  plus  qu'à  y  effectuer  les  intégrations,  après  avoir  mis  pour 
F,,  DlLF,,  F'(  leurs  valeurs  effectives,  ^(1  — <?),  £,  — '(,  dans  le  cas  du 
rectangle  large,  et  |(i  —  ï3),  |,  —  ix'1  dans  le  cas  du  cercle  (avec  une 
approximation  suffisante). 


§  XVII.  —  Mouvement  transversal,  dans  l'écoulement  entre  parois  polies 
et  à  travers  des  sections  d'une  même  forme  quelconque. 

»  56.  Lorsque  o  —  1,  ou  le  coefficient  £\B0,  sont  de  l'ordre  des  quan- 
tités dont  on  néglige  les  carrés  et  produits,  le  rapport  des  vitesses  indivi- 
duelles u  à  leur  moyenne  U  reste  évidemment  voisin  de  1  dans  le  régime 
graduellement  varié,  même  aux  endroits  et  aux  instants  où  un  tel  régime 
se  détruit  ou  s'établit  rapidement,  c'est-à-dire  quand  Les  dérivées  succes- 
sives de  U  ou  de  n  en  x  et  /,  encore  petites,  ou  commençant  à  l'être,  n'ont 
pas  des  grandeurs  décroissantes  à  mesure  que  leur  ordre  s'élève.  En  effet, 
les  frottements,  même  alors,  sont  la  principale  cause  de  Yinëgalitê  de 
vitesse  entre  les  filets  fluides  et  n'agissent  pas  plus  (ou,  du   moins,   pas 
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incomparablement  plusj  pour  l'accroître,  quand  le  régime  varie  que  lors- 
qu'il est  uniforme. 

»  La  différentiation  complète  de  u,  v,  w  en  /,  pour  obtenir  u',  v' ,  w',  s'y 
effectuera  en  faisant  varier  /  de  cil  et  x  de  adt,  ou  même  de  U dt,  mais  en 
laissant  constantes  y,  z;  car  les  dérivées  en  v,  z  non  seulement  de  v,  w, 
mais  aussi  de  u,  v  seront  de  petits  facteurs,  comme  v,  w.  On  aura,  très 
sensiblement, 

/      \         /         du        r,du         dV         T,  r/U  /  ,       /N         dlv,  w)  d(v,w) 

et  la  condition  d'intégrabilité  (i)  deviendra 

,      x        /  <r/         TT  <-/  \  /rfi'        c?«A  dc  (  dv        dw  , 

(70     b>+u^;  b:-Â7  =°»      ou      HAdk-dï  =°' 


\  tft  «t/j? )  \dz        dy  J  dt  \ dz        dy 

équation  où  le  svmbole  dc  désigne  une  différentiation  complète  par  rapport 
au  temps,  effectuée  en  suivant  une  même  particule  fluide  dans  son  mou- 
vement moyen  local. 

»   Donc,  d'une  part,  l'accélération  longitudinale  u'  se  calcule  sans  avoir 

besoin  de  connaître  v,w.  D'autre  part,  la  différence  -r  —  -j-  reste  nulle  tou- 
jours, dans  chaque  particule  fluide,  comme  aux  moments  où  le  régime  est, 
autour  d'elle,  uniforme  ou  très  graduellement  varié;  et  v,  w  sont,  à  l'in- 
térieur de  chaque  section  normale  cr,  les  dérivées  respectives  en  y,  z  d'une 
même  fonction. 

»  D'ailleurs,  ct  étant,  comme  o  —  i,  constamment  voisin  de  zéro,  même 
quand  U  et  a  changent  entre  de  larges  limites,  le  second  membre  de  (53) 
continue  à  être  négligeable,  alors  que  v,  w  ne  le  sont  plus,  et  les  relations 
(54),  (55)  subsistent.  Enfin,  l'on  peut,  dans  les  expressions  générales  (46) 
de  v,  w,  où  tous  les  termes  avaient  déjà  un  petit  facteur,  réduire  9  à  l'unité 
et  supprimer  les  deux  termes  où  figure  /.,  devenu  un  second  pelit  facteur 
(de  Tordre  de  k\'b0)  d'après  (02). 

»  Alors  l'égalité  des  deux  dérivées  de  v  en  z  et  de  w  en  y  donne  immé- 
diatement, comme  équation  indéfinie  en  0,  que  complétera  la  condi- 
tion (55)  au  contour, 

h  d-<i>         a  d*<ï> 

a   dr*         h   dV- 

et  celles-ci,  exactement  pareilles  (vu  9  =  1)  aux  équations  (62)  et  (6r) 
B.  6 
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en  $',  déterminent,  comme  elles,  leur  solution,  qui  est  $  =  o.  Les  formules 
(46)  donnent  donc  simplement 

(72) 


I  dz0  dz0\  / dh        TTrf/( 

dt  dx  )         \dt  dx 


Y 


valeurs  linéaires  en  rt,  *(  ou  en  y,  z. 

»  Ces  expressions  (  72)  de  v,  w,  établies  sans  supposer  les  dérivées  de  U 
et  o  de  plus  en  plus  petites  à  mesure  que  leur  ordre  s'élève,  nous  permet- 
tront d'aborder  l'étude  sinon  des  écoulements  rapidement  variés,  du  moins 
de  ceux  qui  commencent  à  le  devenir  ou  qui  cessent  de  l'être. 


§  XVIII.  —  Distribution  des  vitesses  à  travers  des  sections  semblables, 
dans  les  régimes  graduellement  variés. 

»  37.  Connaissant  par  la  formule  (56)  l'accélération  longitudinale  u 
aux  divers  points  ( y,  z),  ou  mieux  (r„  Ç),  d'une  section  ?,  il  devient  pos- 
sible d'intégrer  le  système  (io)  d'équations  déterminant  la  fonction  F2, 
dans  l'expression  générale  (8)  du  mode  de  distribution  des  vitesses.  Nous 
savons  que  cette  fonction  F2  s'annule  avec  le  second  membre  u'  —  D)1iï  de 
la  première  (10).  Celle-ci  étant  d'ailleurs  linéaire,  il  est  clair  que  F2  se 
composera  d'autant  de  termes  qu'en  comprend  u' ,  et  respectivement  pro- 
portionnels aux  facteurs  indépendants  de  tj,  '(  dans  les  termes  de  u! .  Vu  les 
équations  (5i),  (52),  (48),  (67),  et  enfin  (66)  combinée  avec  la  condi- 
tion T  =  o  (sur  le  contour),  qui  régissent  9,  x.,  y,  <ï>,  nous  aurons  donc, 
dans  toutes  les  sections  semblables,  y  compris  même  les  sections  rectangu- 
laires d'une  très  grande  largeur  constante,  un  résultat  de  la  forme 


F2  = 

(?3) 


AyB0 


'     dU   g(n,K,k\/B0)    ,    dV 


r*.(l.C)] 


k  v/B„  WL  F,  |_  "  dx  l+  k  \lB 0  DK  F ,        dl 

;JB0DTlFJ     l*  dt     -^  dx 


i-hArv/B031LF, 


»  Les  quatre  fonctions  der,,  'C  appelées  §,  ,7,,  #2>  $»*  dont  la  première  et  la 
dernière  seules  dépendent  en  outre  dek\fB~ot  sont  déterminées  par  les 


(76) 
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quatre  équations  indéfinies  respectives 


(74) 


4(F|)   +  |(Fg  =  2(F,-3r,F,K*vB„(F;-or,F;), 

«fy  V     ~dn)  ^  d\ \    "rfç/  —  ~^T  <^î  +  ~^T ^  ~~  J1L  \"rfn" ~cfrt  +  ~dTdl, 
drt\      drj    ;    cK\      dX>)~    dK  lh\         drt    dl  \dK  dn         dr,   dX,)' 


et,  en  outre,  par  les  conditions  définies  communes 

/    ^N  )  r  7 =  o  (au  contour), 

il3)  j  dv 

(  (^,«^1,^0» #3)  =  o(au  milieu  du  fond). 

»  Toutefois,  la  partie  en  §z  de  la  formule (73),  celle  qu'introduit  la  fonc- 
tion $  (ou  la  fonction  T)  et  qui  serait  nulle  dans  les  cas  des  deux  sections 
circulaire  (ou  demi-circulaire)  et  rectangulaire  d'une  grande  largeur  con- 
stante, n'est  donnée  que  pour  le  mouvement  permanent,  auquel  nous  nous 
sommes  bornés  dans  le  calcul  de<ï>;  et  c'est  pourquoi  nous  avons  pu,  le 

débit  Ut  étant  alors  constant,  v  remplacer  le  facteur 7^  par  U-v-jde 

1  a   dx  r  dx 

manière  à  faire  rentrer  ce  terme  dans  le  type  de  celui  d'entre  les  termes 
précédents  que  n'annule  pas  l'hypothèse  de  la  permanence,  et  qui  est  le 
terme  en  §. 

»  58.  La  formule  (11),  caractéristique  du  mode  de  distribution  des  vi- 
tesses, deviendra  donc,  si  l'on  y  substitue  à  u0,  dans  les  petits  termes,  le  quo- 
tient de  U  par  1  4-  i£y  Bo£)ïlF,,  et  sauf  toujours  la  même  restriction  quant  au 
terme  en  $3, 

5  =  1  +  *  fà  01LFf  -h  kj°-  £  [U  g  *(,,  U  y^)  "H  (1  +  h  v/B^OlLFj  §  l  (,,  O] 


dV 


j 


,(d,  0-+-U-t-#3(yj,  Lksf^9/ 


dt  v^'"  *'    '    ~  dx 


»  Dans  les  cas  particuliers  de  sections  rectangulaires  d'une  grande  lar- 
geur constante  et  de  sections  circulaires  ou  demi-circulaires,  l'on  a  §3  =  o, 
les  valeurs  de  J,  ,7,,  32  se  calculent  aisément  par  les  équations  (74)  et 
(73),  où  F,  et  y  ont  les  valeurs  indiquées  plus  haut  (p.  4o),  avec  F  égal 
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soit  à  i,  soit  à  l'inverse  de  ï  ;  et  l'on  obtient  ainsi,  pour  le  rapport  de  u 
à  u0,  les  expressions  données  aux  §§  IX,  X,  XXVI  et  XL  de  mon  Essai  sur 
la  théorie  des  eaux  courantes  (p.  90,  94,  246,  266,  et  5 16  à  $20)  ('). 

§  XIX.  —  Équation  du  mouvement  graduellement  varié, 
aux  degrés  d'approximation  supérieurs. 

»  59.  Le  second  membre  de  (76),  divisé  par  sa  valeur  moyenne  aux 
divers  points  d'une  section  5,  donnera  le  rapport  cp  -+-  n>  de  //  à  U.  On  for- 
mera son  excédent  sur  l'expression  de  o  résultant  de  (5i)  ;  et  cet  excédent, 
réduit  à  sa  partie  linéaire  par  rapport  aux  trois  petites  dérivées  premières 
de  L  en  x  et  de  U  et  g  en  /,  sera  la  fonction  rr  dans  sa  partie  de  première 
approximation,  ou  abstraction  faite  d'écarts  comparables  aux  dérivées 
d'ordre  supérieur  de  U  et  a.  On  pourra  donc  évaluer  les  petits  excès  res- 
pectifs 2<pcr,  3cp2c7  du  carré  et  du  cube  de  cp  -h  cr  sur  ceux  de  cp;  et  leurs  va- 
leurs moyennes  dans  toute  l'étendue  a,  savoir  Dli (2977),  OU  (3 o2 ci),  seront 
les  petites  parties  variables  des  coefficients  1  -+-  r,  et  se,  réduites  à  leurs  termes 
principaux  ou  affectés  des  dérivées  premières  de  U  et  a.  Ce  qui  s'y  trouve 
ainsi  négligé,  étant  d'un  ordre  supérieur  au  premier,  aura  ses  dérivées  en  x 
ou  en  /  d'un  ordre  supérieur  au  second  et,  par  suite,  négligeable,  môme 
à  une  deuxième  approximation  des  lois  du  mouvement  graduellement  varié. 

»  Il  suit  delà  qu'il  suffira,  à  une  deuxième  approximation,  de  substituer 
dans  (23)  et  (26),  aux  dérivées  de  r,  et  x  qui  y  figurent,  les  dérivées  ana- 
logues des  expressions  trouvées  pour  OlL(acpcT)  et  pour  0Tt(3cp2nj).  C'est 
ainsi  que,  dans  le  cas  d'un  cours  d'eau  à  section  rectangulaire  d'une  grande 

(')  On  trouve,  dans  le  cas  de  la  section  rectangulaire  large, 

*-=-  —  u-vy-'^d-Lv^  5^_lr«V        $  =__!(,_  ?»)*=_  I  F2 

12  40    \         3  '        3  "        3  '  /  '  24  6     '  ' 

36o 

et,  dans  le  cas  de  la  section  circulaire,  où  §,  §u  J2  ne  sont  évidemment  fonction 
de  7),  £  que  par  l'intermédiaire  de  v, 

^  =  -£(i-,3)2-^^(i4-33t3+24*G-5^),       £,=-£(i-t»)«  =  -Jfj, 
4o  9.  .\o  (\o  0 


,!• 


2 


"'2  —  ~ 

9'7D 


—  (2  —  9' 
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largeur  constante,  l'on  arrivera  le  plus  simplement  possible  à  l'équation  de 
deuxième  approximation  du  mouvement,  qui  porte  le  n°  367  dans  mon  Essai 
sur  la  théorie  des  eaux  courantes  (§  XXXVI,  p.  437)  et  dont  dépend  la  dé- 
formation plus  ou  moins  rapide  des  ondes  descendantes  ou  ascendantes  (de 
courbure  sensible)  le  long  d'un  tel  courant  ('). 

»  40.  La  différentiation  en  x  de  la  valeur  trouvée  de  nj  fera  de  même 
connaître,  avec  erreur  comparable  aux  dérivées  troisièmes  seulement  de 
U  et  a,  le  second  membre  de  l'équation  (53),  dans  sa  partie  principale,  ou 
du  deuxième  ordre  de  petitesse  ;  et  l'on  pourra,  dès  lors,  aborder  le  calcul 
de}.,  m.,  c,  wt  u',  F2,  gt,  1  H- y),  a  pour  les  termes  de  cet  ordre.  Dans  les 
deux  cas  de  la  section  rectangulaire  d'une  grande  largeur  constante  et  cir- 
culaire ou  demi-circulaire,  e,  w  continuent  évidemment  à  être  les  dérivées 
en  y,  z  d'une  même  fonction  ou  de  "C  ou  de  x  ;  et  les  calculs  n'offrent  guère 
d'autre  difficulté  que  leur  excessive  longueur,  comme  on  peut  en  juger  par 
la  partie  citée  ci-dessus  de  mon  Essai  sur  la  théorie  des  eaux  courantes 
(§  XXXVI).  Je  les  y  ai  effectués,  dans  le  premier  de  ces  cas,  pour  arriver 
à  l'équation  (n°  367)  citée  ci-dessus,  avant  d'avoir  découvert  le  procédé 
qui  la  déduit  de  (23). 

»  Quand  on  emploie  ces  calculs  pour  former  1  -+-•/),  x  et  l'équation  du 
mouvement  à  une  troisième  approximation,  là  où  ils  sont  généralement 
indispensables,  il  ne  faut  pas  oublier  que  les  équations  (25)  et  (26) 
doivent  être  complétées,  comme  on  l'a  vu  au  début  de  cette  Etude,  par 
un  terme  provenant  de  ce  que  la  pression  moyenne  p  ne  varie  plus  alors 
hvdrostatiquement  dans  l'étendue  d'une  même  section  t.  A  partir  de  l'axe 
hvdraulique  où  l'on  se  donne  p  =  p0,  p  s'accroît,  en  effet,  du  terme 

(77)  —  p  /      iy'dy  -t-  w  dz)  =  —  p  I       (av'dr,  -h  hw'dÇ). 

J  i-    -  *Jù  a 


»  Au    terme  —  de  la  première  équation   indéfinie  du  mouvement  (  où 

S  V 

u'  et  p  entrent  par  l'expression 1 -j-  h  il  vient  donc  s'adjoindre, 

g         pg  ax  ) 

quand  on   élimine  de  cette  équation  la  dérivée  de  p  en  x  divisée  par  pg, 


(')  Voir  aussi  les  Additions  à  l'Essai  sur  la  théorie  des  eaux  courantes,  p.  55,  au 
Tome  suivant  XXIV  du  Recueil  des  Savants  étrangers. 
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l'expression 

»  On  a  pu  la  différenlier  sous  le  signe  /  et  même  n'y  différentier  que  le 
facteur  déjà  petit  v'  ou  w\  car  les  dérivées  en  x  des  limites  soit  inférieures 
y0,  z0,  soit  supérieures  r,,'Ç,  ou  celles  dea,7i,  ne  donneraient,  multipliées 
par  les  petites  fonctions  v'  ou  w',  que  des  produits  négligeables. 

»  Par  suite,  dans  nos  équations  (3),  (5),  (7),  (10),  (i3),  (i4).  (ID)> 
(17),  (18),  u',  partout  où  il  figure,  se  trouve  accru  du  produit  de  (78) 
par  g,  et  («a)',  c'est-à-dire  2 U<p u',  se  trouve  lui-même  accru  du  produit 
de  2L9  par  (78)  et  par  g.  Donc,  le  second  membre  de  l'équation  défini- 
tive (18)  ou  (^5)  du  mouvement  devra  être  complété,  en  y  ajoutant  le 
terme 

(79)  -  £»  [(2?  -  i)jf*  («a£*  +  h%^)\ 

où  il  suffira  d'évaluer  à  une  première  approximation,  par  les  formules  (4^) 
de  v,  w,  les  très  petites  accélérations  transversales  v' ',  w  . 


§  XX.  —  Passage  d'un  régime  graduellement  varié  à  un  régime  rapidement  varié, 

ou  vice  versa. 

»  41.  La  longueur  des  calculs  et  surtout  la  complication  des  résultats 
seraient  des  plus  rebutantes,  si  l'on  ne  se  bornait  au  cas  de  parois  assez 
polies,  ou  d'un  coefficient  B0  de  frottement  extérieur  assez  faible,  pour 
réduire  au  premier  ordre  de  petitesse  la  différence  9 —  1  et  aussi,  d'après 
(76),  les  inégalités  relatives  de  vitesse  des  filets  fluides  à  l'état  de  régime 
varié,  où  elles  sont  exprimées  par  <p  -+-  xs —  1.  Alors,  d'une  part,  (9  -h  cr)2 
et  (9  -+-  u>y  sont  très  sensiblement  1  -h  2(9  4-  xs  —  1)  et  1  -h  3(<p  -+-  ct  —  1). 
Leurs  valeurs  moyennes  ne  se  distinguent  donc  plus  de  1  ;  ce  qui  réduit  le 
second  membre  de  (20)  à  ses  trois  premiers  termes,  où,  même,  les  coeffi- 
cients 2<x  — -  1  —  y)  et  1  -1-  2?)  deviennent  l'unité.  D'autre  part,  les  compo- 
santes transversales  v,  w  de  la  vitesse  admettent  les  expressions  simples 
(72),  qui  donnent,  sous  une  forme  symbolique  et  abrégée,  mais  évidente, 
dans  laquelle  U  se  comporte  comme  un  facteur  constant  et  n,  Z,  comme  des 


(47  ) 


facteurs  indépendants  de  x  et  t, 


(80)    d{v'r']  =  a  - 


d  \2  [fdy0  da\        /dz,-,        y  dh 

dx  \dt         "  dx  1       \dx  '  dxj        \dx         **  dx 

»   Le  terme  complémentaire  (79),  à  joindre  au  second  membre  de  (20), 
devient  donc 


(8.) 


d 
dt 


"s)'* 


fdy0  da  rt-\         ,  /  dz0  y        dh  Ç2 

\  dx    '        dx  2  /  \  dx  ^    '     dx  2 


»  Par  exemple,  dans  les  deux  cas  :  i°  d'un  canal  rectangulaire  de  lar- 
geur ia  où  r,  varie  de  —  1  à  1  et,  £,  de  zéro  à  1  ;  i°  d'un  tuyau  circulaire 
où  y0  =  z0  =  o,  a  =  h  =  i  R,  et  où  le  rapport  t  =  i  \/r,2  -+-  C2  de  la  distance  r 
à  l'axe  au  rayon  R  varie  de  zéro  à  1,  ce  terme  (81)  devient  respectivement 


(82) 


dt 


v 


d 
dx 


a  dvn         h  dz0         a  da         h  dh  \ 

R  dR 

i  _) :  _) 1 .  ), 

2   dx          2  <^j?         6  <:/.r         6  rf^c/ 

4    «f  J7 

»  S'il  s'agit,  en  particulier,  d'un  canal  rectangulaire  de  largeur  constante, 
y0  et  a  sont  constants,  la  dérivée  de  z0  en  x  est  l'excédent  de  la  petite 
pente  actuelle  I  de  surface  sur  la  pente  constante  de  l'axe  des  x;  et  le 
terme  (82)  devient  aisément 


(83) 


d         TT    d  V  /r  '  dh 

dt  dx  J    \  3  dx 


»  42.  Par  l'adjonction  du  terme  (81),  (82)  ou  (83)  au  second  membre 
de  l'équation  (25)  et  la  réduction  de  1  -+-■/),  a  à  l'unité  dans  ce  second 
membre,  on  rendra  l'équation  (25)  applicable  à  nn  régime  qui  devient  ou 
qui  cesse  d'être  rapidement  varié,  ou,  encore,  qui  se  maintient,  sur  des  lon- 
gueurs notables,  voisin  d'un  régime  graduellement  varié;  c'est-à-dire,  en 
un  mot,  à  tout  régime  où  les  dérivées  de  g  et  U  en  x  ou  t  sont  petites, 
mais  sans  décroître  de  plus  en  plus  à  mesure  que  leur  ordre  s'élève.  La 
première  formule  (70)  et  les  formules  (72)  employées  dans  la  démonstra- 
tion ne  sont  basées,  en  effet,  de  même  que  les  transformations  opérées 
ci-dessus,  que  sur  l'hypothèse  de  la  quasi-égalité  relative  de  vitesse  des 
filets  fluides  et  sur  la  petitesse  commune  des  dérivées  de  U  et  a. 

»  Dans  le  cas  particulier  d'un  canal  rectangulaire  de  largeur  constante 
pour  lequel  a  été  obtenue  l'expression  (83),  il  suffit  d'observer  que  la 
pente  I  de  surface  égale  la  pente  (constante  ou  variable)  i  de  fond,  moins 
la  dérivée  de  h  en  x,  pour  rendre  cette  expression  (83)  identique  au  der- 


(  48  ) 
nier  terme  d'une  équation  (482)  donnée  pour  le  môme  cas  dans  mon  Essai 
sur  la  théorie  des  eaux  courantes  (p.  524).  On  peut  voir  aux  §§  XX  à  XXV 
de  cet  Essai  comment  l'adjonction  du  terme  dont  il  s'agit  à  l'équation  du 
mouvement  permet  d'étudier  l'état  permanent  d'un  cours  d'eau,  soit  aux 
points  où  un  régime  graduellement  varié  se  détruit  ou  s'établit,  comme, 
par  exemple,  au  pied  et  au  sommet  des  ressauts  brusques  ou  ondulés,  soit 
aux  endroits  où  le  fond  présente  des  ondulations  longitudinales  régnant 
sur  toute  la  largeur,  qui  se  répercutent  plus  ou  moins  à  la  superficie,  etc. 
»   Dans  tous  ces  cas,  l'équation  du  mouvement  permanent  est 

h*  u2  <jn      f  u2  _    d_  /  u»\   _  n-  <£_  /u» 

2    srh  dx-  h  dx  \  2  s  /  6   dx3  \  2  a 


On  le  voit  en  prenant  le  terme  complémentaire  sous  sa  forme  (83),  mais 

dh 

d.v 


en  v  substituant  à  la  petite  dérivée  -r-  (vu  la  constance  du  débit  UA  de 


l'unité  de  largeur  du  lit)  l'expression  —  —  -—  =  —     ,2      -(  —  j.  Et  elle 


devient 


•  ,  {': 1 2  d'[  _  / l  "  _  dh     si  (^1  \  __  —  —  (— 

2    trh   dx-  h  dx         f/.r  \  2  ;r  /  3    dxli  \  2  is 


quand  on  en  élimine,  comme  il  vient  d'être  dit,  la  pente  I  de  surface.  La 
première  forme  convient  surtout  pour  les  courants  sans  courbure  sensible 
de  superficie,  où  l'on  peut  supposer  nulle  la  dérivée  seconde  de  I  en  x  : 
c'est  le  cas  habituel  des  grands  cours  d'eau.  La  seconde  forme  donne 
l'équation  différentielle  en  x  des  variations  de  la  profondeur  h,  dès  qu'on 
met  pour  U  le  quotient,  par  h,  du  débit  q  de  l'unité  de  largeur. 

»  Ces  formules  s'appliquent  môme  sans  que  le  fond  ait  besoin  d'être 
très  poli  ;  car  on  n'a  négligé  dans  leur  établissement  que  les  carrés  et  pro- 
duits des  inégalités  de  vitesse  ou  d'autres  petits  facteurs  ('). 


(')  Équation  plus  approchée  d'un  régime  pseudo-uniforme.  —  La  possibilité 
d'employer  ainsi  l'équation  de  mouvement  obtenue,  même  sans  avoir  à  supposer  extrê- 
mement petit  le  coefficient  B  de  frottement  extérieur,  se  démontre  simplement  dans 
un  cas  particulier  digne  de  remarque.  C'est  celui  du  régime  que  j'ai  appelé  pseudo- 
uniforme, parce  que  la  vitesse  moyenne  U  et  la  profondeur  h  y  sont  constantes,  les 
filets  fluides,  quoique  courbes,  s'y  trouvant  parallèles  entre  eux  à  la  traversé*'  de 
toutes  les  sections.  Il  se  produit  quand,  sur  une  grande  longueur,  le  fond,  coupé  de 
sillons  transversaux  régnant  d'un  bord  à  l'autre,  affecte  une  forme  sinusoïdale,  d'une 
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§  XXI.  —  Du  régime  permanent  graduellement  varié  qui  se  produit  à  l'entrée 
ou  plutôt  dans  la  première  partie  amont  des  tuyaux. 

»  45.  Considérons  enfin  le  régime  permanent  varié,  très  spécial,  qui  a 
été  l'occasion  de  l'étude  actuelle,  savoir,  celui  qui  se  produit  dans  la  partie 
amont  d'un  long  tuyau  rectiligne,  après  l'épanouissement  des  filets  fluides 
consécutif  à  la  rapide  contraction  de  l'entrée,  et  qui  sert  de  transition  au 
régime  uniforme  existant  ensuite  sur  toute  la  longueur.  Dans  cette  ques- 

U 

longueur  complète  d'ondulation  égale  à   2-/1— —=.  (Essai  sur  la  théorie  des  eaux 

\J2gh 
courantes,  p.  233). 

Dans  ce  régime  pseudo-uniforme,  par  le  fait  même  qu'il  est  supposé  voisin  d'un 
régime  graduellement  varié,  la  distribution  des  vitesses  à  travers  les  sections  ne  diffère 
pas  beaucoup  de  celle  qu'exprime  la  formule  (76),  d'ailleurs  réduite  à  celle  d'un 

,  .        d(\],  1)  T  u  .  ... 

régime  uniforme  par  les  hypothèses  —jj- -  =0.   Le  rapport  yy  y  est  donc  sensible- 

ment  cp,  et  1  -+-tj,  a,  ne  différant  nulle  part  beaucoup  des  valeurs  moyennes  constantes 
de  o2,  cp3,  ont  leurs  dérivées  négligeables  dans  l'équation  (25),  dont  le  second  membre 

est  ainsi  réduit  à  6U2->  c'est-à-dire  à  6— •  Il  ne  reste  donc  qu'à  évaluer  son  terme 

a  h 

complémentaire  (79). 

A  cet  effet,  rappelons  que  les  filets  fluides  sont  supposés,  dans  tout  plan  vertical 

w 
longitudinal,  parallèles,  de  la   surface  au  fond  :  le  rapport—  reçoit  donc  partout  sa 

valeur  relative  au  fond,  savoir  i  —  c,  où  i,  c  sont  les  deux  angles  que  le  profil  longi- 
tudinal du  lit  et  l'axe  des  x  positifs  font  respectivement  avec  le  plan  de  l'horizon,  au- 
dessous  de  celui-ci.  On  a  ainsi 

.  TT  dw  di 

t  dw  ..       . 

et,  par  suite,  vu  la  suppression  permise  du  terme  w-r-  non  linéaire  en  w, 


dw  _        d£    2  dw^  __        d2i 

dx  dx  '   '  dx  dx2  ' 


Le  terme  complémentaire  (79)  devient  donc  (v'  étant  nul) 


AU2  dU 
S     dx 


?/'<*? ->(/*'*)*• 


L'hypothèse  simplificatrice  de  l'égalité  de  vitesse  de  tous  les  filets  fluides,  faite  dans 
B.  7 


(  5o) 

tion,  le  changement  des  vitesses  u  avec  l'abscisse  x  des  sections  <j  n'est 
plus  amené  par  des  variations  de  a  ou  de  la  vitesse  moyenne  U,  puisque  c 
et  U  sont  constants;  et  tous  les  termes  qui  prédominaient  jusqu'ici  dans 
nos  équations,  parmi  ceux  qu'introduit  la  non- uniformité  du  régime, 
s'effacent,  pour  laisser  le  premier  rôle  à  d'autres  beaucoup  plus  complexes. 
Nos  démonstrations  et  nos  formules  subsistent,  il  est  vrai,  sans  modifica- 
tion, jusqu'à  (^5)  et  (26)  inclusivement;  mais  les  seconds  membres  de 
celles-ci  perdent  leurs  termes  affectés  des  dérivées  de  U  et  de  g,  c'est-à-dire 
justement  ceux  que  la  simple  connaissance  des  lois  du  régime  uniforme 
permettait  d'évaluer;  et,  par  exemple,   l'équation  (25),  formule  générale 


le  calcul  du  terme  complémentaire  (83),  a  donc  eu  simplement  pour  effet  de  donner 
à  ce  terme  le  coefficient  numérique  l,  au  lieu  du  coefficient  plus  compliqué 


/("-"((V*)*. 


qu'il  aurait  eu  sans  cela  et  qui  se  réduit  bien  à  4  quand  on  prend  <p  ■=  1 . 
Or  on  a,  d'après  (36), 

?  =  H--V^(l-Ç«),  /\2^=r+iV^(ç_r3)4_^(5r_IO,3+9r:i)) 

2  JQ  6  Ibo 

et,  par  suite,  tous  calculs  faits, 

j(<.»-0(j[T*J«=3(.H-^--sdT> 

.Mais  on  voit,  par  les  formules  (36)  et  (3~),  que 
ce  qui,  en  doublant  d'ailleurs,  pour  plus  de  simplicité,  le  résultat  obtenu,  donne  enfin 
Le  coefficient  i  attribué  au  terme  complémentaire  devrait  donc  être  accru  de  la 

o . 

fraction  -7^(1  —  ^y5in)de  sa  valeur.  Or,  cette   fraction   se   trouve   comprise    entre 

4 

|t(  et  f{*j;  car,  vu  la  valeur  de  \Jb  donnée  par  la  première  formule  (3y)  (p.  33)  de 
mon  Etude  de  l'année  dernière  sur  le  régime  uniforme,  on  aura  toujours  k\/b<i  3  et, 
par  suite,  y5ï)  <  1,  1  —  -j^y/Sr,  >  -fa.  La  fraction  dont  il  s'agit,  toujours  positive,  sera 
donc  moindre  que  f  t\  en  valeur  relative,  c'est-à-dire  seulement  de  l'ordre  de  0,01  et 
insignifiante,  sauf  peut-être  dans  quelques  cours  d'eau  à  lit  extrêmement  rugueux. 


(  5.  ) 
du  mouvement,  devient  (vu  d'ailleurs  la  permanence  admise) 

(84)  i  =  ^u2^  +  ^-^(a-I-ï') 

\         ■  /  rr  tr 


g  dx 

m  44.  Toutefois,  celle-ci  n'est  pas  entièrement  suffisante  aux  faibles 
distances  de  l'entrée  du  tuyau  ;  car,  à  partir  de  la  formule  (i5),  nous  avons 
supposé  petites  non  seulement,  comme  la  définition  même  de  la  graduelle 
variation  du  régime  nous  avait  déjà  autorisé  à  le  faire,  l'accélération  longi- 
tudinale u' ,  les  vitesses  transversales  ç,  w  et  surtout  les  accélérations  cor- 
respondantes /,  w',  mais  aussi  la  différence  entre  le  mode  effectif  de  distri- 
bution des  vitesses  que  définit,  par  exemple,  le  rapport  o  -h  cj  de  u  à  U, 
dans  la  section  considérée,  et  le  mode  de  distribution  propre  au  régime 
uniforme,  exprimé  de  même  par  <p.  Or  la  petitesse  de  u'  entraîne  bien, 
d'après  le  système  (10),  celle  de  F2,  mais  non,  dans  la  formule  (1 1),  celle 
du  terme  en  1%,  à  cause  du  petit  dénominateur  y/B0  et  du  numérateur  assez 
grand  k  figurant  dans  le  coefficient  de  ce  terme.  Il  est  donc  possible,  en 
général,  que,  malgré  la  graduelle  variation  de  l'écoulement,  le  rapport  de 
//  à  uQ  donné  par  (r  1)  et,  par  suite,  celui,  <p  -h  u,  de  u  à  U,  s'écartent  très 
notablement  de  ce  qu'ils  sont  dans  le  régime  uniforme,  savoir,  de  <p  pour  le 

dernier  rapport  et  de  —  pour  le  premier,  <p0  désignant  la  valeur  de  o  au  mi- 

lieu  du  fond,  là  où  u  =  u0.  Ainsi,  ne  regardons  plus  comme  petite  la  fonc- 
tion cj,  ni,  par  suite,  sa  valeur  ct0  au  milieu  du  fond,  valeur  qui  est 


(vu  l'égalité  \^ou\0XLf=  bJJ2  dans  le  régime  uniforme);  et  nous  aurons, 
en  multipliant  (85)  par  Uy  \50D\l/, 


(8G)  sjB0ulDTtf  -  sfbU2  =  v'B0orc/Ucj0. 

»  La  différence  des  deux  premiers  termes  de  (16)  sera  donc  compa- 
rable à  chacun  d'eux,  et  l'on  ne  pourra  plus,  dans  (i5),  négliger  devant 
f  unité  le  carré  du  second  terme  entre  crochets;  mais  on  pourra  substituer 
partout  à  F,,  dans  ce  terme,  d'après  la  signification  même,  que  définit  (8), 
de  F,,  le  rapport  de  cp  —  o0  à  o0k  y!30 ,  ou  mieux  de  \  DTLf(y  —  <p0)  à  k  \fb . 

»   La  formule  (i5),  multipliée  par  u0  y b  ,  sera  identiquement 

/o    \  Il tto  rn — ' 7  1  <*  0Tl(ou' — u') 


(  52) 
»   Celle-ci,  comparée  à  (86),  montre  que  le  dernier  terme  de  (87)  égale, 
au  signe  prés,  s/B0DTLfUv50. 

»  Alors,  la  relation  (87),  élevée  au  carré,  donnera,  pour  tenir  lieu  de  (16), 

(88)  B.«:ai/=  6U»  -h  2  '-  M(*u'-»')  _  (bodil/)IP< 

»  Au  second  membre,  9  peut  d'ailleurs,  identiquement,  être  remplacé 
par  le  rapport  de  u  —  Uw  à  U  et  non  plus  simplement,  comme  on  avait 
fait  dans  (16),  par  celui  de  u  à  U;  ce  qui  ajoute  à  3îl(2Mm'),  qui  était  la 
valeur  approchée  de  2U3ÎL(<pw'),  la  correction  2U3ïl(—  au').  A  la  fin  de 
la  formule  (17),  il  faudra  donc  ajouter  l'expression 


2   a 


(89)  -  -  011  (-  uiï)  -  (B0DÏL/)U2< 

o    /. 

»   Par  suite,  cette  expression  se  retrouvera  en  plus,  divisée  par  le  rayon 
moyen,  dans  les  formules  (18),  (13)  de  la  pente  motrice  I. 


;  XXII.  —  Hauteur  motrice  qu'y  dépense  l'établissement  du  régime  uniforme. 

»   45.  L'équation  (84)  du  mouvement,  ainsi  complétée,  devient 

—  1 

dx 


(90)  1  =  bv  »  +  !£  d{*  .'    ^  + 1 ai(-  *«*)  -  (n,,*/)  $  uv;. 


»  La  hauteur  motrice  totale  /  I  dx  dépensée  entre  deux  sections,  abais- 
sement, entre  elles,  tant  de  l'axe  hydraulique  que  de  la  pression  sur  cet 
axe  (mesurée  en  hauteur  du  fluide),  comprend  donc  quatre  parties  :  i°  celle 
qui  provient  du  terme  en  b  ou  qu'absorbe  le  frottement  ordinaire  de  régime 
uniforme;   20  une  autre,  positive  comme  la  précédente  et  également  notable, 

—  /  d(ot  —  1  —  r,)  ou  sensiblement  2  —  j  a\  (vu  a  =■  1  -f-  3r,  à  peu  près), 

employée  à  accroître  les  inégalités  de  vitesse  des  fdets  fluides  et  simple- 
ment proportionnelle  à  l'augmentation  du  coefficient  a.  —  1  —  r,  entre  les 
deux  sections  considérées;  3°  et  4°>  enfin,  deux  petites  parties, 

(91)  -  /  3îl(—  T3u')dx,         —(B0D\L/)^U2JTsldx, 

du  second  ordre  de  petitesse  comme  les  produits — zzu',  cj^,  sauf  sur  une 
faible  longueur  près  de  l'entrée,  où  la  fonction  zs  est  comparable  à  9  —  1, 


(  53  ) 
l'inégalité  des  vitesses  n'y  étant  encore  qu'ébauchée.  La  dernière  (91)  est 
évidemment  négative,  comme  —  gt*.  Quanta  la  précédente,  elle  est  posi- 
tive. En  effet,  d'une  part,  les  filets  périphériques,  d'abord  trop  rapides, 
ou  pour  lesquels  o  est  positif,  se  ralentissent  et  ont  leur  accélération  iï 
négative,  tandis  que,  d'autre  part,  les  filets  voisins  de  l'axe  et  pour  lesquels 
cr  est  négatif,  s'accélèrent  :  le  produit  — vsu!  est  donc  positif  dans  presque 
toute  la  section  et  a  sa  valeur  moyenne  positive. 

»  Si,  pour  prendre  le  cas  le  plus  simple,  on  suppose  l'entrée  du  tuyau 
assez  bien  évasée  pour  que  les  filets  fluides  soient  sensiblement  parallèles 
dès  l'origine  de  sa  partie  prismatique  ou  cylindrique,  la  formule  deD.  Ber- 
noulli  leur  atltribuera  à  cet  endroit,  comme  on  sait,  la  vitesse  commune, 
U,  due  à  la  hauteur  motrice  dès  lors  dépensée  à  partir  des  points  du  réser- 
voir d'admission  où  le  fluide  est  en  repos.  On  y  aura  donc  a  —  1  —  yi  =  o. 
Et,  par  suite,  la  hauteur  motrice  totale  dépensée ,  depuis  le  réservoir  jusqu'aux 
points  où  régnera  le  régime  uniforme,  pour  établir  ce  régime,  ou  en  sus  de 
ce  qu'y  absorbera  le  frottement  ordinaire  de  régime  uniforme,  sera 

(92)  ~   ^-K»-1  — *>)-*-uï/  3it(—  vsu')dx-  gV>0(d}if)~j  f  vs\dx 

»  Les  coefficients  1  +  7),  a  y  désignent  les  valeurs  moyennes  de  <p2  et  o3; 
autrement  dit,  ils  se  rapportent  à  la  limite  supérieure  des  intégrations,  ou 
aux  sections  g  dans  lesquelles  le  régime  uniforme  existe.  L'abscisse  x  de 
celles-ci,  comptée  à  partir  de  l'origine  du  tuyau  prismatique  ou  cylindrique, 
peut  d'ailleurs  être  supposée  infinie,  les  fonctions  sous  les  signes  J  de  (92) 
y  tendant  asymptotiquement  et  assez  rapidement  vers  zéro. 

§  XXIII.  —  Équations  qui  y  régissent  le  mode  de  distribution  des  vitesses. 

»  46.  Abstraction  faite  des  éléments  les  plus  voisins  de  la  limites?  =  o, 
d'une  somme  probablement  insignifiante,  les  deux  intégrales  que  contient 
l'expression  (92)  s'évalueront,  avec  une  assez  faible  erreur  relative,  en 
supposant  la  fonction  gj  de  l'ordre  des  petites  quantités  dont  nous  négli- 
geons habituellement  les  produits.  C'est  donc  dans  cette  hypothèse  simpli- 
ficatrice qu'il  nous  reste  à  déterminer  a,  ou,  ce  qui  revient  au  môme,  F2. 

»   Nous  aurons  pour  cela  le  système  (10)  d'équations,  dans  lequel  u'  sera 
donné  par  la  formule  (4°)»  évidemment  réduite  à 
/on         /       tt/     do  do\         TT„     dm         r,a/    v     do  w    do  dm\ 

(93)  u=v{vjy  +«'aï)  +  u,?E  =  D  (ïoaï  +  ïU3?  +  «,s/ 


(  54  ) 

»   D'ailleurs,  d'après  (46),  où  y0,  z0,  a,  h  seront  constants  et  la  valeur 

(02)  de  /,  nulle,  les  vitesses  transversales^,  sauront  simplement  "X  et  [j.  pour 

quotients  par  aU  et  h\j,  avec  X,  u.  régis  par  l'équation   indéfinie  (53)  et  la 

condition  au  contour  (49)»  sans  compter  la  condition  d'intégrabilité  (1). 

»  47.  Bornons-nous  aux  deux  cas  de  la  section  rectangulaire  très  large, 
de  hauteur  ih,  et  de  la  section  circulaire  de  rayon  R,  où  nous  savons  que, 
par  raison  de  symétrie,  v,  w  ou  X,  y.,  fonctions  impaires  de  y,  z,  sont  les 
deux  dérivées  en  y,  z  ou  en  r,,  £>  d'une  môme  fonctionnaire, soit  de  '(,  soit 
de  r\  et  £  par  l'intermédiaire  du  rapport,  x,  au  rayon  R  =  ia  =  ih,  de  la 
distance  r  =  \y2-\-z-  à  l'axe  :  ce  qui  rend  identique  la  vérification  de  la 
condition  (1)  d'intégrabilité. 

»  Nous  introduirons  comme  inconnue  auxiliaire  une  fonction  co  de  '( 
ou  x,  et  de  x,  dont  la  dérivée  en  x  soit  justement  cetle  fonction  qui  a  X 
et  y.  pour  dérivées  respectives    en  n  et  '(.  Autrement  dit,  nous  poserons 

(94)  x  = 


dx  dr.  { '        dx  dl 


co  comprendra  ainsi  une  fonction  arbitraire  de  r,  et  £.  Nous  en  disposerons 
de  manière  que,  sur  une  première  section  c,  celle  qui  aura,  par  exemple, 
l'abscisse  x  =  o,  et  où  la  valeur  de  a,  évidemment  paire  en  r,,  '(,  sera 
directement  donnée,  l'on  ait 

/    t  \  a  d-w        d2u  dw  ,         ,  N 

(  od)  Aoco-hcr  =  o  ou  — ;-- 4- -^  =  —  nr,  avec  -,  =  o  (  sur  le  contour). 
v  *'    7      "  dr,  -         tf;2  r/v  v  J 

»   Alors,  vu  les  expressions  (94)  de  X  et  ;;.,  l'équation  indéfinie  (53)  et 
la  condition  (/jo,)  deviendront  respectivement  -7-  (A2 a>-t-Gj)=o,  —  ( -^)  =0. 

Cl  Ju  Cl  tic    \Clyt   t 

Donc,  les  équations  (93)  seront  vérifiées  pour  toutes  les  valeurs  dex,  c'est- 
à-dire  sur  toutes  les  sections  a.  On  voit  que,  si  l'on  connaissait  sur  l'une 
d'elles  quelconque  la  fonction  paire  ci  de  r,  et  de  C,  ces  équations  (95)  v 
détermineraient,  à  une  constante  arbitraire  prèsf(x),  la  fonction  co,  égale- 
ment paire  en  r,,  'C.  La  partie/"^)  de  co,  évidemment  étrangère  aux  rela- 
tions (94)»  (95)  de  co  avec  nos  vraies  inconnues  X,  >x,  u,  reste  indéter- 
minée. 

»   Enfin,  l'expression  (g3)  de  u',  en  en  éliminant,  par  (o,5)  et  (94),  o  et 
les  rapports  X,  ;x  de  v,  w  à  a\J,h\J,  puis  observant  que  o  ne  dépend  pas 
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de  x,  et  substituant  enfin  à  9  sa  valeur  tirée  de  (5  i),  deviendra 

,       TTO     d   (do  du>         do  du  . 

dx  \drt  drt         dZ,   dZ,  ' 

AVBJJ'         rf 
1  H-  Ary/BoOl^F,  ^ 


/.__!__   1    p  \  A  ^F,  <^w        <^F,  dia~\ 

A*^  /  "  ~^~  dr<   '     ~~dT~dZ.  I 


»  48.  Telle  est  la  valeur  de  u'  qu'il  faudra  porter  dans  l'équation 
indéfinie  en  F2,  qui  est  la  première  (10)  ou  mieux  la  différentielle  totale 
en  y),  £  de  la  première  (10).  En  effet,  avec  l'adjonction  de  la  condition  au 
contour  correspondante,  la  première  équation  (10)  équivaut  exactement  à 
sa  différentielle  totale  en  yj,  Ç,  que  nous  pourrons  écrire,  avec  un  terme  de 
moins, 

(97)         4l(Ff)+l(Ft)-"'H' 

l'indicée  indiquant  que  la  différentiation  dont  il  s'agit  se  fait  sans  sortir 
d'une  même  section  g.  Car  l'équation  (97)  revient  évidemment  à  la  pre- 
mière (10),  avec  addition,  au  second  membre,  d'une  constante  arbitraire  C. 
Or,  si  l'on  prend  la  moyenne  des  valeurs,  dans  l'aire  n,  de  tous  les  termes 
de  la  première  (10)  ainsi  complétée,  on  trouve  zéro  pour  valeur  moyenne 
du  premier  membre;  car  celui-ci,  multiplié  par  dr,  d^,  et  intégré,  se  change 
évidemment  en  une  intégrale  de  contour  où  figure  sous  le  signe  f  le  pro- 
duit de  F  par  la  dérivée  en  v  de  F2,  nul  d'après  la  seconde  relation  (10). 
Et  il  vient  bien  ainsi,  forcément,  C  =  o. 

»  49.  Cela  posé,  tirons  de  (1  1),  pour  la  substituer  dans  (97),  la  valeur 
de  F2  en  fonction  de  Ft  et  de  gj  ou  w.  A  cet  effet,  observons,  d'une  part, 
que  dans  l'expression  de  F2  fournie  immédiatement  par  (1 1),  le  facteur  li1 
peut  être  remplacé  à  très  peu  près  par  U2<pj);  d'autre  part,  que  le  quotient 
de  a  parw0  est  identiquement  celui  de  9  -hw  par  <p0  -+-gj0  et  excède  sa  valeur 

de  régime  uniforme  i-h  k\/B0¥t,  ou  —  .  de 


Il  viendra 


(98) 


~(oa^  —  o„<p)  =  ~  -  ^(i  +  tyfW,¥,). 


Fs  =  ^U2  U' [*  -  n„  ( .  +  *VB"„  F,  )] 

"  K   yw'ir?r'p^    [*,<■> +  «.d  +  *V'B;F,)]. 
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Les  termes  en  F,  de  celle-ci  s'élimineront  de  la  relation  (97),  à  raison  de 
ce  que  la  première  équation  (9)  y  réduira  leur  somme  à  la  différentielle  da 
d'une  constante,  et,  en  même  temps,  la  substitution  à  iï ,  dans  (97),  de  la 
dernière  expression  (96),  donnera  l'équation  indéfinie  en  u>, 


/     1  „  \  é/a,w        d  fd¥x  du>       dFi  d< 

\k  \fïT0  1  /    dx         dx  \  drt    dr{        <r/Ç    d 


=  0. 


»  La  quantité  entre  accolades  ne  dépendant,  comme  on  verra,  que  de 
x  et  de  £  ou  t,  la  différentielle  da  équivaudra  à  une  simple  dérivation  en  '( 
ou  x,  et  l'équation  (99)  sera  une  équation  aux  dérivées  partielles,  du  cin- 
quième ordre  en  '(  ou  en  t. 

)>  Il  faudra  y  joindre  la  condition  (10)  au  contour,  exprimant  que  s'y 
annule  le  produit  de  F  par  la  dérivée  en  v  de  F2,  c'est-à-dire  de  la  dernière 
expression  (98).  Si  l'on  tient  compte  de  la  seconde  relation  (9)  et  de  ce 
que  G70  sera  justement  la  valeur  de  cr  ou  de  —  A2to  au  contour,  il  viendra, 
en  rappelant  d'ailleurs  une  autre  condition  au  contour  déjà  établie  plus 
haut  pour  co  (formules  93), 


(100)      (au  contour  des  sections)     — = j h  A2co  =  o, 


d\=,  co  ,  du. 


=  o. 


»  D'ailleurs,  la  dernière  relation  (10)  est  satisfaite  identiquement  par 
la  valeur  (98)  de  F2.  Et  la  condition  fvsdn  =  o,  provenant  de  ce  que  s'an- 
nule dans  chaque  section  la  valeur  moyenne  de  gt,  différence  de  celles  de 
o  -+-  ct  et  de  o,  égales  toutes  les  deux  à  l'unité,  n'est  pas  moins  vérifiée  ;  car 
l'égalité  —  tô  =  Aoto,  multipliée  par  dr.dX.,  puis  intégrée  dans  toute  l'aire  a, 
conduit  à  une  intégrale  de  contour  où  la  fonction  sous  le  signe  f  est  la  dé- 
rivée de  co  en  v,  nulle  en  vertu  de  la  seconde  relation  définie  (100). 

»  Si  l'on  veut  ne  considérer,  dans  chaque  section,  que  la  variable  indé- 
pendante Ç2  ou  t,  croissante  de  o  à  1,  il  y  aura  lieu  d'observer  qu'au  centre 
Çf  =  o  ou  t  =  o),  les  deux  fonctions  paires  (en  r,  et  £)  co  et  —  gt  =  A2co, 
maxima  ou  minima  en  ce  point,  ont  leur  différentielle  da  nulle  par  raison 
de  symétrie  et  de  continuité.  On  aura  donc  encore 

(101)  (au  centre  des  sections)     dGi»  =  o,         daA.2u  =  o. 
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§  XXIV.  —  Décomposition  de  ce  mode,  dans  sa  partie  amortissable  ou  de  non-uni- 
formité, en  modes  simples,  plus  ou  moins  lents  à  s'évanouir  aux  distances  crois- 
santes de  l'entrée. 


»  50.  On  vérifiera  les  équations  (99),  (100),  (101),  en  prenant  pour  cd 
une  somme  de  solutions  simples,  dont  chacune  sera  le  produit  d'une  con- 
stante arbitraire  c  par  une  exponentielle  décroissante  en  x,  où  figurera 
un  coefficient  positif  d'extinction  m  constant,  et  par  une  fonction  12  tic  '(  ou 
de  x  seuls.  Bref,  on  posera 

(102)  on  =  2dCe 


,  I  'L 

i-  1 


»  Chaque  terme  de  la  somme  1  devant  séparément  satisfaire  à  (99). 
(100)  et  (toi),  il  vient  en  £2  l'équation  différentielle  linéaire,  du  cinquième 
ordre  et  sans  second  membre, 


(io3)     <IG 


d  f     dl2il 
dr\  \         dr. 


d  /„  r/A,<> 

dz[l^r 


A,<>  dF,  d" 

2       _i_  F    \  o  _  __i  _ 

—  t^  i  1  -il*"  /       1 

saVb0  «*1   «1 


r/F,  dil 
~dl    dC 


=  o. 


Et  l'on  devra,  d'une  part,  adopter  sa  solution  particulière,  déterminée 
(dans  sa  partie  variable)  à  un  facteur  constant  près,  qui  vérifiera,  d'après 
(101)  et  (100),  les  trois  conditions  spéciales 


dQ 


(io/f)     (au  centre)  da£l  =  o     et     r/cA,12  =  o,      (au  contour)  -7-  =  o; 

d'autre  part,  choisir  pour  m  les  racines,   tenues  d'être  toutes  positives,  de 
l'équation  transcendante  fournie  par  la  première  condition  (100), 


(,o5) 


d\,il 


(sur  le  contour)     — = ~ 

J      /.\15«/    ch 


4-  A,£2  =  o. 


»   Comme  c'est  la  partie  u  =  —  A2w,  variable  avec  x,  du  mode  de  dis- 
tribution des  vitesses  que  l'on  désire  connaître,  la  formule  (102)  donnera 


(106) 


.  =  -2 


ce 


»    SI.   Il  restera  à  calculer  les  coefficients  c  de  manière  que,  sur  une  pre- 
mière section  a,  celle  qui  a,  par  exemple,  l'abscisse  x  =  o,  gj  reçoive  cer- 
taines valeurs  initiales  zsh  fonction  de  *(2  ou  de?,  données  dans  toute  l'aire <s. 
La  formule  (10G)  devra  donc  y  devenir  ci,  =  —  2cA2£2.  Mais,  ne  pouvant 
B.  8 
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effectivement  prendre,  dans  la  somme  1  ordonnée  suivant  les  racines  m 
croissantes,  qu'un  nombre  fini  n  de  termes,  on  devra  tolérer  des  erreurs 
rr,  -t-  ïc±.2<2;  et  l'on  rendra  seulement  minima  la  moyenne  de  leurs  carrés 
dans  toute  l'aire  c  ou,  ce  qui  revient  au  même,   l'intégrale 

( 1 06  bis)  f(vti  +  Zc\.2 P.)2  ch. 

»  J'ai  montré  sur  un  exemple,  dans  une  Note  relative  au  même  problème 
de  l'établissement  du  régime  uniforme,  mais  à  l'entrée  des  tubes  fins,  com- 
ment cette  condition  détermine  les  coefficients  (  '  ). 

»  52.  On  voit  que,  l'abscisse  agrandissant,  la  formule  (10G)  fait  évanouir 
successivement  tous  les  termes  de  la  somme  1,  à  commencer  par  les  plus 
éloignés.  Aux  distances  de  l'entrée  excédant  un  nombre  médiocre  de  fois 
le  rayon  moyen,  il  ne  subsistera  donc  que  le  premier  ternie,  affecté  de  la 
plus  petite  racine  positive  m  de  l'équation  (ro.5);  et  u,  d'ailleurs  indépen- 
dant, d'après  (106),  de  la  vitesse  moyenne  U  et  des  dimensions  absolues 
de  la  masse  fluide,  variera  désormais,  sur  chaque  section,  proportionnel- 
lement à  l'expression  correspondante  de  A2!£. 


(')  Comptes  rendus  des  G  et  1 3  juillet  1 8g  1 ,  t.  <  XIII,  p.  9  et  4(J-  La  condition  dont 
il  s'agit,  qui  revient  à  rendre  minima  la  somme  des  carrés  des  erreurs  commises  sur  la 
valeur  de  rn,-  dans  les  divers  éléments  (équivalents)  de  Taire  7,  entraîne  évidemment 
l'annulation  de  la  dérivée  partielle  première  de  l'intégrale  (10G  bis)  par  rapport  à  cha- 
cun des  coefficients  c  à  déterminer.  Si  donc  \.,Llj  désigne  successivement  chacune 
des  n  expressions  obtenues  pour  A2£2,  c'est-à-dire  celles  qui  correspondent  aux  n  pre- 
mières racines  m  de  l'équation  transcendante  (io5),  il  vient,  pour  calculer  les  n  coef- 
ficients c,  le  système  complet  du  premier  degré 

(106  ter)  scf\,ii^,iijch—f{—vyi)iio.jd^. 

Contrairement  à  ce  qui  arrive  dans  les  problèmes  les  plus  usuels  de  la  Physique 
mathématique,  les  inconnues  c  ne  seront  pas  séparées  clans  ce  système,  car  le  produit 
de  deux  expressions  différentes  de  A2<>  n'aura  pas  sa  valeur  moyenne  nulle.  Autrement 
dit,  les  équations  (106  ter)  ne  se  réduiront  pas  à 

Cjf(L%ùj)*d<i=f{-wi)ltQjd*, 

et  les  valeurs  des  constantes  c  varieront,  dans  une  certaine  mesure,  avec  le  nombre  de 
celles  que  l'on  voudra  utiliser  pour  le  calcul  effectif  de  rs. 
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§  XXV.  —  Cas  d'un  tuyau  à  section  rectangulaire  large  :  intégration  en  série. 

»  o5.  Abordons  d'abord  le  plus  simple  des  deux  cas  en  vue  desquels  a 
été  établie  spécialement  la  théorie  précédente  et  où  l'on  af=  i,  B0  =  B, 
savoir,  celui  d'un  tuyau  à  section  rectangulaire  très  large,  de  hauteur  ih. 
Alors  F  =  i,  F,  =  ^(i  —  YJ)  et  1>,  fonction  paire  comme  F,,  ne  dépend 
également  que  de  Ç  :  il  suffira  de  la  considérer  depuis  l'axe  jusqu'au  fond, 
c'est-à-dire  de  ^  =  oà(=i. 

»  L'équation  indéfinie  (io3)  ne  contient  il  que  par  sa  dérivée  pre- 
mière £',  car  A2Î2  est  maintenant  la  dérivée,  £2",  de  celle-ci.  Il  y  a  donc 
lieu  d'adopter  pour  inconnue  la  fonction  impaire  Q',  que  nous  appellerons 
dès  lors  W,  et  dont  A, £2  sera  la  dérivée  W.  D'ailleurs,  le  signe  drj  équivalant 
à  une  différentiation  en  'C,  l'équation  différentielle  (io3)  sera  du  quatrième 
ordre  en  1I*.  Dans  son  terme  en  W",  évaluons  et  mettons  à  part  la  valeur 
moyenne  du  coefficient,  valeur  qui  égalera  l'inverse  de  k\jb,  d'après  la 
première  formule  (37)  de  mon  Étude  de  l'année  dernière  (').  Cette  équa- 
tion en  T,  divisée  par  m,  sera  dès  lors 


(107)  -TIV  + 


1  !   /  1 


T  =  o  ; 


et  les  conditions  définies  (ro/J),  (io5),  où  l'on  aura  a\  =  <7'C  à  la  limite  Ç  =  1 , 
deviendront 

(108)     *F(u)  =  o,        W(ï)  =  o,        T"(o)  =  o,        W"(i)  +  ks!Ïm"(\)  =  o. 

»  Les  trois  premières  (108),  jointes  à  ([07),  suffisent  évidemment  pour 
définir  la  fonction  W,  à  un  facteur  constant  près,  que  l'on  peut  choisir  à 
volonté  puisque  le  résultat  cherché  (106)  contient  déjà  les  coefficients 
arbitraires  c.  Nous  disposerons,  en  général,  de  ce  facteur  constant,  de 
manière  à  avoir  ^"'(o)  =  1  ;  ce  qui  achèvera  de  déterminer  xl\ 

»  5i.  Si  nous  adoptons  pour  cette  fonction  impaire  une  série  procédant 
suivant  les  puissances  entières  et  positives  de  £,  nous  aurons  donc  une  for- 


(')   Comptes  rendus,  t.  CXXI1I,  p.  12  (6  juillet  1896),  ou  p.  33  du  Mémoire  précé- 
dent. 


(6o) 

mule  comme 

(109)  W(K)  =  K  +  A?  -  C'C5  -h  T>V  -Et,9  ■+■..., 

et  les  première  et  troisième  des  conditions  (108)  seront  satisfaites.  Quant 
à  l'équation  (107),  que  le  développement  (109)  devra  identiquement  véri- 
fier, on  reconnaît  qu'elle  laisse  disponible  le  coefficient  A,  mais  qu'elle 
détermine  chacun  des  suivants  C,  D,  E,  ...,  en  fonction  linéaire  très 
simple  des  âeux  qui  le  précèdent,  1,  A,  ...,  multipliés  par  m.  Puis  A  se 
détermine  parla  seconde  relation  (108);  et  la  dernière  (108)  donne  enfin 
l'équation  en  m  ('  ). 


(')  Sur  rétablissement  du  régime  uniforme  le  long  d'un  tube  fin  à  section  rec- 
tangulaire relativement  large.  —  La  loi  de  récurrence  qui,  à  partir  de  C,  relie  chaque 
coefficient  de  la  série  (10g)  aux  deux  précédents,  se  simplifie  un  peu  quand  on  sup- 
pose B  très  grand  et,  par  suite,  k\Jb  égal  à  3.  La  quatrième  condition  (108),  résolue 
par  rapport  à  lr'(i),  se  réduit  en  même  temps  à  W  (1)  =  o.  Alors,  d'après  les  équa- 
tions (8)  ou  (12),  la  vitesse  u0  a  la  paroi  devient  négligeable  à  cùlé  de  la  vitesse 
moyenne  U,  et  le  tuyau,  supposé  ainsi  infiniment  rugueux,  immobilise  par  son  frotte- 
ment la  couche  fluide  qui  le  touche,  tout  comme  le  ferait  par  adhérence  la  paroi  polie 
d'un  tube  fin,  dans  un  écoulement  bien  continu.  D'ailleurs,  celte  vitesse  «0  étant  sen- 
siblement indépendante  de  x  aux  endroits  considérés  où  ttt  est  très  petit,  le  coefficient 

de  frottement  intérieur  e  =  -—  \/Bhu0  s'y   trouve  à  fort  peu  près  constant,  comme  il 

l'est  dans  les  mouvements  bien  continus.  Donc,  tout  se  passe,  quant  aux   équations 

déterminant  m   et  W,  comme  s'il  s'agissait  de  l'écoulement  effectué,  le  long  d'un  tube 

fin,  par  un  fluide  où  le  coefficient  de  frottement  intérieur  aurait  justement  la  valeur 

P  g    1 — 
numérique  de  l'expression  *-y-  y  13  h //0  dans  notre  liquide. 

Ainsi,  la  formule  (106)  de  m,  spécifiée  pour  un  tuyau  infiniment  rugueux,  sera 
celle  qui  représentera  l'établissement  du  régime  uniforme  à  la  partie  amont  d'un   tube 

fin,  si,  remplaçant  -  par  T  et  A2<2  par  W,  l'on  exprime  d'ailleurs,  dans  les  exponen- 

tielles,  k*  en  fonction  de  z,  h  et  U.  L'on  aura,  pour  cela,  l'équation  (3i)  de  mon  pre- 
mier Mémoire  (p.  27),  équation  qui,  vu  la  valeur  \  de  Ole  F,  et  la  grandeur  supposée 

de  \l\i,  donnera  k\/]iu0  =  3U,  dans  l'expression  —■  kyBu0h  de  e.  De  celle-ci,  devenue 

dès  lors  z  ■=  — ^7-; — >  on  tirera  —  =  - —  ■  Donc,  la  formule  représentant  la  partie 

variable  avec  x  du  mode  de  distribution  des  vitesses,  dans  le  problème  de  l'établisse- 
ment du  régime  uniforme  Je  long  d'un  tube  fin  à  section  rectangulaire  large  et  de  hau- 
teur 2  h,  sera 

m  £  .r 


(6,   ) 


§  XXVI.  —  Solutions  simples  en  termes  finis,  quand  les  parois  sont  lisses. 

»  oo.  Les  résultats  deviennent  très  simples  quand  on  suppose  la  paroi 
assez  lisse,  ou  plutôt  b  assez  petit,  pour  que,  dans  (107),  le  coefficient  de  W", 
alors  très  grand,  puisse  être  remplacé  par  l'inverse  de  k\Jb,  sa  valeur 
movenne,  dont  il  ne  s'écarte  qu'entre  les  limites  fixes  |  et — |  comptées 
respectivement  au  delà  et  en  deçà.  Le  troisième  terme,  W,  de  (107)  se 
trouve  dès  lors  encore  plus  négligeable  que  la  petite  partie  supprimée  du 
terme  en  W";  car  W,  s'annulant  aux  deux  limites  £  =  o,  £  =  1,  est,  en 
moyenne,  dans  l'intervalle,  au  plus  de  l'ordre  de  grandeur  de  sa  dérivée 
première  W,  qui,  tenue  elle-même  de  s'y  annuler  une  fois,  y  esta  peine 
comparable  à  sa  dérivée  W".  Donc,  à  plus  forte  raison,  W  disparaîl, 
dans  (107),  devant  la  partie  principale  du  terme  en  W".  Et  l'équation  (107), 
dès  lors  binôme,  donne,  à  un  facteur  constant  près,  vu  l'avant-dernière 
condition  (108),  W"  =  —  xsinaÇ,  si  x2  (avec  x  pris  positif)  désigne  le  quo- 
tient de  m  par  k\[b.  Or,  de  cette  valeur  de  W",  multipliée  deux  fois  par  d'Ç 
et  inlégrée  chaque  fois,  il  résulte,  en  tenant  compte  des  deux  premières 
conditions  (108), 


(  1 1  o)  *L 


ou 


\     ks/l 


»   Enfin,  la  dernière  condition  (108)  devient  l'équation  en  m  ou  en  x, 
tang  x  =  — —  — •  Supposons  le  quotient  positif  de  k\Jii  par  x  assez 

a       r+/î  \/B 

petit  pour  qu'on  puisse  négliger  son  carré;  en  sorte  que  ce  quotient  soit 
lui-même  voisin  de  zéro  et  x  voisin  d'un  multiple  positify-  de-.  Le  petit 

.  ,      ,  •  ,  ,  .       ,    .   ,  ,      A-y/B        ks/B. 

excédent  x — j-  aura  donc  pom    tangente,  a  très  peu  près,  — - — ou — ? — , 

a  j  .. 

où  m  et  ir'  se  calculeront  comme  s'il  s'agissait  d'un  tuyau  de  même  forme  à  paroi  infi- 
niment rugueuse.  Quant  à  la  partie,  tp,  de  régime  uniforme  ou  indépendante  de.r,  elle 
sera,  en  faisant  y/B0  infini  dans  (5i), 

F,  3 


3ILF, 


(i-O- 


Malgré    les    simplifications    indiquées   des   formules  {107)    et    (108),   la    série   qui 
exprime  W  et  l'équation  en  m  restent  encore  d'un  calcul  fort  laborieux. 


et  il  viendra  successivement 


(ni) 


(- 


m 


(62) 


*>lfî  =  (j>**+2il/B)tjb, 


W  =  cosx'C 


=  cosj-'Ç 


/.yB  /cos/' 


(sinyV£ 


m  Les  valeurs  de  m  sont  bien,  comme  il  le  fallait,  toutes  réelles  et  posi- 
tives; celles  de1!",  qu'on  substituera  à  A.2Q,  s'écartent  peu  de  cos  jrX, 
dont  elles  diffèrent  d'une  même  quantité  aux  deux  limites  "£2  =  o,  Ç2  =  i. 

»  56.  Bornons-nous  au  terme  de  u  fondamental  ou  affecté  de  l'exponen- 
tielle la  plus  lente  à  s'évanouir,  celle  qui  correspond  à/  =  i.  Admettons 
de  plus  une  entrée  du  tuyau  parfaitement  bien  évasée,  qui  y  donnerait 

u  =  U  ou  crJ-=  I  — ip=-i^K2-J)  d'après  les  formules  (37),  précé- 
demment citées,  de  mon  Etude  de  l'année  dernière.  Le  coefficient  c,  devant 
rendre  minimum  le  carré  moyen,  /  (ccos-Ç  -+-  vs^dZ,  environ,  de  l'éeart 
entre  cW  et  —  uif  se  déterminera  par  la  condition  approchée 

/    (ccosr'C  4-  gj,-)  COS-'(r/'C  =  o. 

Il  sera  sensiblement  le  quotient,  par  la  valeur  moyenne  -  de  cos2t:'C,  de 
l'intégrale   /    ( —  ct/)cosx£</£,  que  donnera  presque  immédiatement  l'inté- 

gration  par  parties  effectuée  avec  sin-'C  ou  cos-'C  comme  facteur  intégré. 
Et  il  viendra  finalement,  vu  (106), 


(112) 


ik\fb   - 

tj  =  — :—  e 


y  k  \/B    ,  y       .  y  , 

COStcC  H zr~  (f  —  "C  sin-Q 


»   Les  vitesses  u  sont  inférieures  à  celles  de  régime  uniforme  dans  la 
partie  centrale   des  sections  qui  s'étend  jusqu'à   l'ordonnée  relative,   '(, 

dont  la  valeur  absolue, -—  ( j,  annule  W;  elles  leur  sont,  au 

contraire,  supérieures  dans  tout  le  reste,  qui  comprend  un  peu  plus  que  la 
moitié  des  sections.  Aussi,  l'écart  —  cr  est-il  plus  fort  au  centre  Ç  =  o  qu'à 
la  paroi  £  =  1,  comme  il  était  naturel  de  le  penser,  l'action  régulatrice  du 
frottement  devant  le  plus  vite  se  faire  sentir  au  voisinage  des  parois. 


(  63  ) 

§  XXVII.  —  Longueur  nécessaire  pour  l'établissement  approché 
du  régime  uniforme  dans  un  tel  tuyau. 

»  57.  Pour  la  paroi  en  ciment  fin  sur  laquelle  M.  Bazin  a  fait  ses 
récentes  observations  (dans  un  tuyau  circulaire,  il  est  vrai,  et  non  rectan- 
gulaire large),  on  avait  b  =  0,000166,  \[b  —  0,0129,  et,  par  suite  (avec  la 
valeur  48,6  de k),  y/B  =0,0175,  k  VB  =o,85i  (')  ;  d'où  résulte,  dans  la  sec- 
tion rectangulaire  large,  encore  d'après  les  mêmes  formules  (32)  de  l'Étude 
citée,  &  =  0,000186,  \7>  =  o,oi36  et  k\fb  =  o,663.  L'inverse  de  k\fb  y 
est  donc,  non  pas  très  grand,  mais  seulement  égal  à  i,5o8.  Toutefois,  cha- 
cune des  deux  catégories  de  valeurs  que  prend  le  coefficient  de  W  dans 
(107),  les  unes  (pour  £2  <î)  en  excédent,  les  autres  ( pour  'Ç2  >|)  en  déficit 
sur  cette  moyenne  i,5o8,  s'en  écarte  relativement  assez  peu,  surtout  en 
movenne,  pour  qu'on  puisse,  dans  une  première  étude,  lui  substituer  la 
valeur  constante  î,jo8;  et  l'expression  (110)  de  "F,  où  a.  est  de  l'ordre 
dey-,  réduit  encore  le  terme  W,  dans  (107),  à  une  fraction  presque  négli- 
geable du  précédent  i.SoSW,  au  yr  environ.  Enfin,  dans  l'expression  (1 1 1) 
de  m,  le  terme  ik\  l>,  de  deuxième  approximation,  n'est  guère,  même 
poury  =  1,  que  le  £  du  terme  principal/2-2  :  preuve  que  les  formules  obte- 
nues continuent  à  être  applicables  avec  quelque  approximation. 

(')  Ces  valeurs  se  déduisent  de  la  première  formule  (60)  de  mon  travail  Je  l'année 
dernière  (voir  le  Mémoire  précédent,  p.  46),  où  k  =  48,60  d'après  Tune  des  for- 
mules (07)  du  même  travail,  et  où  il  suffit  d'introduire  en  outre  le  résultat 
&  =  0,000166,  fourni  directement  par  l'observation  des  débits  du  tuyau  circulaire 
expérimenté. 

Si,  dans  l'étude  du  régime  uniforme  à  l'intérieur  des  tuyaux  circulaires,  l'on  se 
bornait,  comme  je  le  fais  ici  dans  l'étude  des  régimes  variés,  à  cette  première  approxi- 
mation où  le  coefficient  de  frottement  intérieure  est  supposé  valoir  celui  du  cas  de  la 
section  rectangulaire  large  pour  même  rayon  moyen  et  même  vitesse  à  la  paroi,  mul- 
tiplié par  l'inverse  de  la  distance  relative  t  à  l'axe,  la  valeur  de  À  la  plus  propre  à 
réduire  autant  que  possible  les  écarts  entre  la  théorie  et  les  récentes  observations 
de  M.  Bazin  serait  encore,  d'après  les  pages  \!\  et  45  du  Mémoire  précédent, 
A  =  48,60.  Et  les  premières  formules  (43),  (37)  du  même  Mémoire  donneraient, 
l'une  (en  y  faisant  b  =  0,000166)  y^B  =  0,0172,  l'autre,  ensuite,  \lb  —  o,oi35  pour 
la  section  rectangulaire  large.  Ce  seraient  donc,  à  fort  peu  près,  les  valeurs  de  y'B,  sjb 
adoptées  dans  le  texte;  et  l'on  aurait  aussi,  par  suite,  sensiblement  les  mêmes  valeurs 
que  ci-dessus  pour  k\/B  et  k\^b. 


(64) 

»  L'expression  (i  12)  de  —  u  donnera  donc  une  idée  encore  assez  juste 
du  phénomène  étudié.  Le  coefficient  c  figurant  devant  l'exponentielle  y 
est  environ  o,i34,  vu  la  valeur  o,6G3  de  k\Jb.  Mais,  pour  un  tuyau  non 
muni  de  la  bouche  parfaitement  évasée  que  nous  avons  admise,  et  où  se 
produira  toujours,  après  la  brusque  contraction  des  filets  fluides,  un  épa- 
nouissement rapide,  avec  frottements  notables  qui  ébaucheront  déjà  l'iné- 
galité des  vitesses  dans  le  tuyau,  l'écart  initial  rzh  sur  la  section  où  le 
régime  commencera  à  varier  graduellement,  aura  des  valeurs  absolues 
moindres  que  leur  expression  supposée  1  — 9,  et,  par  suite,  le  coefficient  c 
ne  devra  guère,  ou  pas,  excéder  0,1.  Il  suffira  donc  que  l'exponentielle  se 
réduise  elle-même  à  0,1,  ou  que  son  exposant  égale  au  moins,  en  valeur 
absolue,  2,3026,  ou  enfin  que  a?  atteigne  la  valeur  72, 3h  environ,  pour  que 
l'écart  cj  soit  partout  inférieur  à  -^  et  insensible.  Ainsi,  le  régime  uniforme 
sera  établi  après  un  parcours  x  d'environ  72  rayons  moyens,  ou  3G  fois  la  hau- 
teur 2  h  de  la  section,  à  partir  de  V  endroit  où  les  filets  fluides  commencent  à 
être  sensiblement  rectilignes  et  parallèles. 

§  XXVIII.  —  Cas  d'un  tuyau  circulaire  :  intégration  en  série. 

»  58.  Abordons  enfin  le  cas,  plus  pratique,  mais  beaucoup  moins 
simple,  d'un  tuyau  à  section  circulaire.  Alors  la  fonction  F  peut  être  réduite 
à  l'inverse  de  t  avec  une  approximation  suffisante  :  F,  est,  par  suite, 
|(i  —  t3),  et,  Q  ne  dépendant,  comme  F  et  F,,  que  de  la  variable  ■<  =  Vr'2  + 


^  » 


1    d.til 


&.,Q  a  l'expression  7 '-r—  D'ailleurs,  dans  l'équation  indéfinie  (io3),  les 


4^ 


deux  derniers  termes  du  quadrinome  entre  parenthèses  ont  évidemment 
pour  somme  —  }F',  Œ',  ou  ^x2Q';  de  sorte  que  cette  équation  indéfinie  ne 
contient  £2  que  par  le  produit  tS2'.  Celui-ci,  ou  plutôt  son  quart  \iQ! ',  sera 
donc  notre  inconnue  auxiliaire.  Nous  l'appellerons  encore  W,  en  posant 
ainsi 

/      q\  a    ^        '  dxv        w' 

(ii3)  A2_Q  =  --_  =  _, 

x         y  x   dx  x 

»  Enfin,  le  signe  </0  équivalant  à  une  dérivation  ent,  cette  équation  (io3), 
divisée  parm,  deviendra  presque  immédiatement,  en  y  mettant  en  évidence, 
comme  dans  (107),  la  valeur  moyenne,  sur  toute  l'étendue  de  la  section, 
d'un  coefficient  variable,  valeur  qui  s'exprime  simplement  au  moyen  du 


(65  ) 
coefficient  b  propre  à  la  section  circulaire  ('  ), 


a  — 


<»4)  ^\^)+[wb+ï{î-*)]-^+>*=°- 


»   Les  conditions  définies  (io/i)  et  (io5)  deviennent  en  même  temps,  vu 
que  (h  (ou  d ^t\2  ■+-  X?  suivant  la  normale  au  contour)  est  idt  à  la  limite 

t  =  i, 


(n5) 


V(o)  =  o,        V(i)  =  o,        j£^j=o(pourï-  =  o), 


I  cl    /1"\  W  ,  v 

H =  o  (pourt  =  i). 


k\/B 


dx 


»  o9.  La  fonction  W  (t)  se  développera  en  une  série  procédant  suivant 
les  puissances  entières  de  t.  Cette  série,  d'après  la  première  condition  (  1 1 5), 
n'aura  pas  de  terme  indépendant  de  t.  D'ailleurs,  la  valeur  (i  i3)  de  A2Œ 
devant  rester  finie  au  centre,  W  contiendra  le  facteur  t  ;  et  le  terme  du  pre- 
mier degré  manquera  dans  W.  Ceux  des  troisième  et  quatrième  degrés 
manqueront  également;  car  un  terme  en  t,3  ne  vérifierait  pas  la  troisième 
condition  (n5),  et  un  terme  en  t4,  porté  dans  (u4)>  )r  en  donnerait  un 
ent"2,  incapable  de  se  réduire  avec  aucun  autre.  Les  deux  premiers  termes 
seront  donc,  l'un,  ent2,  l'autre,  ent5.  Après  quoi,  viendront  t8,  t11,  t,14, . . .  ; 
car  toute  expression  de  la  forme  Mta,  substituée  dans  (n4)>  y  donne  trois 
termes,  respectivement  affectés  de  ta_G,  t*-3,  ta,  essentiellement  différents 
de  zéro  tous  les  trois  pour  a  >  5  et  dont  les  deux  premiers  ne  pourront  se 
réduire  qu'avec  d'autres  issus  de  même  des  deux  termes  de  W  où  a  était 
moindre  soit  de  3,  soit  de  6  unités.  Ainsi  la  différence  des  divers  expo- 
sants a  est  toujours  un  multiple  de  3;  et  si,  W  n'étant  déterminé  qu'à  un 
facteur  constant  près,  l'on  prend  —  i  pour  second  coefficient,  il  viendra 

(i  16)  W(i)  =  At2  -  t'°  4-  Ct8  -  Dr*  '  -h  Et' 4  -  . . .. 

(l)  Voir  la  première  des  formules  (43)  de  mon  Étude  de  l'année  dernière  (p.  34,  ou 
Comptes  rendus,  t.  CXXIII,  p.  77).  Cette  formule  donne 

1  2  1 


k  VB        3        k  sjb 

Dans  le  cas  limite  où  y/B  =00,  qui  sera  considéré  ci-après  en  note,  k\Jb  atteindra  sa 
plus  forte  valeur  f . 

B.  9 


(66) 

»  Une  loi  de  récurrence  assez  simple,  fournie  par  la  vérification  iden- 
tique de  r équation  indéfinie  (i  i4),  permettra  d'évaluer  chaque  coefficient, 
à  partir  de  C,  en  fonction  linéaire  des  deux  coefficients  précédents  multi- 
pliés par  m.  Puis  la  deuxième  condition  (ii5)  déterminera  A,  et  la  qua- 
trième (ii 5)  deviendra  enfin  l'équation  en  m  ('). 


(')  Sur  l'établissement  du  régime  uniforme  dans  un  tube  fin  à  section  circu- 
laire, et  sur  la  convergence  des  séries  rencontrées  dans  la  partie  actuelle  de  ce 
travail.  —  Pareillement  à  ce  qui  arrivait  pour  un  tuyau  à  section  rectangulaire  large 
et  à  paroi  très  rugueuse,  la  loi  de  récurrence  et  la  dernière  relation  (n5),  équation 
en  m,  se  simplifient  quand  B  croît  indéfiniment,  cas  où  k  \/b  tend  vers  f.  Mais,  alors, 
malgré  l'immobilisation  relative,  par  le  frottement  extérieur,  de  la  couche  fluide  con- 
tiguë  à  la  paroi,  la  fonction  W  obtenue  ne  s'applique  pas  au  problème  de  l'établisse- 
ment du  régime  uniforme  le  long  d'un  tube  fin  à  section  circulaire,  parce  que  le  coef- 

og   ,—  R 
ficient  s  de  frottement  intérieur,  de  la  forme  ^-vB—  u0F(t),  ne  pourrait  être  rendu 

A  2 

constant  qu'en  prenant  F  =  const.  ;  ce  qu'on  n'a  pas  fait. 

Pour  trouver  l'expression  de  ra  qui  convient  au  cas  d'un  tel  tube  fin,  il  faut  donc 

poser,  par  exemple,  F  =  i.  Dans  celte  hypothèse,  le  système  (9)  donne  Fj  =  i  —  t2, 

0)L F,  =  |;  et  il  résulte  de  la  formule  (5i),  prise  avec  y/B0  infini,  cp  —  2(1  —  t2).  L'on 

a  aussi,  d'après  la  formule  (3i)  du  premier  Mémoire   (p.   27),  A:\/Bm0=2U  :  d'où 

z  —  r—  UR;  et  la  valeur  de  —  à  porter  dans  les  exponentielles  de  (106)  est       "    -• 

K  k  û  g  U  r\ 

De  plus,  dans  l'équation  (io3),  où  les  deux  derniers  termes  du  quadrinome  entre  pa- 
renthèses donnent  toujours  pour  somme  —  |FjO',  soit,  actuellement,  \tQ',  la  dérivée 


en  t  de  celte  somme  est 


d.t  il' 


,  c^est-à-dire,  identiquement,  2iA2i2.  Dès  lors,  ii  ne 

figurant  dans  (107)  que  par  son  paramètre  A2>  c'est  celui-ci  qu'il  y  a  lieu  de  prendre 
pour  fonction  de  t  à  déterminer.  Appelons-le  n,  ou,  autrement  dit,  réduisons  la  for- 
mule (106)  à 


et  l'équation  indéfinie  (io3)  devient 


d  V  1     d  f   dW\ 
di,  L  4  *  dt\     dt  J 

1t\-t(' 

2  dt  [_  2 1  dv  \ 


4-  m  (  1  —  t2  )  n 


2  m  t  II  —  o , 


d\\  Y 


m  (1 


t2)-T-     —     O. 

dt 


du 


Prenons-y  pour  fonction  inconnue  la  quantité  t—  =  *,  et  pour  variable  indépen- 
dante le  carré  t2  =  s,  aire  relative,  comparativement  à  la  section  entière  j  =  -R2,  de 
la  courbe  aitt  d'égale  vitesse  constituant  le  lieu  des  points  considérés.   Alors  cette 


(«7  ) 


§  XXIX.  —  Simplification  des  intégrales  quand  la  paroi  est  polie. 

»  60.  Mais  bornons-nous,  comme  nous  l'avons  fait  pour  la  section  rec- 
tangulaire larçe,  au  cas  d'une  paroi  assez  polie,  ou  plutôt  d'une  valeur 
de  Passez  petite,  pour  que,  dans  (i  il\),  le  coefficient  du  second  terme  soit 
réductible  à  sa  valeur  moyenne  inverse  de  y/6,  alors  grande  comparati- 
vement à  son  écart  (variable  entre  —  et  —  -  j  d'avec  cette  moyenne.  Le 


équation  différentielle,  s'abaissant  au  second  ordre,  acquiert,  multiplié  part,  sa  forme 
la  plus  réduite, 

(a)  s— —  -t-/n(i  —  s)<P  =  o. 

Il  s'y  adjoint,  en  vertu  de  la  seconde  relation  (104),  la  condition  *  =  o  pour  s  z=  o, 
c'est-à-dire  <ï>(o)  — o,  et  celle-ci  détermine,  à  un  facteur  constant  près  dont  on  peut 
disposer  de  manière  que  *'(o)  =  1,  l'intégrale  particulière  *  à  choisir.  Il  vient  d'ail- 
leurs, en  tenant  compte  de  (100)  où  yB0  a  crû  indéfiniment, 

f\dt  1    /*'  ^  ds 

J  x,  IIS 

et,  puis,  vu  la  formule  11  =  —  —  l  t  — —  |  définissant  n,  vu  aussi  la  première  condi- 
1  i\-c  dt\    di  J  * 

tion  (104), 

do.        rx  rs 

t.-j-=4        Utdt  =  2        Hds. 
Alors  enfin  la  dernière  relation  définie  (io4)  devient  l'équation  en  m  : 

Uds  =  o. 


v  n 


Mais  toutes  ces  formules  ont  eu  beau  se  simplifier  ainsi   à  un   degré   inespéré,  le 
développement  de  <£,  savoir 

(J3)  *  =r  S  —  A2,Ç-+  A353 —  A4S*-f-  A355  — 

n'en  appartient  pas  moins  au  type  complexe  des  séries  qui  s'étaient  déjà  présentées 
dans  nos  trois  autres  problèmes  sur  l'établissement  du  régime  uniforme  à  l'entrée  de 
tuyaux  ou  de  tubes,  soit  rectangulaires  larges,  soit  circulaires  (p.  60  et  66)  :  la  loi  de 
récurrence  résultant  de  l'équation  (a)  y  rattache  chacun  des  coefficients  (pris  en 
valeur  absolue)  A2,  A3,  A4,  .  .  .  aux  deux  coefficients  précédents  o,  1,  A2,  A3,  .  .  .,  et 


(68) 

troisième  terme,  iW,  pourra  encore  être  supprimé,  comme  étant,  pour  les 
mêmes  raisons  que  dans  la  section  rectangulaire  large,  tout  au  plus  compa- 
rable à  la  partie  ainsi  négligée  du  second  terme;  et  si  l'on  introduit,  pour 
abréger,  une  fonction  u  et  une  constante  K  définies  par  les  relations 

d  fW'\  ^m 


dt\x,J  k\Jb 

l'équation  indéBnie  (n4)>  s'abaissant  au  second  ordre,  deviendra 

non  à  un  seul  comme  dans  les  séries  les  plus  usuelles.  En  effet,  l'on  trouve  ici 
m   .  m 


(Y) 


A2= (i+o),  A3=— 5(A,  +  l), 

1.2  2.0 

Ak=^j (A,  +  A,),       A3  =   ™  (A*  +  A3), 


La  convergence  de  ces  sortes  de  séries  peut  se  démontrer  comme  il  suit  : 

Prenons  pour  exemple  (P).  Soit  {M  le  plus  grand  de  deux  coefficients  consécutifs 

m  , 

Ap_!,  Ap,  supposés  assez  éloignés  dans  la  série  pour  que  le  rapport  — et,  a  plus 

forte  raison,  les  rapports  suivants — — r>   -, r-, çt»    •  •  •>  n  excédent 

F1  (/?  +  l)(/>+2)        (/>+2)(/>  +  3) 

pas  une  fraction  donnée  £X,  inférieure  à  {.  On   aura,   d'après  la  loi  de  récurrence, 
d'abord 

Ap+i<  -(Ap-h  Ap_,)  <*^  (d'où,  a  fortiori,  Ap+.,<  -Mj, 


AP+2<  -(Ap+,  +  Ap)<  A  — , 


et,  ensuite, 


Ap+s  <  -  (  Ap+2  +  Ap+,  )<  X2  —  (d'où  Ap+3  <  X  —  j 


Ap+4  <  -  (  Ap+,  -1-  Ap+2  )<  X2  — , 

Ainsi,  les  coefficients,  groupés  deux  par  deux  à  partir  de  Ap_,,  admettront,  dans  les 

.-',,..  .  •  M     ,  M     ..M     ...M  n  .        , 

groupes  successits,  les  limites  supérieures  — ■>    X — >   kL  —  >   h* — ■>    ••••    uonc,  dans  la 

série  (  p),  où  s  varie  de  o  à  i ,  la  somme  absolue  des  termes  croîtra,  à  partir  du  p  —  ilème, 
plus  lentement  que  ne  fait  celle  d'un  nombre  de  termes  moitié  moindre,  dans  la  pro- 
gression géométrique  précédente  doublée,  savoir  M  f-  MX  -h  MX2  +  MX3  -+- .  .  . , 
comptée  à  partir  du  premier  M;  et,  à  plus  forte  raison,  la  série  ((3),  composée  de 
termes  moindres  à  signes  alternant,  sera-t-elle  convergente. 


(69) 
»   Portons-y  l'expression  de  \j.  résultant  de  (i  17)  et  (116),  savoir 

(119)  |t=  -3.5t2  +  6.8Ct5-...: 

l'équation  (118)  déterminera  immédiatement  chacun  des  coefficients 
C,  D,  E,  . . .  en  fonction  du  coefficient  précédent  j,  C,  D,  ...  ;  et  la  for- 
mule (1 16)  sera 

Kt8  KH"  K3tu  K4t17 

v  /  6.8        6.8.9.U        6.8.9.1 1. 12.14       6.8.9.11.12.14-15.17 

)>   Portons-y  la  valeur  de  A  résultant  alors  de  la  seconde  condition  (1 15), 
savoir 

(iai;     A  — 1       6.8  +  6.8.9. u  "     6. 8. 9. 11. 12. 14  +  6. 8. 9. 11. 12.14. i5. 17 

et  la  dernière  condition  (1  id),  multipliée  par  —  ik\[h,  puis  divisée  par  3.5, 
deviendra  l'équation  en  K  (ou  en  ni)  : 

K  K2  K3  K4 


3.5        3.5.6.8        3.5.6.8.9.11        3. 5. 6. 8. 9. n. 12. 14 
(122)    / 

|  _f.WË/I_^_h_^ «K» 

(  5      V  6.8^6.8.9.11         6. 8. 9. 11. 12. 14  / 

»  Conformément  à  ce  qu'on  avait  prévu,  elle  n'a  pas  de  racine  néga- 
tive; car  tous  les  termes  de  son  premier  membre  sont  essentiellement 
positifs  pour  les  valeurs  négatives  de  K. 


§  XXX.  —  Longueur  alors  nécessaire  pour  l'établissement  du  régime  uniforme 
et  autres  particularités  intéressantes. 

w  61.  Si  le  coefficient  \k\fË>  était  assez  petit  pour  qu'on  pût  négliger 
dans  (122)  la  série  où  il  figure,  quelques  tâtonnements  donneraient,  comme 
valeur  de  la  plus  petite  racine,  K  =  25,64-  Mais,  cette  valeur  rendant  la 
série  dont  il  s'agit  positive  (égale  à  0,272),  il  faudra  prendre  K  un  peu 
plus  grand.  Pour  la  valeur  k\B  =o,85i,  qui  convient  à  une  paroi  en 
ciment  lissé,  quelques  essais  donnent  assez  exactement  K  =  3o;  et  il  vient 
ensuite,  d'après  (121),  A  =  o,5342.  On  aura  donc  pour  le  paramètre  m,  vu 
(1 17),  7,5£  y/6,  et,  dans  l'exponentielle  correspondante  de  la  formule  (106) 

de  zs,  l'exposant  sera,  en  valeur  absolue,  — ~—  0  ou,  très  sensiblement, 


(  7°  ) 

o,o3q  —  (vu  que  \/b  =  °>OI29)-  Comme  les  valeurs  initiales  ct,  de  l'écart  cj, 

dans  le  mode  de  distribution  des  vitesses,  ne  dépasseront  guère  encore  un 
dixième,  ou  tout  au  plus  un  dixième  et  demi,  ces  écarts  ts  se  réduiront  à  des 
quantités  de  l'ordre  de  0,01,  et  seront  insensibles,  quand  l'exponentielle 
n'excéder.i  paso,i,  ou  quand  l'exposant  atteindra  2, 3026  en  valeur  abso- 
lue; ce  qui  aura  lieu  pour  x=  5qR,  à  très  peu  près.  Un  parcours  d'environ 
3o  diamètres,  après  V épanouissement  des  filets  fluides  consécutif  à  la  contrac- 
tion de  V entrée,  suffira  donc  pour  établir  le  régime  uniforme. 

»  62.  En  comptant  4  ou  5  diamètres  en  plus  depuis  l'entrée  jusqu'à  la 
section  où  l'épanouissement  est  ainsi  effectué  et  où  le  régime  commence 
à  varier  graduellement,  on  voit  que  l'établissement  du  régime  uniforme 
dans  un  tuyau  de  conduite  à  parois  polies  demandera,  au  maximum,  une 
longueur  de  35  à  4°  f°is  1°  diamètre  du  tuyau.  Conformément  aux  obser- 
vations récentes  de  M.  Bazin,  ce  régime  devait  donc,  dans  ses  expériences, 
exister  après  un  parcours  de  5o  diamètres,  mais  non  après  un  parcours  de 
23  diamètres,  où,  seulement,  l'expression  (106)  de  l'écart  u  était  évidem- 
ment réduite  à  son  terme  principal.  Or,  celui-ci  est,  vu  les  formules  (  1 1 3  ) , 
(120)  de  A2  Œ  et  de  ^'(t-),  et  les  valeurs,  3o,  o,5342,  de  R  et  de  A, 

[  -0,039  ~ 

\  cr  =  —  ce 
(I23)  x  (r,o685-5i3-f-5t6  -2,o833ï9  +  o,4735t'2 

'  —  o,o6y6t15  +  0,0066V8  —  o,ooo5t21  -h  .  .  .). 

»  La  fonction  de  t  entre  parenthèses,  à  laquelle  l'écart  cj  est  propor- 
tionnel dans  chaque  section,  décroît  de  i,o685à  —  0,6028,  quand  x  grandit 
de  zéro  à  r,  c'est-à-dire  quand  on  va  du  centre  au  contour;  elle  est  donc, 
en  valeur  absolue,  plus  grande  sur  l'axe  qu'auprès  de  la  paroi,  dans  Je 
rapport  de  1,773  à  1.  Et  elle  s'annule  pour  t  =  0,659;  ce  T11*  est  aussi 
d'accord  qu'on  pouvait  l'espérer  avec  l'expression  particulière  de  cr, 
observée  par  M.  Bazin  pour  l'abscisse  x  =  5oR  et  constituée  par  les  diffé- 
rences respectives  des  deux  séries  dénombres  rapportées  au  commence- 
ment de  cette  Étude.  On  y  voit,  en  effet,  que  la  valeur  de  %,,  pour  laquelle 
se  produit  l'égalité  des  nombres  des  deux  séries,  est  voisine  de  f  =  0,625, 
légèrement  moindre,  toutefois,  et,  par  conséquent,  un  peu  inférieure  à 
o,65q.  PRr  SIHte»  **i  nL,l  en  moyenne,  ayant  ainsi  le  champ  de  ses  valeurs 
négatives,  qui  constitue  la  région  centrale  des  sections,  sensiblement  ré- 
duit, celui  de  ses  valeurs  positives,  constitué  par  la  région  périphérique, 


(7'   ) 
se  trouve  accru  d'autant;  et  les  valeurs  absolues  de  zj  constatées  vers  le 
contour  sont  encore  plus  faibles  que  les  valeurs  calculées. 

»  La  raison  de  ces  écarts  est  évidemment  dans  l'insuffisante  petitesse 
du  paramètre  k\b,  qui  ne  justifiait  pas  tout  à  fait  les  simplifications  aux- 
quelles nous  avons  soumis  l'équation  différentielle  (n4). 

»  63.  Observons  à  ce  propos  que  la  valeur  de  t  pour  laquelle  gj  s'annule 
et,  par  conséquent,  l'étendue  de  la  région  centrale  où  u>  est  négatif,  grandis- 
sent lorsque  y^B,  \[b  tendent  vers  zéro.  Car,  si  l'on  suppose  \/B,  \b  infini- 
ment petits,  il  vient,  comme  on  a  vu,  R=  23,  64  (d'où  m  =  6,4 1  k  \[b, 
A  =  o,585o);  et  la  formule  (i23)  est  remplacée,  d'après  (106),  (i  i3), 
(121)  et  (120),  par  celle-ci  : 

\  cy  =  —  ce  A    R 

(I24)  \  x  (1,1701  -  5t3+4,2733t6-  [,52i8t9  4-o,2956t,:2 

(  —  o,o36r  t15  -h  o,oo3ot18  —  o,ooo2t'-1  ■+- . ..). 

»  Pour  t  croissant  de  o  à  1 ,  la  fonction  de  t  entre  parenthèses  décroit  de 
1,1701  à  —  0,8161;  de  sorte  qu'en  valeur  absolue  elle  est  seulement,  sur 
l'axe,  ifols,434  (et  non  plus  ifois,773)  sa  valeur  à  la  paroi.  Aussi,  par  com- 
pensation, la  racine  t  qui  l'annule  est-elle  0,6^36,  c'est-à-dire  un  peu 
supérieure  à  0,039  :  ce  qui  rend  la  région  centrale,  où  gj  est  négatif,  égale 
à  la  fraction  o,6~362,  ou  aux  4^4  millièmes,  de  la  section  totale  g,  au  lieu 
de  la  fraction  o,65o/,  ou  des  434  millièmes  (').  » 


(*)  La  plus  grande  partie  de  ce  Mémoire  a  paru  en  mai,  juin  et  juillet  1897  dans  les 
Comptes  rendus  de  l'Académie  des  Sciences  de  Paris  (t.  CXX1V,  p.  1196,  1261, 
1327,  i4'  1,  1^92,  et  t.  CXXV,  p.  6,  69,  1/42  et  2o3). 
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ADDITIOxN  A  LA  NOTE  DES  PAGES  66  A  68. 

Sur  la  manière  d'embrasser  dans  une  même  analyse  les  deux  cas  de  l'écoulement 
le  long  des  tubes  fins  et  de  V écoulement  tourbillonnant. 


En  dehors  du  cas  d'un  régime  uniforme  ou,  tout  au  plus,  d'un  régime  quasi  uniforme 
considéré  seulement,  comme  aux  notes  des  pages  6o  et  66,  dans  ses  circonstances  à  peu 

près  indépendantes  des  légères  variations  de  u0  et  de  -  qui  peuvent  s'y  produire,  les 

deux  hypothèses  F  =  i,  B0=oc  deviennent  insuffisantes  pour  rattacher,  comme  cas 
limite,  les  écoulements  bien  continus  entre  parois  mouillées,  aux  écoulements  tour- 
billonnants; car  la  proportionnalité  admise  du  coefficient  e  de  frottement  intérieur  à 

«0  et  à  -  fait  varier  -.  avec  x  et  t.  Pour  embrasser  dans  une  même  analyse  les  deux 

modes   d'écoulement,  il  convient   alors   d'introduire,  dans   les   expressions  de  £  et  du 

/  <T  Y" 

frottement  extérieur  vei  un  même  nouveau  facteur,  de  la  forme  (  —  I  m",  où  m,  n 
sont  deux  exposants  constants.  Autrement  dit,  l'on  posera 

ce  qui,  d'une  part,  ne  change  pas  l'expression  du  rapport  —  et,  d'autre  part,  entraîne 

les  deux  conditions  £  =  const.,  u0~  o  (caractéristiques  des  mouvements  bien  continus 
entre  parois  mouillées)  quand  on  prend  non  plus  m  =  O,  n  =  o,  mais  bien  m  —  —  i, 
n  —  —  i ,  en  même  temps  que  F  =  i,  ^B0  —  oo,  k  =  ce. 

J'avais,  du  reste,  observé  déjà,  dans  le  Mémoire  de  1878  cité  plus  haut,  p.  20,  et 
publié  au  Journal  de  Mathématiques  pures  et  appliquées  (3e  série,  t.  IV,  p.  344); 
après  y  avoir  donné  directement  la  théorie  des  écoulements  bien  continus  et  graduel- 
lement variés  les  plus  simples,  que  cette  généralisation,  très  naturelle  pour  certains 
modes  d'écoulement  intermédiaires  observés  dans  de  petites  sections  (p.  32  du  Mé- 
moire précédent,  de  1896),  ne  compliquait  pas  notablement  les  formules. 

Et,  en  effet,  dans  les  relations  (1)  à  (i5)  ci-dessus,  elle  n'introduit  pas  d'autre  chan- 

-  .  /  <j  \  1~'" 

gement  que  de  substituer  aux  deux  facteurs  -.■>    u\,  partout  où  ils  figurent,  (  -  1 

et  u\~n .  Rien  n'est  donc  changé,  en  particulier,  aux  systèmes  (9),  (10)  déterminant  F, 
et  F2,  ni,  par  suite,  à  la  formule  de  9,  donnée  par  (5i),  ni  à  l'expression  (i4)  de  DÏLF2, 


(73  ) 


non  plus  qu'à  celle,  4  /  — — — —  j  de  —  dans  le  régime  uniforme.  Mais  l'équation  (5) 


du  mouvement  devient 


(5  bis) 


i  b.cv=(;)'-(.-^), 


i=w(!l    ""(-5)   + 
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Or  la  formule  (i5)  modifiée,  élevée  à  la  puissance  —  (2  -f-  n)  et  traitée  comme  dans 


le  texte,  fournit  une  expression  approchée  de  (  —  )       ,  et  donne  pour  cette  équation 
du  mouvement,  au  lieu  de  (18), 


18  bis)      \  —  b 


b 
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relation  où  les  valeurs  de  Oïlu',  DTl(u-)'  résultent  ensuite  de  (24.).  Puis,  <p  n'ayant 
pas  changé,  rien  n'est  modifié  (p.  27  à  ^2)  à  la  théorie  du  petit  mouvement  trans- 
versal, non  plus  qu'aux  expressions  de  l'accélération  longitudinale  u'. 

Dans  l'hypothèse  n  ■=.  —  1 ,  qu'on  devra  faire  s'il  s'agit  de  mouvements  bien  continus, 
le  dernier  terme  de  (18  bis),  en  DR.  u' ,  disparaît. 

En  même  temps,  cette  équation  (18  bis)  devient  applicable  sans  que  l'écart  m  du 
mode  de  distribution  des  vitesses  soit  petit  devant  tp,  à  la  seule  condition  d'ajouter  au 

second  membre  le  terme ••  Car  l'équation  (i5)  modifiée,  multipliée  par  u0, 

o 

devient  linéaire  en  u0,  comme  l'est  alors  déjà  (5  bis),  et  il  suffit  de  substituer  <p  —  cp0 
à  A  t  /  Ft  (comme  au  bas  de  la  page  5i),  pour  qu'elle  donne 


h=uVî57 


gB0DKf\X 


dK(yu'—  u') 


DK,(u*)' 
2L1 


Uïlu'-h3ÏL(—  mu') 


valeur  de  u0  à  porter  simplement  dans  (5  bis). 

Par  suite,  la  hauteur  motrice  totale  à  dépenser  pour   établir  le  régime   uniforme, 
dans  un  tuyau  ou  tube  évasé,  sera,  au  lieu  de  (92), 


(92  bis) 


U2 

8 

L2 


- — -  4-  ~   f    3ïi(—  wu')dx\ 
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(  74) 

Le  nombre  a  égale  2  quand  le  tube  est  circulaire;  en  soi  te  que,  si  l'on  néglige  dans 
un  premier  calcul  le  petit  produit  ( — m)u',  cette  hauteur  de  charge   dépensée  prend 

U2 

alors   la  valeur   approximative   simple Les   principes    exposés    dans   la    note   des 

S 
pages  66  à  68  permettent  d'ailleurs  d'évaluer  —  m  et  a'  en  série  :  après  quoi,  une  qua- 
drature numérique,  effectuée,  par  exemple,  au  moyen  de  la  formule  de  Thomas  Simpson, 
donne  le  dernier  terme  entre  crochets  de  (92  bis).  C'est  ainsi  qu'ont  été  obtenus  les 
résultats  résumés  à  la  note  de  la  page  20. 

Enfin,  quel  que  soit  l'exposant  /i,  les  théories  des  §§  XXlil  et  suivants,  relatives  à 
la  distribution  des  vitesses  dans  la  première  partie  amont  des  tuyaux  ou  des  tubes, 
subsisteront,  sauf  le  remplacement,  dans  les  équations  (99),  (102)  et  (106),  du  facteur 

g  y  ?<Un  // Y"'"         ...  gU" 

-  par  n  '  "  —  1-1        ,  qui  revient  a 
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ERRATA 


Mettre  en  note  au  bas  de  la  page  8,  avec  signe  de  renvoi  à  la  deuxième  ligne  du 
n°  h,  après  les  mots  axe  hydraulique,  la  phrase  suivante  : 

J'appelle  axe  du  courant,  ou  axe  hydraulique,  comme  on  verra,  du  reste,  à  la 
page  10,  l'axe  même  du  tuyau  contenant  le  fluide,  quand  celui-ci  coule  le  long  d'un 
tuyau  qu'il  remplit,  et,  dans  le  cas  contraire,  d'un  courant  liquide  limité  supérieure- 
ment par  une  atmosphère,  une  ligne  tracée  sur  la  surface  libre  et  suivant  sa  longueur, 
comme,  par  exemple,  à  égale  distance  des  bords,  ou  encore  suivant  l'axe  du  tuyau 
dont  la  section  comprendrait,  avec  celle  du  canal  découvert  considéré,  sa  symétrique 
par  rapport  à  la  coupe  transversale  de  la  surface  libre. 

A  la  page  58,  après  la  formule  (106  bis),  ajouter  : 

»  11  est  clair  que,  à  la  limite  où  le  nombre  n  des  termes  deviendrait  infini  et  où  la 
somme  2 cA2Q  serait  susceptible  d'exprimer  la  fonction  arbitraire  — wh  l'intégrale 
(106  bis)  admettrait  la  valeur  minima  zéro,  que  donne  l'égalité  alors  possible 
lcA2<2=: — T7T;  annulant  chacun  de  ses  éléments  ;  en  sorte  que  la  recherche  de  cette 
valeur  minima  fournirait  les  vrais  coefficients  c.  Ainsi,  ces  coefficients  c,  si  on  les 
détermine,  môme  quand  /*  est  fini,  par  la  condition  de  rendre  minimum  l'intégrale 
(106  bis),  doivent  bien  tendre  vers  leurs  vraies  valeurs  limites  à  mesure  que  leur 
nombre  n  s'accroît. 
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